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Preface 



This volume contains a selection of invited papers, often in expanded form, that 
are based on presentations made at the 6th Conference on Clifford Algebras and 
their Applications in Mathematical Physics, May 20-25, 2002, in Cookeville, 
Tennessee [1]. The organizers of the conference, Rafal Ablamowicz (Tennessee 
Technological University), and John Ryan (University of Arkansas) grouped all 
conference presentations into five sessions in Clifford analysis, geometry, mathe- 
matical structures, physics, and applications in engineering. Thus it seemed natu- 
ral to organize the book into five parts under the same titles. 

The 6th Conference on Clifford Algebras continued a 16-year sequence of in- 
ternational conferences devoted to the mathematical aspects of Clifford algebras 
and their varied applications in mathematical physics, and, more recently, in cy- 
bernetics, robotics, image processing and engineering. Previous meetings took 
place at the University of Kent, Canterbury, U.K., 1985; University of Montpel- 
lier, Montpellier, France, 1989; University of Gent, Gent, Belgium, 1993; Univer- 
sity of Aachen, Germany, 1996; and Ixtapa, Mexico, 1999. Three edited volumes 
appeared after the Ixtapa conference, published by Birkhauser, Boston [2-4]. 

Chapter 1, in this volume, is devoted exclusively to topics from Clifford anal- 
ysis and range from the Morera problem, inverse scattering associated with the 
Schrodinger equation, through discrete Stokes equations in the plane, a symmetric 
functional calculus, Poincare series, to differential operators in Lipschitz domains, 
Paley-Wiener theorems and Shannon sampling, Bergman projections, and quater- 
nionic calculus for a class of boundary value problems. Among geometry topics, 
not so visibly present at previous conferences that gave rise to Chapter 2, are 
spin structures and Clifford bundles, eigenvalue problems for Dirac and Rarita- 
Schwinger operators, differential forms on conformal manifolds, connection and 
torsion, Bochner identities on Riemannian manifolds, and noncommutative ge- 
ometry in physics. Chapter 3 is devoted to mathematical structures such as Grass- 
mann algebras, Lie superalgebras, Grassmann-Hopf algebras, symplectic Clifford 
algebras and graded central algebras. Applications in physics, collected in Chap- 
ter 4, cover a wide range of topics from classical mechanics to general relativity, 
twistor and octonionic methods, electromagnetism and gravity, elementary parti- 
cle physics, noncommutative physics, Dirac’s equation, quantum spheres, and the 
Standard Model. Chapter 5 includes papers on Clifford geometric algebras and 
designing co-processors in computer engineering, applications in the description 
of an image space using Cayley-Klein geometry, and pose estimation. Below the 
reader will find brief introductions to all chapters in this volume. 

While the papers collected in this volume require that the reader possess a solid 
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knowledge of appropriate background material, as they lead to the most current 
research topics, the fundamentals of this background were presented prior to the 
6th Conference in the form of six lectures. The lectures were delivered by promi- 
nent specialists in the field (including Pertti Lounesto to whom this volume is 
dedicated). They were aimed at graduate students and newcomers to the field. 
These lectures will appear in a separate volume from Birkhauser in 2003 [6]. 

Papers from regular sessions at the 6th Conference included here were selected 
by the conference session organizers who oversaw the refereeing process: 

o Clifford Analysis: Marius Mitrea (University of Missouri-Columbia, Co- 
lumbia, Missouri), and Mircea Martin (Baker University, Baldwin City, 
Kansas); 

o Geometry: Tom Branson (University of Iowa, Iowa City, Iowa), and Ugo 
Bruzzo (International School for Advanced Studies, Trieste, Italy); 

o Mathematical Structures: Ludwik Dabrowski (SISSA, Trieste, Italy), and 
Bertfried Fauser (Universitat Konstanz, Konstanz, Germany); 

o Physics: William Baylis (University of Windsor, Windsor, Canada), and 
Giovanni Landi (Universita di Trieste, Trieste, Italy); 

o Applications: Jon Selig (South Bank University, London, England), and 
Gerald Sommer (Christian- Albrechts-Universitat Kiel, Kiel, Germany). 

Papers presented as plenary talks were selected by the Editor. The main aim of 
this volume, which guided the selection process, was to collect the best papers 
in each of the five representative areas ranging from pure mathematical theory 
to a wide array of applications in physics and computer engineering. Following 
the spirit of the conference, this volume unlike previously published volumes, in- 
cludes a separate chapter on geometry. 

Dedication 

I am dedicating this book to the memory of Pertti Lounesto, friend and colleague, 
whose plenary lecture at the 6th Conference proved to be his last major contribu- 
tion to the field of Clifford algebras - the field he loved so much. 



PART I. CLIFFORD ANALYSIS 

Berenstein, Chang and Eby: In this paper the authors discuss the Morera prob- 
lem for functions / taking values in a Clifford algebra Cl n . As a replacement 
for the Morera theorem from the complex plane C, which gives necessary and 
sufficient conditions for a continuous function / in an open subset fl of C to be 
holomorphic in they formulate integral conditions whose vanishing on a Lip- 
schitz boundary of every bounded domain secures that / be left regular. They 
establish a connection, similar to the one from the Euclidean space, between the 
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Pompeiu and Morera problems in Clifford analysis. Finally, the Morera theorem 
under Mobius transformations of a disk and the moment problem in the context 
of Clifford analysis is considered. Then, the authors discuss functions defined on 
the Heisenberg group H n and investigate conditions under which these functions 
are CR-functions. In particular, they study moment and conjugate moment condi- 
tions that, if satisfied by functions in L 2 (H n ), L p ( H n ), and eventually L°°(H n ), 
guarantee that these functions are CR-functions or conjugate CR-functions, re- 
spectively. 

Bernstein: The author applies Clifford analysis methods to a study of the inverse 
scattering problem for an n-dimensional Schrodinger- type equation. In particular, 
she relies on a Borel-Pompeiu formula in C n and a representation of holomor- 
phic functions from complex Clifford analysis. This Clifford algebra setting is 
particularly useful in explaining the so-called “compatibility conditions for in- 
verse scattering”. A method of reconstructing the potential from a scattering data 
by quadratures with the help of the Borel-Pompeiu formula is provided. It is ap- 
plied to a regularized Schrodinger equation in C n and to time-dependent and the 
time-independent Schrodinger equation in R n . In fact, for a data T fulfilling a 
certain operator equation determined by an operator V defined by the author in 
analogy to the 3 operator, the compatibility conditions, and decaying properly as 
|fc| — > oo, the scattering potential is computed by taking the inverse Fourier trans- 
form of the scalar part of the Borel-Pompeiu formula that involves T. 

Giirlebeck and Hommel: The authors apply the theory of discrete analytic func- 
tions to the solution of the Dirichlet problem for the Stokes and Navier-Stokes 
equations. Finite difference operators are used to approximate the Cauchy-Rie- 
mann operator. A discrete T-operator and a discrete Cauchy integral are studied. 
In particular, authors discuss that the T-operator is a right-inverse of the discrete 
Cauchy-Riemann operator. A discrete Borel-Pompeiu formula is established and 
discrete Bergman projections are defined via a suitable orthogonal decomposition 
of the underlying / 2 -space. Adapted finite difference schemes are used to approx- 
imate solutions to the boundary value problems. 

Jefferies: The author discusses a symmetric functional calculus / /(A) for 

bounded monogenic functions / defined in sectors in M n+1 that take values in the 
Clifford algebra C£{ C n ). Here, A = (Ai, . . . , A n ) is a finite sequence of closed 
linear operators, whose real combinations have spectra contained in a fixed sector 
of C. If / has decay at zero and infinity, f(A) is defined even if A \ , . . . , A n do not 
commute. In the case when n — 1 , the author obtains the usual Riesz-Dunford 
functional calculus for a single operator. In particular, if A — (Ai, . . . , A n ) com- 
mute and satisfy certain square functions estimates, the author uses a recently 
proven bijection between bounded monogenic functions defined on a sector in 
M n+1 and bounded holomorphic functions defined on a corresponding sector in 
C n to obtain a functional calculus / h-* /(A) for bounded holomorphic func- 
tions / defined in a sector in C n . The method is applicable to the Dirac operator 
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on a Lipschitz surface E in R n+1 and leads to another proof of the boundedness 
of the Cauchy integral operator on E. 

Krausshar: This paper is a continuation of a series of earlier papers by the same 
author in which began developing a theory of monogenic and polymonogenic 
hypercomplex-valued modular forms using methods of Clifford analysis. Here, 
the author provides nontrivial examples of such forms related to discrete sub- 
groups of the Vahlen group. This is accomplished by constructing nontrivial mon- 
ogenic and polymonogenic Eisenstein and Poincare series. 

Marmolejo-OIea and Mitrea: The authors study a large class of Dirac-like op- 
erators D : C°°(M,£) — > C°°(M J !F) where £, T — > M are two smooth, Her- 
mitian vector bundles on a smooth, compact, boundaryless, orientable, Rieman- 
nian manifold M. They provide a variety of examples of such operators, both 
for real and complex manifolds M. The main goal is to show that there is a 
rather rich function theory associated with these operators, derivable by means 
of real-variable techniques. The starting point is a generalization of the classical 
Pompeiu formula for some suitably defined Cauchy type integral operators. Later, 
they study the interior regularity and boundary behavior of these Cauchy type op- 
erators, and prove results related to Hardy spaces associated with such operators 
in Lipschitz domains. Then they specialize their results by discussing separately 
Clifford monogenic functions and octonionic analysis as manifestations of a gen- 
eral, unifying function theory. Finally, they observe that many of the versions of 
the classical Bochner-Martinelli-Koppelman formula currently scattered in the 
literature are, in fact, particular cases of their generalized Pompeiu formula. 

Qian: The author discusses Paley-Wiener theorems for entire functions in one 
and several complex variables, and their generalization to left-monogenic func- 
tions on Ry = {x = xo + x | xq € R, x E R n }. The R n version uses Alan 
McIntosh’s generalized exponential function e(x,£). The author remarks that a 
proof of the classical theorem in C can be adapted to a proof of the generalized 
version due to the availability of Laurent series, Fourier multipliers and convolu- 
tion integrals of Clifford monogenic functions. An analysis of both proofs follows. 
As an application, the Clifford monogenic sine function is studied. A theorem is 
stated giving an extension of Shannon sampling to functions that are left-Clifford 
entire with exponential bounds of the Paley-Wiener type, and whose restriction 
to R n belongs to L 2 (R n ). A complete proof of the theorem is given. 

Ren and Malonek: The authors prove a theorem that gives conditions when the 
Bergman projection operator P a from the Banach space L P (B, Cl o, n > dV a ) onto 
a subspace of monogenic functions A p ( B, Cl o, n , dV a ) is continuous. The proof 
follows a study of the Clifford algebra valued Bergman kernel. 

Sprossig: Methods of quaternionic operator calculus developed earlier by the au- 
thor in cooperation with K. Giirlebeek for the solution of time- independent bound- 
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ary value problems are extended here to a study of Galpern-Sobolev equations 
with a variable dispersive term. The time derivative is approximated by forward 
finite differences. The author discusses stability and convergence of the time- 
discretization method. 



PART II. GEOMETRY 

Bartocci and Jardim: The authors study a generalization of the Nahm trans- 
form for instantons on minimal resolutions of R 4 /G where G is a finite sub- 
group of SU(2). Their approach allows them to define the Nahm transform on 
hyper-Kahler ALE 4-manifolds. Unlike the usual Nahm transform that associates 
instantons on T n x R 4-n , n = 1, 2, 3, 4 to (possibly singular) solutions of the 
dimensionally reduced anti-self-duality equations on T n , the torus dual to T n , 
the Nahm transformed instanton connection has no singularities. Under certain 
conditions described by the authors this new Nahm transform is invertible. Dif- 
ferential geometric properties of this transform are discussed. 

Grantcharov: The author studies the hyper-Kahler geometry with torsion. The 
geometry of a connection with totally skew-symmetric torsion and holonomy in 
Sp(n) is referred to as an HKT-geometry. He finds that a hyper-Hermitian man- 
ifold admits an HKT-connection if and only if for each complex structure there 
is a holomorphic (0, 2)-form. This allows him to study the potential theory for 
the HKT-geometry. In particular, he constructs a family of inhomogeneous HKT- 
structures on compact manifolds including S 1 x S 4n+3 . Then he discusses the re- 
duction theory for HKT- structures. He eventually shows local existence of HKT 
metrics on any hypercomplex manifold and points out that global existence in 
general does not hold. 

Homma: The author studies the enveloping algebra and Casimir elements of 
so(n). He gives relations among them corresponding to Bochner identities. He 
then defines the principal symbols of gradients called the Clifford homomorphisms 
and relates them to the enveloping algebra through conformal weights. He calcu- 
lates the eigenvalues of Casimir elements. In the last section, Homma gives all 
the Bochner identities for the gradients explicitly. Examples of geometric first- 
order differential operators which can be realized as gradients on Riemannian 
manifolds include the exterior derivative, its adjoint, and the conformal Killing 
operator and on spin manifolds, the Dirac operator, the twistor operator and the 
Rarita-Schwinger operator. 

Hong: The author gives explicit formulas for the eigenvalues of Dirac and Rarita- 
Schwinger operators on M = S 1 x S n ~ l where S 1 is endowed with a Lorentzian 
or Riemannian metric and a nontrivial spin structure, while S' 71-1 (n even) carries 
Riemannian metric and a standard spin structure. 
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Ugalde: Author presents partial results on differential forms associated to even- 
dimensional compact conformal manifolds. He associates to any pseudodifferen- 
tial operator S of order 0 and acting on sections of a bundle B on a compact 
manifold without a boundary M, a differential form h) of order n acting 

on C°°(M) x C°°(M). The form is uniquely given as the Wodzicki 1 -density 
of the product of commutators and can be expressed in terms of the total sym- 
bol of S. The author discusses a family of differential forms 9 ) of order 6 
with /, h £ C°°(M) on a 6-dimensional, conformal, oriented, and compact man- 
ifold M without boundary. He shows that the integral f M /o^6(/i, / 2 ) defines 
a Hochschild 2-cocycle over the algebra Conformally flat and non-flat 

cases are considered. 

Varilly: The author discusses connections between noncommutative geometry 
and quantum physics. As examples that have influenced development of mathe- 
matical theory, he brings the Standard Model (SM) of fundamental interactions, 
noncommutative quantum field theory and Feynman graph combinatorics. The 
paper begins with a review of reconstruction of the SM due to Connes and Lott, 
and mentions the Mainz-Marseille approach to the SM. Then the author dis- 
cusses spectral triples and their role in constructing spin manifolds. He reviews 
the spectral action principle introduced in order to include gravity in the early SM 
model and then moves on to discuss noncommutative field theory. Mathematics 
and physics are closely linked: “Seiberg-Witten map” and gauge theory, Moyal 
product and NC theory, phase space approach to quantum mechanics, noncommu- 
tative spaces and quantum gravity, Hopf algebras and QFT calculations. Extensive 
references are included. 



PART III. MATHEMATICAL STRUCTURES 

Brini, Regonati and Teolis: The authors consider supersymmetric letterplace al- 
gebras Super[C\V] over a pair of Z 2 -graded alphabets C and V of letters and 
places. These algebras can be viewed as bimodules under the action of general 
linear Lie superalgebras pl(C) and pl(V) (and their enveloping algebras) associ- 
ated to the alphabets C and V. Following A. Capelli’s idea, they do so by embed- 
ding Super[C\P) into a larger “virtual letterplace algebra”, obtained by adjoining 
new “virtual symbols” to the former. They describe how “virtual operators” may 
be used in place of operators from the universal enveloping algebras U(pl(C)) 
and U ( pl(V )). As an example of this process they treat generalized Capelli oper- 
ators as products of virtual polarizations. They extend these ideas to the so-called 
“symmetrized bitableaux” that give a basis for Super[C\P ], and “Young-Capelli 
symmetrizers” that act on these bitableaux. They give several complete decompo- 
sition theorems for these modules and algebras. As an application, they consider 
Berele-Regev representation theory of the symmetric group S n over Z 2 -graded 
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tensor spaces T n [V 0 ® Vi] built from Z 2 -graded finite dimensional vector spaces. 

Eastwood: Various filtered and graded algebras arise, like the Clifford algebra, 
as quotients of the tensor algebra over a finite-dimensional real or complex vector 
space V, the latter possibly equipped with a nondegenerate symmetric or skew 
bilinear form. Examples presented include the exterior (Grassmann) algebra, the 
symmetric algebra, the symplectic Clifford algebra (also known as Weyl algebra), 
the enveloping algebra of a (real or complex) Lie algebra, and Laplace algebra that 
arises from symmetries of the Laplacian. These algebras can be recovered via a 
general construction that involves the tensor algebra over an irreducible represen- 
tation W of a semisimple Lie algebra g and suitably chosen homomorphisms of 
g-modules: $ : W -> R or C and ^ : W <g> W — > W. In general, three 
classes of algebras can be obtained that are determined by the symmetric product, 
skew product or Cartan product. The author discusses three conjectures about the 
Cartan and tensor products of irreducible representations of g and shows how they 
apply to the general construction and the specific cases. In an appendix he gives 
an outline of his talk “Symmetry and Differential Invariants” that was presented 
at the 6th Conference and that gave a motivation for the presented study. 

Fauser: A modern categorical approach is necessary to properly formulate grade- 
free product formulas for almost all algebra and coalgebra products in Grassmann, 
Grassmann-Cayley, and Clifford algebras. That is, formulas that apply to general 
multivector polynomials instead of only generators or homogeneous elements. 
General Rota-Stein’s cliffordization process — a categorically formulated gener- 
alization of Chevalley deformation — yields a combinatorial, instead of Cheval- 
ley’s recursive, computationally efficient formula for the most important of these 
products — the Clifford product. Starting from a Grassmann-Hopf algebra, the au- 
thor derives by categorial duality a self-dual Grassmann-Cayley double algebra 
with formulas for meet and join, comeet and cojoin, left/right contraction and co- 
contraction, Clifford and co-Clifford products. Chevalley formulas are then shown 
to be a specialization of this general approach while the co-Chevalley formulas are 
new. The method is based on algebraic invariant theory. Implementations of these 
products by mostly combinatorial algorithms prove to be computationally effec- 
tive. 

Hahn: This is a survey of the theory of quadratic forms q over a commutative 
and associative ring R with unity (often a field of characteristics not 2 (because 
there are integral results as well)), the theory of quadratic modules (M, q) (even- 
tually free modules with finite bases), and the theory of Clifford algebras C(M) 
of quadratic modules. A classification theorem of Clifford algebras as central sim- 
ple algebras is given. The author shows that any finite-dimensional central simple 
algebra A with involution over a field R is “Brauer equivalent” to a Clifford al- 
gebra C(M) for some quadratic space M of even rank and discriminant 1. He 
then moves on to discuss analogous results when R is a global field, or an arith- 
metic Dedekind domain. For quadratic forms over a field F of characteristics 
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not 2 , the author reviews the role of the Witt ring W(F) of F. He shows that 
three maps dim, disc, and Cliff are sufficient to classify quadratic spaces when 
F is a local field, while a total signature map sig is needed when F is an alge- 
braic number field. Then, he moves on to discuss the spin groups Spin (M) as the 
central extensions of the commutator subgroup of the orthogonal group O(M). 
Generalized even Clifford algebras are Clifford algebras C(A, cr, /), where A is a 
finite-dimensional central simple algebra over a field F, char F ^ 2 , and (<r, /) is 
a quadratic pair on A. These algebras are defined as factor algebras of the tensor 
algebra T(A). The associated spin groups are central in the classification theory 
of algebraic groups. Finally, in historical retrospective, the author discusses the 
role that Clifford algebras play in Adams’ solution to the existence problem of 
the largest number p(N + 1) — 1 of linearly independent vector fields on the 
sphere S N where p(N + 1) is the Radon-Hurwitz number of iV + 1. 

Helmstetter: In this first paper the author compares (orthogonal) Clifford alge- 
bras and symplectic Clifford algebras, also called “Weyl algebras”. Following 
Chevalley, one can introduce the Clifford product on the exterior algebra /\ (M) of 
a vector space M equipped with a bilinear form (3 : M x M — > F, (here K = M 
or C, and the form (3 is assumed to have a nondegenerate symmetric part </>). 
Analogously, when the dimension dim M = m is even, M can be equipped 
with a bilinear form 7 that has a nondegenerate skewsymmetric part and the 
finite-dimensional exterior algebra may be replaced with the infinite-dimensional 
symmetric algebra S(M). Indeed it is possible to introduce a 7 -dependent Weyl 
product on the vector space S(M) so that S(M) equipped with this new product 
becomes isomorphic to the Weyl algebra of (M, ^), also referred to as symplec- 
tic Clifford algebra. The author then moves on to state the Lipschitz Theorem for 
the orthogonal Clifford algebras. The theorem relates an inner automorphism of 
/\(M, /?), which is generated by an exterior exponential of an element of /\ 2 (M) 
and that leaves M invariant, with an ”induced-by-it” orthogonal transformation of 
M without an eigenvalue of 1 . In the main part of this paper the author explains 
difficulties and misconceptions related to extending the Lipschitz theorem to sym- 
plectic Clifford algebras: the need for a proper enlargement S(M) of S(M) and 
the need to realize that Weyl multiplication of elements of S(M) often fails due 
to divergent infinite sums. In order to construct a “local” symplectic Lipschitz 
group, the author proves a symplectic counterpart of the Lipschitz theorem. Then, 
he presents ideas as to how to introduce a “global” symplectic Lipschitz group. 

Helmstetter: It is well known that Clifford algebras C£(M, Q) over a vector 
space M with a nondegenerate quadratic form Q are examples of graded cen- 
tral simple (G.C.S.) algebras provided with an involution called the “reversion”. 
In general, finite-dimensional G.C.S. algebras A over a field K are isomorphic to 
matrix algebras (with matrices of suitable sizes) over some graded central divi- 
sion algebra D. Thus, these algebras can be classified by the isomorphy classes 
of D which form the Brauer-Wall group BW (K) under a twisted tensor product 
® as the group multiplication. The author studies G.C.S. algebras A = Aq © A\ 
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over a field K of characteristic •=/=■ 2, provided with an involution p such that 
p(Aj) = Aj for j = 0,1. Using a “complex divided trace” of p, which is an 
element of the group Rg of the eighth roots of 1 in C, he obtains a classifica- 
tion of these involutions. Then, the pairs ( A , p) can be classified by suitable pairs 
(6, r) G BW (K) xRg.As an example, he gives a classification of such pairs over 
R and, as a particular case, he recovers the classification of the Clifford algebras 
C£(p, q ). The last section of this paper contains two specializations of a bijective 
mapping A ® Z 9 (A, M) — > M derived from a graded bimodule M over A and 
the graded centralizer Z 9 (A , M) of A in M. 

Marks: The author describes a binary labeling scheme for powers of generators 
(p = 0, 1, . . . , n — 1) in the Clifford algebra C£(p , q) (n — p + q) in terms of 
prescripts, postscripts, subscripts, and superscripts. While the subscripted (super- 
scripted) indices denote covariant (contravariant) generators, prescripted (post- 
scripted) indices denote monotonically decreasing (increasing) order of the gen- 
erators. Thus, a typical (Clifford) basis monomial in a preferred labeling will look 
like m e = = (7n-i) mn_1 ... (7/x) m " ••• (7o) m ° • One advantage 

of this approach is that the powers can be treated as bits of a binary number 
m = • Then, the grade of the monomial is just the sum of the bits, 

while the actions of reversion, grade involution, conjugation, and duality trans- 
formations on m e can be expressed as simple binary operations on the labels. An 
important formula is, of course, Clifford multiplication of two basis monomials 
ie and m e, which is similar to a formula due to P.-E. Hagmark and P. Lounesto 
that uses Walsh functions and is implemented in CLIC AL. 

Schmeikal: Traditionally, in Euclidean spaces R n with a positive definite metric, 
transpositions of base units are carried out by Weyl reflections. Those reflections 
represent the lattice of root spaces A n connected with the classical Lie algebras 
su(n+ 1, C) and the other real forms of sl(n, C). Schmeikal develops that concept 
further by defining and investigating transposition involutions of Clifford algebras 
generated by spaces with indefinite metric. He applies that tool to the real Clif- 
ford algebra Cl^\ of the Minkowski spacetime to derive the normal real forms of 
the algebra su(3, C) and its Lie group. It becomes evident that a definite algebra 
representation can cover different forms of a group. The original quark multi- 
plet found by Gell-Mann and Zweig is reconstructed. The generators of the inner 
symmetries of strong interaction turn out as Lorentz transformations of inhomo- 
geneous differential forms. The six copies of the normal real SU(3)-form act on 
six color spaces which are isomorphic with the fourfold ring of real numbers 4 M. 
The flavor SU(3) is an isotropy (sub)group at the point of a fixed lepton, the 
most prominent element of a color space. Those spaces are the appropriate spinor 
spaces for fermions with baryon number | . Possibly this concept opens new math- 
ematical schemes to treat the internal symmetries of particle physics.This can also 
lead to a new perspective on the physical meaning of these symmetries. 
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PART IV. PHYSICS 

Baylis: The aim of this paper is to use the spinorial formalism of the Clifford 
algebra of “classical physics” to explore connections between the classical and 
quantum theories. In particular, the author uses his paravector representation of 
spacetime in Ct<$ to discuss spin-| systems. One advantage of selecting Cl 3 as 
the main algebra is that reversion (called “reversal” by the author) in Cl$ coin- 
cides with hermitian conjugation since the orthonormal basis vectors (generators 
of CI 3 ) are represented by hermitian matrices (normally there is no such connec- 
tion). This gives many classical expressions their corresponding quantum form. 
The author discusses spin transformations of paravectors, eigenspinors and their 
time evolution in accelerating and rotating systemS, classical de Broglie waves, 
electromagnetic gauge potential, spin and g-factor, magnetic moments, Dirac’s 
equation, and the Stern-Gerlach experiment. 

Bonechi, Ciccoli and Tarlini: This paper discusses noncommutative geometry of 
Podles 2-sphere and its generalization, the 4-sphere, proposed earlier by the au- 
thors. In particular, they analyze local “quantum charts” that can be used to define 
local trivialization for the monopole and instanton vector bundles. 

Budinich: The author discusses “simple” spinors of E. Cartan (renamed as “pure” 
by C. Chevalley) and shows how the Cartan’s equations defining them may be 
identified as the spinor-field equations in ten-dimensional spacetime R 1,9 after 
restriction to R 1,3 . In this way the author obtains equations of motion for fermion 
multiplets where charges and internal symmetry groups originate from the known 
correlation of Clifford algebras with the three complex division algebras: com- 
plex numbers (J7(l) of charges), quaternions (isospin 517(2), and electroweak 
SU (2) l ) and octonions ( SU (3) flavor and color). Furthermore in this framework, 
momentum spaces, bilinear in simple or pure spinorsn, are compact and also sev- 
eral problematic aspects of recent elementary particle theoretical physics might 
obtain simple geometrical explanations. 

Chen, Nester and l\ing: This paper addresses problems associated with a proper 
mathematical formulation of gravitational energy-momentum using spinors. An 
important achievement in this field (described in this paper) is Witten’s spino- 
rial positive energy proof, which can best be understood in terms of the associ- 
ated Hamiltonian density. The latter is a candidate for a “truly physical” energy- 
momentum density for a gravitational field. The authors discuss the role of the 
Hamiltonian boundary terms, the issue of positivity, (quasi)localization, and the 
various roles the spinor field can play. Using a certain spinor-curvature identity 
they recover Witten’s proof and the Dougan-Mason quasilocalization. A simi- 
lar identity involving SU(2) spinors gives an alternate reparametrization of the 
Hamiltonian, allowing the authors to provide an alternate positive energy proof 
and quasilocalization. Yet another identity allows for the inclusion of the spinor 
field into a “quadratic spinor Lagrangian” as a dynamic physical field. Although 
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they find that spinors give nice results for energy-momentum, they note a serious 
limitation of this approach to gravitational energy-momentum based on spinor 
Hamiltonians: it seems that it cannot give all of the other physically conserved 
quantities of asymptotically flat spacetimes. 

Daviau: The author discusses a chiral relativistic wave equation for neutrinos with 
a nonzero mass term. Unlike in the classical complex formalism, where the Dirac 
equation is formulated in the Clifford algebra C¥ 3 of the space, two new vector 
K( i) = (f)(i i (f )^ , 2 ) = and two new bi vector = </>cri</>, 5(2) = </><7 2 </> 

Lorentz-co variant tensorial densities appear that are related to Pauli’s o\ and <j 2 
directions. They are not electric gauge invariant. Here </> G C£$ is a general multi- 
vector that the author eventually expresses in terms of two Weyl spinors £,77. With 
the K(j) and tensors, the author has 36 densities, and not only 16 as in the 
classical approach. Then, the author undertakes a detailed study of a chiral wave 
equation with mass and charge terms in the spacetime algebra setting. He solves 
this equation for the chiral wave for the hydrogen atom using H. Kruger’s method. 

Dray and Manogue: This paper is a continuation of earlier work by these au- 
thors on providing an octonionic description of the massless Dirac equation in 
(9 + 1) Minkowski space. Using 2-component octonionic spinors they previously 
showed that three generations of leptons combine into such a description. The 
three generations live in three distinct quaternionic subalgebras of the octonions 
whose intersection is a complex subalgebra generated by a preferred octonionic 
unit. This paper begins with a review of the earlier work. The authors then extend 
their method to a 3-component formalism used by J. Schray for the superparticle. 
They hope to formulate a description for particles in the language of a special 
exceptional Jordan algebra of 3 x 3 octonionic Hermitian matrices (here called 
“Jordan matrices”). These matrices satisfy the usual characteristic equation of or- 
der 3. Of interest are those Jordan matrices V that under the Jordan product “ 0 “ 
(symmetric commutator) are primitive idempotents of trace 1 (Moufang plane). 
The authors exploit the fact that each such matrix can be dyadically factored as 
V = vv t where v is a 3-component octonionic column vector with quaternionic 
components. They then consider the Jordan eigenvalue problem A o V — AV for 
Jordan matrices A and V. 

Lasenby, Doran and Arcaute: The authors apply spacetime algebra (STA) for- 
malism to problems in electromagnetism, gravitation and multiparticle quantum 
systems. As examples of applications they discuss solutions to full Maxwell equa- 
tions in the presence of a conductor, a spectrum of bound states created by a 
black hole, conformal models of space and spacetime, and multiparticle systems. 
They give a detailed study of a minimally-coupled relativistic wave equation for 
fermions in a spherically-symmetric black hole. Since this equation reduces to a 
free-field equation with a single interaction term Hj that lacks Hermiticity, ener- 
gies of the gravitational bound states contain an imaginary component that causes 
states to decay. They display sample energy spectra of Si and 2 Pa states. Their 
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goal is to include multiparticle effects in the quantum physics of the black hole. 
Then, they move on to conformal models of Euclidean geometry and their appli- 
cations in computer graphics. In particular, they discuss a stereographic represen- 
tation of R n versus a projective one. The latter approach is especially conveniently 
implemented via blades in a geometric algebra. Conformal geometry enters with 
rotors that, initially, give rotations and translations (Euclidean transformations) 
about the origin in the Euclidean spaces, and then can be used for rotations about 
any point. The last part of the paper is devoted to applications of a geometric 
algebra of the 4n-dimensional relativistic configuration space (called the “multi- 
particle spacetime algebra” or MTS A for short) in entanglement and multiparticle 
wave equations. At the end, the authors return to the conformal geometry in the 
context of R. Penrose’s twistors treated as elements of STA and classically ex- 
pressed in terms of one dotted and one undotted Pauli spinors. Authors recall that 
a physical point in the spacetime can be (classically) constructed in terms of an- 
tisymmetric bilinear forms in twistorial components 1 and show how this can be 
implemented in MSTA. 

Majid: The author explores the idea that spacetime itself may have noncommuta- 
tive coordinates and gives three such models based on Lie algebras viewed as 
“noncommutative spacetime”. Other examples of noncommutative spaces dis- 
cussed are the famous “noncommutative torus”, 0-spaces, and, unusually, Clif- 
ford algebras themselves with the familiar 7 -matrices viewed as noncommutative 
coordinates. In all these cases “quantization” may be achieved by means of a 
Moyal-product type cocycle deformation from a commutative space, in the Clif- 
ford algebra case as quantization of the finite lattice (Z^) 71 . As well as several 
mathematical results, of which the main one is a full account of how differential 
forms on such noncommutative spaces can be obtained by a cocycle twist of the 
classical differential forms, the article also contains a fully worked example in 
which noncommutative geometry is used to quantize 17(1) Yang-Mills theory on 
a finite lattice of 4 points in a path integral approach. 

Owczarek: This paper is a survey of difficulties that need to be overcome on the 
way to unify S. L. Woronowicz’s approach to the noncommutative geometry on 
quantum groups with A. Connes’ noncommutative geometry. One problem is to 
provide a proper description of spectral triples on quantum groups. Another is 
to find a way to define Dirac operators on quantum homogeneous spaces. The 
author outlines difficulties of defining the Dirac operator with Clifford algebras. 
He then moves on to compare several (similar) definitions of quantum homoge- 
neous spaces and related problems, e.g., problems with differential calculi on such 



1 Editor’s comment: There exist two additional projection formulas that give a spacetime point in 
terms of symmetric (bilinear) and hermitian forms, respectively, in two twistors. The original idea goes 
back to J. Rzewuski (Wroclaw, Poland) while the actual formulas, derived entirely using coordinate- 
free geometric algebra formalism, appeared in [10]. For another view on the Penrose’s projection, see 
J. Rzewuski [11]. 
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spaces. He reviews results by P. Podles on this issue. As an example, he discusses 
the quantum a 2-sphere of Podles as a homogeneous space of the quantum SU (2) 
of Woronowicz and shows that a natural Dirac operator on such 2-sphere does not 
satisfy one of Connes’ axioms. At the end he gives an outlook on possible future 
developments and concludes that axioms of Connes should be modified. 

Pozo and Parra: An “r-fold Clifford algebra” r Cl p , q of r-fold multivectors is 
defined by the authors by tensoring r Clifford algebras Cl PA . Grassmann and 
Clifford products, grade involution, reversion, conjugation, and the Hodge dual- 
ity can be extended to r Cl Ps q by applying them to each tensor slot. They contrast 

V 

C£ Pyq with a multiparticle geometric algebra (MPGA). Then, they discuss the V 
differential operator on r-fold multivectors and apply their formalism to superen- 
ergy tensors, which are important objects in General Relativity. 

Trayling and Baylis: The authors review a recent approach to the standard model 
in terms of the geometric algebra Ci 7 . The Dirac algebra is then embedded in Ci 7 . 
In particular, they demonstrate how the exact gauge symmetries of the minimal 
standard model are related to a Poincare group in a linear space with four extra 
spacelike dimensions. In this model, the gauge symmetries arise from the new 
algebra as unique symmetry groups of the current. The authors show how spinors 
used to represent all fermions of a single generation can be combined into one al- 
gebraic spinor of Cl 7 ; how the currents can be computed from such spinors, and 
how contributions from individual spinors can be identified. They further discuss 
the rotational symmetries of these spinors and derive from them the gauge sym- 
metries of the standard model. Finally, they show how a minimal Higgs field with 
the correct weak hypercharge assignment emerges in this approach. 



PART V. APPLICATIONS IN ENGINEERING 

Gebken, Perwass and Sommer: Due to expanding uses of (orthogonal) Clifford 
algebras Cl p , q in applied fields, i.e., robotics, computer vision, image process- 
ing, in which coordinate-free computations with multivectors (blades) replace 
traditional matrix computations, increasing computational efficiency with these 
algebras is of paramount importance. The authors present a first-ever design of 
a Clifford (geometric) algebra co-processor and discuss its implementation on 
a Field Programmable Gate Array (FPGA). The design allows for soft chang- 
ing the signature (p, q), the dimension n = p + q, and the numerical accuracy. 
High calculation speeds are claimed via a pipeline architecture. They also review 
briefly existing software packages, both numeric and symbolic, for computations 
with Clifford algebras. 2 Internally, the co-processor represents multivectors as 



2 Editor's comment: See also [4, 5] and, more recently [6, Appendix] (and references therein) for 
an extensive review of various software packages for these algebras. 
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lists of basis Grassmann monomials with associated scalar factors. The geometric 
product of two basis monomials is computed directly rather than retrieved from 
pre-computed multiplication tables, as this is computationally superior and more 
efficient. The authors then present their hardware design of the co-processor. The 
main algebraic operation implemented in the co-processor is, of course, the ge- 
ometric product in C£ Piq . Others are: the inner and the outer products, addition, 
subtraction, and operations between two multivectors. This presents demands on 
the arithmetic logic unit (ALU) that the authors discuss. Then, they present tech- 
nical aspects of FPGA, how the direct computations are performed on a chip, and 
then go on to the basis blade pipeline architecture . 3 

Koenderink: The author begins with a premise that image processing techniques, 
for example, in medical imaging alone, are “inconsistent” with basic principles of 
physics because images are surfaces in E 3 (called “image space”) where two di- 
mensions give a localization of a pixel on a surface, while the third dimension 
serves as the “intensity dimension”. As such, the latter is not interchangeable 
(commensurable) with the first two. The author proposes to look at the intensity 
dimension z as a (full) affine line (fiber) attached to the picture (Euclidean) plane 
(x, y). This construction yields a trivial fiber bundle in which sections are images, 
or, more precisely, Monge parameterizations of the “intensity surface”. Since the 
author considers images to be invariant under action of a group of movements and 
similarities, it is natural now that he postulates the image space to be a homoge- 
neous space. This rules out the Euclidean group as the group of proper motions. 
A pixel then is the base point of a fiber while the latter is just a line parallel to the 
intensity dimension. The result is that a bundle of parallel lines is preserved by a 
unique three-dimensional homogeneous space, one of the twenty seven Cayley- 
Klein geometries that has a single isotropic dimension. An eight-parameter group 
Gg of similarities acts on the image space preserving the fibers. Two of the eight 
parameters are moduli of similarities of the first and the second kind, with the 
remaining six parameters needed for a group of proper motions Gg. Similarities 
of the image space, upon projection onto the image plane, give the Euclidean 
similarities in that plane. Then, in order to understand the Gg similarities, the au- 
thor gives a detailed description of how one-parameter subgroups of Gg act in 
the image space: some appear like rotations, and some like translations in the pic- 
ture plane. Others act trivially in that plane while acting as, for example, “gamma 
transformations” in the fibers. Since the Cayley-Klein geometry is dependent on 
one isotropic vector and the image plane is the Euclidean plane, one obtains, in 



3 Editor’s comment: In Clifford algebras C£(B) of an arbitrary bilinear form B (not covered by 
this hardware design), computational complexity is of many orders higher than in C£ Piq . A substantial 
improvement in speed of symbolic computations in C£(B) has been achieved by using non-iterative 
Rota-Stein eliffordization based on Hopf algebra theory (see [8, 9] and the help pages of the “Bigebra” 
package [7] for further literature). It would be of great benefit to quantum physics to implement the 
Clifford product from C£(B) via a hardware design similar to the one presented in this chapter for 
Cl p , q . 
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fact, a Clifford algebra degenerate in one dimension C^2,o,i- The author describes 
consequences of this approach to the image space and pixels, explains why the 
similarities of the image space have two moduli, and how dual numbers (a suitable 
subalgebra of the Clifford algebra) are related to the so-called “normal planes”. 
He discusses the geometry of the latter and then moves on to discuss the geometry 
of the entire image space. Finally, he describes the differential geometry aspects 
of images and points out the relevance of his approach in the image presentation 
business. 

Rosenhahn and Sommer: Pose estimation is a task of determining a position 
and an orientation of a 3D object in a reference camera coordinate system from 
its observed 2D image. Since mathematical structures involved are the projective 
plane, projective space, kinematic space and Euclidean space, the authors use a 
conformal geometric algebra (CGA), as it contains all geometric algebras model- 
ing these spaces as subalgebras. In this chapter the modeled objects are cycloidal 
curves (circles rolling on circles or lines). As a (CGA) they use a Clifford algebra 
Ci 4,1 (denoted here as Ga,i) over the vector space R 4,1 = R 3 ® R 1,1 that con- 
tains an identified hyperbolic plane R 1,1 with a null basis. The latter is needed 
to identify a point at infinity and the origin. Products of Clifford exponentials of 
bivectors in Ct^ \ (rotors and translators) called “motors” give rigid transforma- 
tions (elements of SE{ 3)) in R 3 . The authors show how cycloidal curves can be 
generated in CGA by means of twists and discuss how pose estimation of these 
curves can be accomplished in that algebra. In particular, they investigate 2-twist 
and 3-twist curves. The chapter ends with a display of experimental results on 
pose estimations of convex and non-convex objects. 
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The Morera Problem 
in Clifford Algebras and 
the Heisenberg Group 

Carlos A. Berenstein, Der-Chen Chang, 
and Wayne M. Eby 



ABSTRACT In an open subset Cl of the complex plane C the Morera theorem gives 
a simple looking necessary and sufficient condition for a continuous function / to 
be holomorphic in Cl. Namely the vanishing of all the integrals f f(z) dz , where 7 
is an arbitrary Jordan curve in Cl whose interior also lies in Cl. The Morera problem 
consists in finding relatively small families T of Jordan curves such that the van- 
ishing of the corresponding integrals still ensures that the conclusion of Morera’ s 
theorem still holds. In this lecture we discuss a number of recent results obtained in 
the case one considers functions defined in two fairly different settings, the Clifford 
algebras and the Heisenberg group. 

Keywords: Morera theorem, Pompeiu problem, Clifford algebras, Heisenberg group 



1 . 1 Introduction 

Let us summarize briefly what is known about the Morera problem in the case of 
the complex plane C. We refer to the excellent survey article by L. Zalcman [26]. 
(See also [27].) The theorem of Morera asserts that if / is a continuous function 
in C, then / is an entire function, i.e., / E H( C), if and only if for every Jordan 
curve 7 , 



/ f(z)dz = 0 . ( 1 . 1 ) 

The Morera problem consists in finding “small” families T of Jordan curves such 
that if (1.1) holds only for 7 E I\ one can still conclude that / E H( C). 

Let T = { 7 } be a family of Jordan curves in the plane C and assume that T 
is invariant under the action of the group A of all affine, orientation, and area 
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preserving transformations of C, and let / be a C 1 function that we are testing for 
holomorphicity. Then, by Stokes’ formula 

f f(z) dz = f jjz(z)dz Ndz. (1.2) 

J'T J int7 

So that if / passes the test, i.e., the left-hand side of (1.2) vanishes for all 7 £ T, 
then the corresponding area integrals of the continuous function g — || vanishes 
also. 

Assume that T is generated over the group A by a single curve 70 , and let 
fio = int(7o) be the corresponding Jordan domain. Then it is elementary that 
what one needs to verify is the Pompeiu property. That is, if g G C(R 2 ) and 
satisfies 




for all 



c t E A, 



(1.3) 



then g = 0 . 

In fact the Pompeiu property for Qo and the Morera property for <9^o are really 
equivalent. What is really useful to decide whether dfio has the Morera property 
or not is the following equivalence with Schiffer’s overdetermined boundary value 
problem, as defined in the following proposition. 

Proposition 1. Let C — fi, where is a bounded open set , C c is connected and 
dC is (at least ) Lipschitz . Then £1 fails to have the Pompeiu property if and only if 
there is an eigenvalue a > 0 and a function u on £) satisfying the overdetermined 
Neumann problem 



( A u + au = 0 , in 

\« = i, & = 0 , on an. (1 ' 4) 

See for instance [4, 5]. The following corollary allows us to find a plethora of 
regions satisfying the Pompeiu property. 

Corollary 1. Let Q be as above. Then if Q, fails to have the Pompeiu property ; it 
must be that dQ, is real analytic. 

For instance, any triangle, rectangle, or polygonal region has the Pompeiu 
property. Similarly, a simply connected domain with an infinitely flat point in 
the boundary will also have the Pompeiu property. It is easy to see that in the 
case where fio is a disk, then there is an infinite number of eigenvalues for the 
overdetermined Neumann problem; simply look for those that correspond to ra- 
dial eigenfunctions. Therefore a disk does not have the Pompeiu property. Let us 
also remark that one can prove that if there are infinitely many overdetermined 
eigenvalues, then fio is necessarily a disk [4]. The question of whether disks are 
the only domains for which there is at least one such eigenvalue, sometimes called 
Schiffer’s conjecture, remains open. 
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Although a single disk, with its translations, cannot help detect whether a con- 
tinuous function in C is identically zero if all its averages over the disks of that 
radius vanish, one can decide this question using two different disks as long as 
the quotient of the radii is not in a certain countably infinite set. 

Theorem 1 (Two Radii Theorem [24]). Let f E C(M 2 ) and let r 1? r 2 > 0. 

Suppose these radii satisfy the condition 

— Q(Jn(x) ) — {- • Ji(rA) = Ji(sA) = 0 for some A E C}. 

r 2 5 

Then f — 0 if and only if for r — r i and r 2 , one has 

I /(x) dm(x) — 0 for all y E M 2 . 

JD r { y) 

The condition on the quotients of the two radii is equivalent to the condition 
that the Fourier transforms of the measures associated to the two disks have no 
common zeros. The non-existence of common zeros of the Fourier transforms 
of the associated measures is a necessary condition for the Pompeiu property, 
though not in general sufficient. Note that the disk is the only shape for which 
integrals over the rotations of that shape do not provide any extra information. 
When S belongs to a large class of domains, including those with a corner, it has 
been demonstrated that the collection {p a ^ S )}ae rotations does not possess any 
common zero [16]. However it is intriguing that disks are a special case, for the 
following reasons. This collection then degenerates to one element, ps, and it is 
known that for any domain S', ps will have zeros [16]. Thus we see why two 
disks are required for the Pompeiu property and have explained the source of the 
conditions for the radii in terms of zeros of the associated Fourier transforms. 
Note that circles, taken with Lebesgue measure, present a situation comparable 
to disks. Later a different method of avoiding the common zeros of the Fourier 
transforms of the associated measures will appear: the use of moments. 

It is natural to ask whether one can study the Pompeiu and the Morera prob- 
lems on a local region rather than requiring information on the whole plane. That 
is, suppose that / is a continuous function in the open unit disk D, is a Jordan 
domain in D. Then we would like to test whether / is zero (respectively holo- 
morphic) based on the vanishing of integrals of / over all copies of Q that remain 
within the disk D. Clearly, one cannot consider translations only because one will 
then really only test points near the center of D. The natural thing is to consider 
all Mobius transformations of D, that is the group M ~ SU{ 1,1) instead of the 
group A , as done above. For instance one obtains a result similar to the above for 
the Morera problem. 

Theorem 2. Let Q C D be a Jordan domain of class (7 2,e for some e > 0 and 
suppose that the Jordan curve V — dQ is not real analytic . Assume f E C(D) 
satisfies 

j f(a(z))dz = 0 for every o E M. 

Then f is holomorphic on D. 
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There is also a Pompeiu type theorem in this setting; one only needs to replace 
the Lebesgue measure by the .M-invariant measure in D, 



dp = 



dz A dz 



Investigating the Pompeiu problem on bounded open subsets of the ambient space, 
there is a related property, called the local Pompeiu problem, in which only a sub- 
set of the group (Euclidean or Mobius) is used. This question is discussed, for 
instance, in [11, 12]. In fact, similar ideas allow one to prove two-radii Pom- 
peiu theorems for all non-compact irreducible symmetric spaces of rank 1 [14]. 
In the same spirit, it is natural to consider these questions for variants of the 
Cauchy-Riemann equations. One such variant appears in Clifford analysis, re- 
lated to the study of the Dirac operator; another appears in the context of the 
boundary Cauchy-Riemann operator, that is, in the Heisenberg groups. It turns 
out that when one studies the “natural” Morera problem in Clifford analysis, it is 
easily reduced to the Euclidean case, as noted in [22]. 

On the other hand, for a holomorphic function / in the complex plane, we 
clearly have that for any r > 0 and all zq G C, 



1 

27 r 




+ re iB )e inB d0 = 0 



for any n > 0, not just n — 0. In other words, all the non-negative moments of / 
vanish on the circle. The converse is also correct. In [25], Zalcman investigated 
whether the converse of this result holds when only a few of the moments are 
considered and obtained the following theorem. 

Theorem 3 (Two Moments Theorem [25]). Let f G L 1 1 oc (R 2 ) and let r > 0 be 

fixed. Suppose there exist integers m , n such that for almost all w G C, 

/•27T r2lt 

/ f{w+re iB )e inB d0 = 0, / f{w + re lB )e imB d0 = 0 . 

Jo Jo 



Then 

(a) if 0 < n < m, f agrees almost everywhere with a solution of 

4 )”/ = °; 

(b) if 0 > n > m, f agrees almost everywhere with a solution of 

4 > W/ = 0; 

(c) ifn>0>m,m^ —n, f agrees almost everywhere with a solution of 
the pair of equations (^) n / = 0, (^) |m| / = 0. Thus in this case f is 
(essentially) a polynomial. 
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For instance, the case 1 — n < m characterizes when / is a holomorphic 
function. Similarly, 0 = n < m characterizes when / = 0. 

It is important to discuss the concept of homogeneity in relation to these mo- 
ments. This concept is directly related to the moments themselves, as well as the 
conclusions which may be drawn from them. For example, the two moments con- 
dition in Theorem 3 is equivalent to the condition that, for all w E C, 

f z rn f(w + z)da(z) = 0 and f z n f(w ^ z) da{z) — 0 
J\z\=r J\z\—r 

when n, m > 0, letting a be area measure on the circle in C. Or, if both n, m < 0, 
the conjugate moments turn up in the statement: 




z m f(w + z) da(z) = 0 



and 



\z\=r 



z n f(w + z) da(z) = 0. 



Note the distribution given by 



<f,Tm>= / z m f(z)dc{z) 
J\z\—r 

is of homogeneity m, while f m , given by 



< /, Tm >— / z m f(z)da(z ), 

J | z\ — r 

is of homogeneity -m. Note the conclusion of the theorem matches the homo- 
geneity of the smaller moment. When considering the moments on circles in C 
represented by the distributions T m and T n , the Fourier transforms are 

Ai(0 = r-/m(r|C|) and £ 2 (0 = f 1 J n (r\£\). 

When \n\ 7^ |m|, say \n\ > |m|, the Bessel functions an d ^n(r|£|) 

have no common zeros. Therefore the only common zero between fii and fa is 
£ = 0, with a multiplicity of m. These facts lead to the conclusion (J^) m / = 0. 
The connection between common zeros, moments, homogeneity, and the Pompeiu 
problem on spheres in C has now been described. 

Notice that when the lower moment is m = 0, then the two moments theorem 
gives the conclusion / = 0, similarly to the two radii theorem. However the 
two moments theorem has one advantage, namely this theorem is valid for any 
sphere of one fixed radius. In contrast, the two radius theorem does not work for 
many pairs of radii. In fact, the set of bad radii is dense. Further, by adjusting the 
lower moment, the two moments theorem can be utilized to detect more subtle 
properties of / : namely, 

It is of considerable interest that differential properties, such as these just men- 
tioned, may be tested by use of integral conditions. In closing let us remind the 
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reader that there are other examples of theorems using moments and yielding 
differential relations for a function. Zalcman has an earlier version of a theorem 
involving moments [24], and instead of getting the results by use of two moments, 
it uses one moment on two spheres of distinct radii satisfying certain conditions 
comparable to those of the two radii theorem. This theme of moments and differ- 
ential conditions will be explored further in the coming sections. 



1 .2 The Morera problem and moment problem in 
Clifford analysis 

In this section we present an overview of results on the moment problem in Clif- 
ford analysis first presented in [21]. Before discussing these problems in the Clif- 
ford analysis setting, it is necessary to introduce the Clifford algebra Cl n and 
present results relevant to the Morera problem, such as the version of the Morera 
theorem which appears in Clifford analysis. This material is adapted from [22], 
and further information may be found in the references [15, 20]. Let us first re- 
call that the Clifford algebra C£ n is generated from R n by the basis elements 
ei, . . . , e n of R n and the anti-commutation relations 

“l - Qj&i 25ij. 

The algebra Cl n will have basis elements 

1 , G i , . . . , G n , . . . , G j x • • • G j r , • • • , G i G n 

for all possible collections of indices 1 < ji <••<>< w and 1 < r < n. An 
element of Ci n may then be written in terms of the basis elements as 

a — 

a 



where o Q Gl and e a = e jl • • • e jr corresponding to a = (e jl , . . . , e jr ). 

There are two involutions which are defined by their actions on the basis ele- 
ments. These are given by 

rs -'. Cl n * C£ n . Qj l G j r > G j r • • * G j l , 

where a is written for ~ (a), and 

• C£ n > • Gji Gj r > ( 1) Gj r • • ’ Gj^ , 

where a will be written for —(a). 

The Clifford algebra Ci n becomes a Hilbert space and a Banach algebra with 
the inner product given by 



< a, b >= 






a 



a 
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for any a, b E C£ n . The space of functions to be considered will be the space 
of smooth C^-valued functions £(R n ,C£ n ). This space is also a C£ n module 
under pointwise multiplication. More generally there is the corresponding space 
of continuous functions C(R n , Ci n ). 

The notion of left and right regular functions replaces the concept of analytic 
functions as used for complex space. To define these we need to define the fol- 
lowing differential operator known as the Dirac operator. 



D =£ 



d 



e 



! dx. 



A function / € C l (R n , Cl n ) is called left regular if 



d f_ 

dxi 

i= 1 i= 1 a 






OXi 



and similarly g £ C l (R n , Ci n ) is called right regular if 



i= 1 



i—l a 



Note that the involution ~ yields Df — —fD so that / being left regular is 
equivalent to / being right regular. 

Note that D 2 = — A, so that regular functions are harmonic. Also consider the 



function 



G(x) 



1 —x 

LO n ||x|| n 



1 X 1 

U) n llxjl^ 2 ’ 



where LJ n is the surface area of the unit sphere in M n . This function is both left 
and right regular, and furthermore it takes the place of the Cauchy kernel. This 
observation leads to the following theorems [22]. 



Theorem 4. Let M be a bounded domain with Lipschitz boundary. Then for f E 
C l (U,Ce n )andxeM, 



/(x) = / G(y — x)n(y)f(y)dS(y) — / G(y - x)Df(y) dv(y). 

JdM J M 

The Cauchy integral formula for the Clifford algebras is a consequence of this 
identity, as the second term on the right-hand side vanishes when / is left regular. 
The following takes the place of the Morera theorem. 

Theorem 5. Iff is a Clifford algebra-valued continuous function on the domain 
U such that 

[ n(y)f(y)dS(y) = 0 
JdM 

for every bounded domain M in U with Lipschitz boundary, then f is left regular. 




10 Carlos A. Berenstein, Der-Chen Chang, and Wayne M. Eby 



With the Morera theorem in place, it is desirable to present what form the Mor- 
era problem takes in Clifford analysis. The following result [22] establishes that 
the relationship between the Pompeiu and Morera problems in Clifford analysis 
is the same as it is in Euclidean space. Let {Sj} be a collection of Jordan surfaces 
and let Mj = int Sj . 

Theorem 6. {Sj} has the Morera property in C£ n if and only if {Mj} has the 
Pompeiu property in R n . 

Therefore the surfaces which possess the Morera property are the same whether 
in the underlying Euclidean space or in the setting of Clifford analysis. 

It is possible also to consider the Morera problem with respect to Mobius trans- 
formations for Clifford analysis. They are the functions f acting on x E C£ n by 

<j)(x) — (ax + b)(cx + d)” 1 , 

where a, 6, c, d E Ci n and are assumed to satisfy the two conditions: 

(i) ac, cd, db , and da E R n , 

(ii) ad — be = ±1, 

and the transformation </> defines then a Mobius transformation formula over 
R n U {oo}. The corresponding 2x2 matrices, are called the Vahlen 

matrices. 

The following theorem [23] then describes the change of variable for these 
Mobius transformations in Clifford analysis. 

Theorem 7. Suppose that y = </>(x) = (ax + b)(cx + d)~ l is a Mobius transfor- 
mation and that f and g are Clifford algebra-valued functions. If S is a closed , 
bounded and oriented surface, then 



[ 9(y)n{y)f(y) dD( y) 

Js 



j 

J4, 



-HS) 



g(4>{x))J {(j>, x)n(x)J(4>, x)/(<£(x)) dS{x), 



where 



J(<A,x) 



cx + d 
cx + d\\ n * 



Finally it is possible to consider the Morera theorem under Mobius transforma- 
tions of the disk in the setting of Clifford analysis. 



Theorem 8 ([22]). Let f be a continuous Clifford algebra-valued function defined 
in the unit ball B in R n and let S be a Jordan surface inB. If 



L 



J {4>, x)n(x) x)f({<j>(x)) dS(x) = 0, 
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for every (f) £ M, where M is the group ofMobius transformations of the ball , 
then f is left regular if and only if M = int S has the Pompeiu property with 
respect to M . 

Finally the moment problem will be considered in the context of Clifford analy- 
sis. A set of moment conditions is presented which characterize those continuous 
functions which are left regular. 

Theorem 9 ([22]). Let f : R n — > C£ n be a continuous Clifford algebra-valued 
function. Let r > 0 be fixed. If for each x E R n , 

f n{y)f(y) d s{y) = o, 

JdB{*,r) 



and 



for 



/ dB(x,r) 



Pi(y - x)n(y)/(y) dS( y) = 0, 



Pfx) = eiXj + eiX X , i — 2, . . . ,n, 
then f is a left regular function. 



1.3 The moment problem in the Heisenberg group 

Let us look at some of the recently established results for the moment problem 
in the setting of the Heisenberg group. But first, as a reference point, it would be 
useful to describe the known results about the Pompeiu and Morera problems on 
the Heisenberg group which do not involve moment conditions. A good reference 
both for these results and for the methods used on the Heisenberg group is the 
book [10]. 

Formally the Heisenberg group H n may be defined by the collection of points 
{[z, t] G C n x R} under the group law 

[z, t] • [w, s] = [z + w, t 4- s + 2 Im(z • w)] . 

However, it is easier to visualize by identifying it with a manifold in C n+1 (of real 
codimension 1). This manifold {(z, z n +i) G C n+1 : Imz n+ i = |z| 2 }, which is 
the boundary of the Siegel generalized upper half space 

9, n - {(z ,z n+ i) e C n+1 : Imz n+ 1 > |z| 2 }. 

There is a natural group action of H n on dCt n giving an identification between 
these two objects, both as spaces of points and as groups. The map that identifies 
one to the other is given by 

$ : H n — > dCl n , [z, t\ i— ► (z, t + i\z\ 2 ). 
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Thus, there is an effective method to study H n as the boundary of the Siegel 
generalized upper half space dCl n . The set of left-invariant complex vector fields 
Zj , Zj , and T, defined by 






_d__. d_ „ _ T _^_ 

dzj lZj dt ’ 3 dzj lZ] dt ’ dt 



forms a complex basis for the Lie algebra f) n of left-invariant vector fields on H n . 
These vector fields have the following commutation relations: 



[Zj, Zj] — 2 i5jkT, 

and all the other commutators among Zj,Zj and T vanish. These vector fields 
endow the space H n with a CR structure. At some places in the paper we will use 
the notation Z to represent Z\ • • • Z n , so that Zf represents successive applica- 
tion to the function / of the vector fields Z\ out to Z n . Similarly, Z m represents 

A class of functions on H n known as CR functions play a role comparable 
to the holomorphic functions on the Euclidean space C n . Similar to the Cauchy- 
Riemann equations f = 0 for j = 1 , . . . , n satisfied by holomorphic functions 
on C n , a function / on H n is said to be a CR function when it satisfies the equa- 
tions Zjf = 0 for j = 1, . . . , n. CR functions are closely related to holomorphic 
functions as can be seen in the following result. 

Proposition 2. A function f on H n is a CR function if and only if there exists a 
function F holomorphic on the upper half space and continuous to the boundary, 
F £ H(Q n ) n C(Q n ) such that / o $ = F\dn n . 

This relation could be stated as: CR functions are the boundary values of func- 
tions holomorphic on the upper half space Cl n C C n+1 . The Morera and moment 
problems on H n will deal primarily with CR functions and means of characteriz- 
ing them by the vanishing of certain integrals. Similarly, the moment problem can 
be used to characterize functions satisfying certain differential relations, which 
are now expressed in terms of the complex vector fields Zj and Zj . 

On the Heisenberg group the function spaces being considered become more 
of a matter of importance. The first result proven in this direction was a Morera 
type result and was proved for the function space L 2 (H n ). 

Theorem 10 ([2]), Let f £ C l ( H n ) 0 L 2 (H n ). Then f is a CR function if 
and only if f satisfies the following integral conditions: For all g £ H n and 
k = 1 , . . . , n, 

[ L g /(z,0)w fc (z) = 0 

J\z\=p 

with tUfc(z) = dz\ A • • • A dz n A dz\ A • • • A dzk A • • • A dz n . 

The methods here are based on decomposition into a Laguerre series. The meth- 
ods are much more complicated than the corresponding results for L 2 functions 
in Euclidean space, which require only the Paley-Wiener theorem. This result 
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requires only one sphere of a single radius, but this is no different from what hap- 
pens in Euclidean space when the function space is also restricted to L 2 . Spheres 
of suitably chosen radii are again required when considering the function space 
L°°(H n ). Note that the conditions now require avoidance of zeros of Laguerre 
functions, as well as zeros of Bessel functions. 



Theorem 11 ([2]). Let f be a bounded continuous function on H n which satisfies 
for every g £ H n , 



/ (L g /)(z,0)doy(z) = 0 

JS(r) 

for two radii r \ , Here a r is the area measure ofS n (r). Assume that the above 
radii r\ and V 2 satisfy the following two conditions: 



1. 

2. % i Q(%#). 

Then f = 0. Conversely, if one of the conditions (1) or (2) fails to hold, then there 
is a bounded continuous function f ^ 0 on H n satisfying the integral conditions. 



Remark 1. A Morera version of this same theorem was proved in [1]. 

Let us now consider the more recent results about moment conditions. In mov- 
ing to a multi-dimensional setting, it is important to establish what kinds of mo- 
ments should be used to generalize those from C where the moments e irn6 d9 and 
e ~ ime dd correspond to terms z m ~ l dz and z rn ~ 1 dz. In C n the moments will be 
monomials z m = z™ 1 • • • zff ri of homogeneity m and z m = zff x • • • zff n of 
homogeneity — m, corresponding to the negative moments. These lead to conclu- 
sions comparable to the one-dimensional setting. 

Theorem 12 ([8]). Given a function f £ L 1 1 oc (C n ) satisfying 



/ /(w + z)z m do(z) — 0 

J |z| =r 



for all w £ C n , and for j — 1 , ,n and 1 < aj £ Z + , 

[ /(w + z)z m+ajej do(z) = 0. 

J\z\ = r 



We may then conclude 



gM 



In particular, for the integral conditions 



/ = o. 



' | z | —r 



/( w -f Z )zj do(z) =0, f /( w + z)z a j J d<j{ z) = 0 
J |z| —r 
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for j = 1, . . . , n and aj as before, we conclude that f is a holomorphic function, 



i.e., 



W> 



''dz. 



-f = o. 



It is also interesting to consider this kind of moment defined for a more general 
radial measure than the area measure on the sphere. The following result is in that 
direction. First, let us recall the Euclidean Laplacian A can be written 



A = 4^ 



3 = 1 



d 2 

dzj dzj 



Denote S(n, |m|, p) as the set of quotients of (nonzero) zeros of 



* i 

/i,n,|r 



i(*) = E 



(-i) A 



^ k\(k + £)\ 



M-\ z 



k 2 k+i 



k = 0 



where £ = n + Iml — 1 and 



A k = J t 2k+2e+1 



Now we may state the following result. 

Theorem 13 ([17]). Let f £ L 1 1 oc (C n ) and suppose there exist distinct positive 
real numbers r \ , V 2 such that 



rj 

J 0 J\z\=p 



z 



/( w + z) do( z) dp r (p) = 0 



(3.1) 



for almost every w 
on the sphere. If 
function g satisfying 



£ C n and r = ri,r 2 . Here do represents surface measure 
(f: S(n, |m|,/i), then f agrees almost everywhere with a 



0l m l 

— A s 

9(z,z) m 



9 = 0. 



Remark 2. In special cases, such as the area measure on the sphere or disk, 
the functions m | reduce to Bessel functions, whose zero sets are thoroughly 
understood. In these cases it is possible to study criteria which involve only one 
radius and multiple moments, such as the two preceeding theorems do. 

Returning to the Heisenberg group and the study of moment conditions, let us 
consider the setting of L 2 (H n ). There is a modification of the integral conditions 
considered above which allow us to conclude that a function is a CR function. In 
fact, moment conditions turn out to hold also. Further, with the introduction of 
conjugate moments (and corresponding differential forms) arises the possibility 
to conclude that a function is a conjugate CR function, that is, it satisfies the 
equations Zjf = 0 for j = 1, . . . , n. 
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Theorem 14 ([8]). Let f G C^H") n L 2 ( H"), and let m G (Z+) n . Then f is 
a CR function if and only if f satisfies the following integral conditions: For all 
g G H n and k = 1 , . . . , n 



/ z m L g /(z,0)w fc (z) = 0 
J\z\=p 

with cj/c(z) = dz\ A • ■ • A dz n A dz\ A • • • A dzk A • • • A 

In the case where — m G (Z+) n , / w a conjugate-CR function (it satisfies 

the equations Zjf = 0 for j = 1, . . . , n) if and only if f satisfies the following 
integral conditions: For all g G H n and k = 1 , . . . , n : 



z\=p 



z m i g /(z,0)a;/ c (z) = 0 



with u;/c(z) = cizi A • • • A dzk A • • • A A dzi A • • • A dz n . 

Beyond expanding the integral conditions to include moments, it is also pos- 
sible to consider more general radial measures, corresponding to the exploration 
of general radial measures in Theorem 13 for C n . In shifting from the differential 
forms of the above theorem to certain radial measures, we shift from the realm of 
Morera to the realm of Pompeiu. Other Heisenberg integral conditions involving 
radial measures include 

(1) Jj z|=r z m L g / (z, 0) do(z) — 0, 

(2) Jj z | <r z m L g /(z, 0) dfi( z) = 0, also written as 




|z|=s 



z m L g /(z, 0) do(z) ) s 2n 1 ds — 0, 



/n<|z|<r 2 zmL s /( z > °) d F( z ) = °> also written as 




z m L g /(z, 0) do(z) 



s 2n 1 ds — 0, 



(4) Jj zl=r z m L g / (z, 0) do(z) = |^"-i|/(g). 

The above integral conditions all consider area measures on a given radial surface, 
with the last one describing a mean value. The following two conditions introduce 
radial polynomial weights onto the measures. 

(5) J| z | <r z m L g /(z, 0)p(|z| 2 ) dfi( z) = 0, also written as 
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( 6 ) /r 1 <|z|<r 2 zmL s/( z > °)p(l z l 2 ) M z ) = 0, also written as 




|z|=s 



1 L g /(z, 0) da(z) p(s 



2\ g 2n-l 



ds = 0. 



Theorem 15 ([17]). Choose any one of the sets of integral conditions (l)-(6) 
and denote these integral conditions as (*). Let f G C l ( H n ) D L 2 (H n ), and let 
m G (Z + ) n . Then f = 0 if and only iff satisfies the integral conditions (*). 

Before considering the similar problem in the space L°°(H n ), we pause to 
consider the spaces L p ( H n ) for 2 < p < oo. We already know from Theorem 1 1 
that two radii are required for the space L°°. Recall that associated to the sphere 
{(z, 0) G H n : |z| = r} is an area measure <r r . Observe the following con- 
volution relationship in H n between the function and the measure ay. The 
T n -spherical Laguerre function ^ v is defined as 

n 

V£>,f) = e 2 - At e- 2w|A||z|2 n4?(4^|A||^| 2 ) (A,*/) G R* x (Z + ) n . 

j = 1 

The convolution relation is given by 

(V£ v (M) * <7r( z)) (w) = Cip^(r, 0) • ^„(w ,t). 

By choosing Ao so that 'ipfy(r, 0) = 0, the convolution becomes 

^(M) *CTr(z) = 0. 

In this situation we may say that is mean-periodic with respect to the mea- 
sure ay. As a consequence of this relation, and noting that ^^(z, t) G L°°( H n ), 
we may conclude a single sphere fails the Pompeiu property in L°°(H n ). Similar 
results are valid for a sphere with a single moment factor. The measures associated 
with these moments on the sphere are written as 

<(/>,a r m >= [ z m </>(z) do(z). 

J\z\=r 

Then the convolution becomes 

(i/^„(z, t) * a”(z)) (w) = /(A, r) • V£ + m,„(w, t). 

Thus, we can choose Ao so that /(Ao,r) = 0. In effect this yields a mean- 
periodicity relation, as before 

* cr™( z ) = 0. 

So when moving to consider moments in L°°(H n ), it will not be enough to con- 
sider only one sphere with a moment. 
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In both cases, whether with or without a moment factor, it is important to ask at 
what point it is necessary to move from one sphere to two. The following theorems 
answer that question, and interestingly, a main point in the proofs consists of the 
above convolution relations. 

Theorem 16 ([6]). Let r > 0 and let f E L p ( H n ) where 1 < p < oo. Then 
f = 0 on H n if and only if 



/ L g f(z,Q) da(z) = 0. 

J |z|=r 

The proof here is based on the Abel summability of the decomposition [18]: 

oo 

/ = lim ^2 r ' k| /*V’k,k 

r— > 1~ — 

| k |=0 

in L p - norm, for all / E L p ( H n ). This summation formula is used in conjunction 
with the relation 



(V£„(z, t) * a r { z)) (w) = V£„(r, 0) • t). 

Since the ^ reappears after convolution with a r , the same form of the summa- 
tion formula appears in the convolution / * a r , which is known to be zero. The 
summation formula is built from spectral projections, and it is then possible to 
show that each spectral projection is zero so that the function itself is zero. 

A newer result including moment factors is the following 

Theorem 17 ([19]). Let f E L p (H. n )for 1 < p < oo, and suppose r > 0. Then 
f = 0 on H n if and only if 



[ z m Lg/(z,0)dcr(z) = 0 for all g E H n . 

J I z I =r 

Similarly for the conjugate moment we may conclude f = 0 on H n if and only if 

[ z m L g /(z, 0) do(z) = 0 for all g E H n . 

J |z|=r 



The proof follows the same outline as that of the previous theorem, but there 
are some new difficulties in that the relation 

(v£,„(z, 0 * <C( Z )) ( W ) = /( A - r ) • V^ + m,„(w, t) 

introduces a shift in the indices. There is an additional problem with the Abel 
summation formula which is no longer valid for these functions i n place 

of the However this is remedied by considering convolution with the area 
measure on the sphere weighted by the conjugate moment. This additional con- 
volution shifts the index back to ground level and allows the proof to proceed. 
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Finally, when we consider bounded functions in H n , the need for conditions 
on two radii returns. Furthermore, at this level the moments introduce differential 
conditions into the conclusions, as they had in the Euclidean case. All of the 
remaining results are from [19] and are based on methods displayed in [1] and 
further described in [10]. First consider a moment problem with spheres of two 
distinct radii. 

Theorem 18 . Let f £ C(H n ) n L°°(H n ) and suppose that for all g £ H n , 

I z m L g /(z, 0)do(z) = 0 for i — 1,2. 

Js(n) 

Assume that r\ and r<i satisfy the following properties: 

(1) (£) 2 i e(£[f l+n_1) (*)). for alive (Z + )», 

(2) ^ ^ Q(Jn+\m\-l) = {f£ : ^n+|m|-l( 5 l) = ^n+|m|-l( 5 2) = 0}- 
Then we may conclude (Z 2m /)(z, t) = 0. 

Note that while the differential conditions are comparable to those for C n , there 
is a change from m to 2m. It is believed that this change is due to a limitation in 
the method, and further work is currently underway to establish the same results 
for m. The preceeding result can be specialized to a specific set of moments 
corresponding to CR functions, although a few additional steps are required in 
the proof. 

Theorem 19 . Let f £ C(H n ) D L°°(H n ) and suppose that for all g £ H n , 

[ ziL g /(z,0)da(z) = • • • = [ z n L g f (z, 0) do(z) = 0 fori = 1,2. 

Js(n) JS(n ) 

Assume that r\ and r 2 satisfy the following properties: 

(1) (£) 2 t Q(L$(x)). for alive (Z + ) n , A G R*. 

(2) % t Q(Jn) = {% : Jn(si) = Jn(s 2 ) = 0}. 

Then we may conclude f is a CR function. 

One of the advantages of the two moments theorem was the elimination of the 
requirement for two radii satisfying special conditions. The following results use 
n + 1 moments, as was the case for the moment theorem for C n . However the 
conclusions here are weakened to Z 2m /, as it happened in the previous theorems 
for H n . 

Theorem 20. Let f £ C( H n ) n L°°(H n ) and suppose that for all g £ H n , 
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and 



[ z m+ej L g / (z, 0 ) da(z) = 0 for j = 1 ,..., 

JS( r) 



n. 



Then we may conclude (Z 2 m /)(z, £) w a CR function. 

Once again there is a specific test for CR functions, requiring n + 1 specific 
moments. 

Corollary 2. Let f G C(H n ) D L°°(H n ) and suppose that for all g G H n , 

[ L g /(z, 0 ) do(z) = 0 
JS(r) 

and 

/ ZjL s f(z, 0) da(z) = 0 for j = l,...,n. 

Js(r) 

Then we may conclude f is a CR function. 



1.4 Additional questions 

Many possibilities arise when considering what may still be done in any of these 
directions. For the moment problem in Clifford analysis, note that so far only one 
set of moments has been developed. It is interesting to consider whether certain 
other sets of moments could be used to obtain conclusions regarding other differ- 
ential equations involving the Dirac operator. Likewise, it should be possible to 
consider the moments on more general radial surfaces. 

In the Clifford algebra setting, work has been done on the Morera problem 
when the surface of integration is a general Jordan surface 7 rather than only a 
disk or sphere [22]. It would be interesting to determine whether the same result 
holds for the Heisenberg group, though at present, the methods used are not well- 
adapted to non-radial surfaces. In the Clifford algebras it was also possible to 
consider Mobius transforms and the Morera theorem which appears in this setting. 
Of course, it would be interesting to investigate the same kind of issues in the 
Heisenberg setting. 

Finally, there are many other results known in the Euclidean case, and it is of 
considerable interest to find out whether an extension of them to the Heisenberg 
group or to Clifford algebras is possible. For instance, could the kind of local the- 
orems described in the paper [1 1] be valid if investigated in these settings? More 
particularly, for the Heisenberg group it would be nice to be able to investigate 
whether any result comparable to Proposition 1 could be established. This kind of 
result may help to determine whether the Pompeiu or Morera properties hold in 
the Heisenberg group for a larger variety of submanifolds. 
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Multidimensional Inverse Scattering 
Associated with the Schrodinger 
Equation 

Swanhild Bernstein 



ABSTRACT There is a deep connection between the theory of several complex 
variables and complex Clifford analysis. We will use a Borel-Pompeiu formula in 
C n and the representation of holomorphic functions obtained in the context of Clif- 
ford analysis to study the inverse scattering problem for an n-dimensional Schrodin- 
ger-type equation. Equations are found for reconstructing the potential from scatter- 
ing data purely by quadratures. The solution also helps elucidate the problem of 
characterizing admissible scattering data. Especially we do not need a “miraculous 
condition”. 

Keywords: Clifford analysis, multidimensional inverse scattering transform 



2. 1 Introduction 

2.1.1 Opening remarks 

The connection between nonlinear differential equations and certain linear spec- 
tral problems has motivated the study of various inverse problems. As it is shown 
in [2] and [3] by recasting the notion of scattering data as 9-data, it is possible 
to generalize the theory to higher dimensions. The use of the one-dimensional 
generalized Cauchy-formula in [1] which leads to n different reconstruction for- 
mulas causes the so-called “miraculous condition”, i.e., the inverse data must be 
of that kind that all reconstruction formulas give one and the same potential. The 
9 approach, which goes back to Beals and Coifman [2, 3] as well as Nachman 
and Ablowitz [8], does not need such a condition; neither does our Clifford anal- 
ysis based approach. A first attempt to use Clifford analysis was made by Dix [6]. 
His approach replaces “holomorphic function” by “monogenic function” which 
is equivalent to the replacement of the 9-operator by the Dirac-operator D. But 
the inverse scattering problem is deeply connected with complex analysis which 
cannot be replaced. Therefore complex Clifford analysis, which is connected with 
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(real) Clifford analysis as well as complex analysis in C and C n , is a good basis 
for studying inverse scattering problems. The product in Clifford analysis com- 
bines the scalar (or inner) product and the exterior (or wedge) product for vectors. 
For vectors a, b the Clifford product ab = -a • b + a A b. This combination 
explains in an easy way the crucial role of compatibility conditions for the in- 
verse scattering problem. The heart of the Clifford analysis approach is a general- 
ized Cauchy-Green formula also called the Borel-Pompeiu formula which in fact 
contains the Martinelli-Bochner formula as well as the compatibility conditions. 
Essentially, this formula will be used for reconstruction purposes and the condi- 
tions under which the reconstruction is correct. We give a rigorous treatment of a 
regularized Schrodinger equation using a Borel-Pompeiu formula in C n as well 
as the time-dependent and the time-independent Schrodinger equation where we 
use the Borel-Pompeiu formula in R n . 

It should be emphasized that we do purely formal manipulations to demonstrate 
the way our method can be used. What we must have in mind is that the operator 
identities are valid only on certain domains of functions. This problem should be 
discussed elsewhere. 

2 . 7.2 Real and complex Clifford analysis 

We want to start with the real Clifford algebra C£o,n which is generated by the 
elements ei, e 2 , . . . , e n . These generating elements fulfill the n< 9 n-commutative 
multiplication rule 

GjGj 4“ — ‘IdijCOi i,j = 1,2,..., 71, 

the unit element of this algebra is denoted by eo and commutes with all gen- 
erating elements, i.e., eoG* = e^eo, obviously e 1 = 1. Thus we can identify 
eo = 1. An arbitrary element of the algebra is now a linear combination of all 
possible products of generating elements, more exactly, a = a a e a, where 
Gq- — Gaia 2 ...o:fc = . . . e Q , k and 1 ^ ol\ ol 2 ol ^ n, olj G 

{1,2,..., n}. We identify G 0 = eo. The part Sc(a) = aoeo = ao is called the 
scalar-part and Vec(a) = a — Sc(a ) is called the vector-part of a. 

We will identify elements of the Euclidean space with elements of the Clifford 
algebra in the following way: 

R n 9 (xi,X 2 , . . . ,X n ) = X = ^2 X 3 e 3 ^ 

3 = 1 

Especially for a, b G W 1 we have 



ab = —a b + aAb, 

n 

where a • b = b • a = ^ ajbj denotes the usual scalar product of two vectors 
3 = 1 

and aAb = -bAa = ^(ab - ba) denotes the wedge or exterior product. We 
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introduce a first-order differential operator, the Dirac operator, which plays the 
role of the 8 operator in complex analysis: 



3 = 1 



d 

docj 



Functions / with Df = 0 will be called (left)-monogenic. Due to 



DD 



£' 



3 dxj 
3 ■ 



0 = 1 



dxj 



71 fi 2 



. , dxj 

3 = 1 J 



monogenic functions are harmonic. A fundamental solution of D is 

i n i 

— 1 XjQj — 1 X 



^n-l . = l | x r 



^n — 1 x f 



where cj n -i denotes the surface area of the unit sphere in M n . Using the funda- 
mental solution we are able to state a (general) Cauchy formula. 

Theorem 1. Let G be a bounded domain in R n with boundary T. Then 



UJ, 



,-i Jr l x 



(x-y) 



+ 



n{y)u(y)dT y 
y I " 

i f (x-y) 
Wn-l Jo l X - y|” 



(Du)(y)dy 



u(x), xgG 
0, x G M n \{0} * 



In case ofG = M n we have 

— — ^-{Du)(y)dy = tt(x) - tt(oo), (1.1) 

x - y| n 

where u( oo) = lim u(x). 

|x|— >oo 

Suppose u(x) is a scalar-valued function and equation (1.1) is equivalent to 
the system: 






- 

i-l JlT 



u(oo) 

and 



^n -1 JlRrx 



1 



<^n- 1 




(x 



(X 



- y) • W(y) 

|x - y|" 

— y) A (Du)(y) 
|x - y|" 



dy = u(x), 
■dy = 0, 



( 1 . 2 ) 

(1.3) 



where (1.2) gives the reconstruction formula and (1.3) describes conditions which 
have to be fulfilled for the reconstruction. Especially, (1.3) is equivalent to the 
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compatibility conditions. We have 



u n- 1 J 



(x — y) A (Du)(y) 
|x — y| n 



dy 



(n-2) 



/ 

) w n- 1 7]Rr 



D 



|x — y|” -2 



A(Du)(y)dy 



(n-2) 



1 f 1 A / d 2 u d 2 u \ 

)w„-i i R n |x — y |"~ 2 dxkdxj) ejGfc y 

3 — ±,3<.k 



and thus 



d 2 u 

dxjdxk 



d 2 u 

dxkdxj ’ 



Vj, A:, j ± k. 



For our considerations we will also need complex Clifford analysis. We introduce 
an analogue to the 3 operator 



T) — Z) x -j - iDy, 

n n 

where Z) x = ej and D y = e j . Then a holomorphic function 

3 = 1 J _ 3 = 1 J 

/(x, y) in C n fulfills the relation Vf = 0 because 

(l“ + i |-) /<x ' y) = 0 ’ i = 1 ' 2 -- n - 

We will identify 



(x,y) G M 2n ~ x + iy G C n 

~ (X,y) G ff 0 ,n(R) X C^ 0| n(R) 

~ x + iy G C£o,n(C). 

This allows us to combine methods from real, complex and Clifford analysis. 
More information about Clifford analysis can be found in [5] and [7]. The con- 
nection between (complex) Clifford analysis and the theory of several complex 
variables is investigated in [9]. 



2.2 The ^-operator 

Even though the operator V seems only to be a new “structure”, we will see that 
the view of holomorphic functions in several complex variables analysis will be 
considerably improved. The operator V is not an elliptic operator but if we restrict 
our consideration to scalar, i.e., complex- valued, functions, we are able to obtain 
a generalized Cauchy formula in C n (also called Borel-Pompeiu formula). The 
proof is given in [4]. We have 
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Theorem 2 . Let G E R 2n be a bounded strongly Lipschitz domain with boundary 
T and exterior unit normal n = n x 4- in y (we identify R 2n with C n , thus our 
normal has a real part in the x- direction and an imaginary part in the y -direction) 
defined almost everywhere on T and m E C 1 (G). Then we have for (x, y) E G, 



m(x,y) 



-!-[ 

<*>2n-l JC 



{Vm){^ r))(Z - x) - i(j] - y)(Dm)(Z, rj) 
(1C — x | 2 + \r) — y| 2 ) n 



dr] 






W2n-i Jr ' (l^-x| 2 + |r?-y| 2 ) n 

where (n^ + in^) dT = n dT . A/id in case ofG = R n we have with 

= , a Jim m(x, y), 

|x| 2 +|y| 2 -^oo 



m(oo) H- 



1 f (Pm)(£, 77) (£ - x) - z(t? - y)(X>m)(C, ??) 
W 2 n-1 i R 2n (|C — x| 2 + |r? — y| 2 )” 



dr] 



= m(x,y). (2.1) 



In fact (2.1) is a system of equations where the scalar part (2.2) gives the recon- 
struction formula itself and the vector part (2.3) gives the conditions under which 
the reconstruction can be done, i.e., (2.1) is equivalent to 



m(oo)+ 

— / 

^2n-l JR 



]) ■ [(£ - x) - i(ri - y)] 



(1C 



+ \v-y\ 2 Y‘ 



d£, dr) = m(x, y) (2.2) 



and 



-L- [ 

^2n-l JR 2" 



(Vm)(£, r?) A [(£ - x) + i(r] - y)] 
(1C — x | 2 + \r] — y| 2 ) n 



d£ dr) = 0. (2.3) 



Moreover the conditions (2.3) can be rewritten as the compatibility conditions 

— —TO - ——TO l ± ' 



2.3 History and principle of the Inverse Scattering 
Transform 

The historical background of the problem is deeply connected with the connection 
of overdetermined linear systems and non-linear equations as it will be explained 
in the following example. 

U xx (x, t, z ) + [q( x, t) + z 2 ]u(x , t, z) - 0, 

Ut(x, t, z) — [ z 2 -f s(x, t)]u x (x, t, z) — [ iz 3 -h r(x, t)]u(x, t , z). 
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Here z is a real or complex parameter; the coefficients q, r, s vanish rapidly as 
\x\ — > oo. Cross-differentiating, equating coefficients of powers of z, leads to 
the compatibility or integrability conditions linking the coefficients, which can be 
resolved for r and s in terms of q leaving the Korteweg-de Vries equation 

qt + \q x xx + \qq x = 0. 

With the asymptotic condition u(x, t, z) ~ e lxz as x — > -oo, u has the form 
u(x, t, z) ~ a(z, t)e lxz + b(z, t)e~ lxz as x — > -f oo 
and the time-evolution of the coefficients a and b is given by 

a t (z, t) = 0, b t (z , t) = —2 iz 3 b(z, t). 

Thus: The scattering transform q ^4 (a, b) from the Schrodinger potential q to the 
asymptotic data (a, b) linearizes the KdV evolution. 

This behavior is not exclusive for the KdV equation. Similar results are known 
for the following equations: 

(i) Modified Korteweg-de Vries equation 

Ut — 6u 2 u x + u xxx = 0. 

(ii) Cylindrical Korteweg-de Vries equation 

Ut + 6uu x + u xxx + ^ u = 0. 

(iii) Klein-Gordon equation 

u xt = e 2u — e~ u . 

(iv) Sine-Gordon equation 

u xt = sinu. 

(v) Boussinesq equation 



^tt ^xx H" (u ) xx zb u xxxx — 0. 

(vi) Non-linear Schrodinger equation 

iu t + u xx ± 2|u| 2 u = 0. 

All these equations have the common properties that the non-linear equation is an 
integrability condition for the coefficients of an overdetermined linear system. If 
solutions of that linear system are normalized by their asymptotic behavior, then 
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one can expect the scattering map — the relations satisfied by the asymptotics — 
to linearize the associated evolution. The analytic problem is to find a suitable 
normalization for solutions and to determine the scattering map and its inverse. 
To treat the analytic problem we consider the direct problem for a linear equation 
depending holomorphically on a complex parameter z : 

P z u = q(x)u , 

where P z is a differential operator in x. We look for Jost-type solutions with some 
kind of boundary condition: 

u(x, z) = m(x, z)e z (x), P z e x — 0, lim ra(x, z) = 1. 

X 

Then the corresponding equation for m is: 

P z m — q(x)m(x , z). 

Let G z be a fundamental solution to P z and rewrite our problem as an integral 
equation of Lippmann-Schwinger type: 

m(-,z) = 1 + G z [qm( - ,z)] 



which has the formal solution 

m(-,z) = {I -G z q)~ l \. 

Now the inverse problem consists in determining the 9-data 

dm = [5, (I - G z qy l \l = (/ - G z q)~ 1 [d, G z }q(I - G z q)~ 1 l 

= (I -G z q)~ l \8,G z )qm. (3.1) 

Because P z depends holomorphically on z, dm satisfies the same differential 
equation. If the ra( ■ , z) are a complete set of solutions, then dm should be ex- 
pressible as a combination of the m, i.e., 

dm — Tm , (3.2) 

where T is a multiplication or integral operator which depends (non-linearly) 
on q. In this case T should be computable from (3.1). With the asymptotic condi- 
tion m(x,z) — > lasz—> oo and the Cauchy transform C we get 

m = 1 + CTm or, respectively, m = (/ — CT ) _1 1. 

T commutes with P z and with left multiplication by functions of x. This ensures 
in practice, that T is determined by a functions of z alone. 
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Now, define functions m by m = (/ - CT ) 1 1 ; we get 

P z m = [P z , ( I - CT)~ l }\ = (I - CT)~ l [P z , CT}(I - CT)~ l l 
= {I-CT)~ 1 {P z ,C}Tm. 

Because P z depends in an affine way on z, from the form of C we see that the 
commutator [P Z ,C] maps to functions which are independent of z. Thus if we set 

q(x) = [P z ,C]Tm 



we get 

P z m = (I- CT)~\ = (/ - CTy'lql] = q(I - CT)~ l l = qm. 



Thus on a formal level the 9-approach provides an invertible correspondence be- 
tween the operators P z — q and T. 

How can this procedure be generalized to higher dimensions? Again we start 
with P £ u = q(x)u where P? is a partial differential operator in x and z G C n is 
again equivalent to P^m = q(x)m for Jost-type solutions 

u(x, z) — m(x, z)e?(x) 

for x = (xi,X 2 , . . . ,x n ). Then with m, also dm is a solution. Again we can 
conclude if the solutions m form a complete set, then the solution dm should be 
expressible as 



dm = Tm <=> m — 1 = dm = CTm m = (/ — CT) 1 1. 

If T commutes with P? and with left multiplication by functions of x, we are able 
to define m = (/ - CT ) -1 l and q(x) = [Pg , C]Tm in the same way as before 
and we obtain P?m = qm. 

What we need is an appropriate analogue to the Cauchy transform in complex 
analysis which is given by the Cauchy transform in complex Clifford analysis. In 
the next section we will demonstrate how the procedure works in higher dimen- 
sions. 



2.4 Inverse scattering associated with a 

multidimensional generalized Schrodinger 
equation 

In our considerations we follow [1]. We look for Jost-type eigenfunctions m of 
the form ve lk ^ x +l k l 2t / <r where k € C n , k • x = \ then m 

fulfills 



amt + Am + 2ik • Vm — um — 0. 



(4.1) 
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According to [8] we use the Green’s function 



G(t,x,k) 



i(tr+x-£) 



1+1 ioT - |£| 2 - 2k • £ 



d £ dr. (4.2) 



A solution m of (4.1) can be obtained by solving the integral equation 

m — 1 + G[um ), (4.3) 

where G denotes the convolution operator with kernel G(t — t' — x' ,k). We as- 

sume that (4.3) has no homogeneous solutions. Differentiating the integral equa- 
tion for m with respect to kj gives another solution of (4.1), which can be ex- 
pressed nonlocally in terms of m using the symmetry property: 

e -i/3(t,x, k«,0 G ( t) X) kj?) k/ ) = G(t, x, Z, k/) (4.4) 

2 2 

on s(Z) = Z H — -k i - kfi H — -k/ = 0 
or or 



with (3(t , x, kij, k/,£) = ( x + 2^4-k / j • (£ - k#) . We obtain the dj derivative 



dm 



f e i/3(t,x,k ^^_ kRi) 5 {s ^))T(k R ,k,,Z)m(t,x,Z,ki)dZ (4.5) 

J 



where the scattering data are found to be 



T(k/*,k/,£)= f e x )ra(£, x, k) dt dx (4.6) 

which defines the scattering data T. Here, following [1] we have used the d 
method introduced in [2]. Now, our improvement comes by using Clifford analy- 
sis. Applying the general Cauchy formula of Theorem 2 together with m — > 1 as 
\k\ -4oowe get 



m(*,x,k fl ,k/) = 1 + 



^2n-l JlR2r 



(Dm)(k' R ,k))(k' fl - k fl ) - i(k' f - kj)(X>m)(k' R , k' 7 ) , , 

(Ik' -k^ + Itt-k/R" * 



Using Theorem 2 we get a new motivation to follow the idea of M.J. Ablowitz and 
A.I. Nachman to derive the “characterization equations” from the compatibility 
conditions 

d d _ d d 
WiWj m ~ WjW™' * 3 ' 




32 Swanhild Bernstein 



Following [1, 8] we get a non-trivial non-linear equation for T : 

r m-(c l x ( dT , 1 dT \ i ^f dT 1 dT 

C..,\T\ - ft k Rl ) j a( J ~ ft - <=«,) 

= f [( r)j - k Rj )(ti - Tji) - ( r)i - k Ri )( - r]j)}5(s(r])) 

J R n 

x T( kfi, kj, 7?)T(k K , k/, £) [T] . 



To avoid redundancies we consider only i = 1 and introduce new variables 
(z, w, wo) G C n_1 x R n x R, where z = (z 2 , . . . , z n ) : 



=J2 w i ZR j - 

3 = 2 



1 j (JiWqWi , i (JlWQWj 



2 w 2 



k h = Y j w j z Ij +\ 
j=2 



<J R WqWi 
2 W 2 



1 a R w 0 Wj 

kh = -wi zj + x s-A 

J J z w z 



c I 1 jCT/Wo^x J , (JjWQWj 

^ = Zw + 2^1-2 — w 2 > 6 = + 2™J - 2 ,„ 2 • 



j=2 



Ur 



If wi ^ 0, there exists a one-to-one mapping (k#, k /, £) — > (z, w, u>o) such that 

2 9 

w = £ - k/e, wo = — k/-(£ — kjj), = Aj, 

a R dzj 

and, for i — 1, j = 2, . . . , n, we obtain 

^7 =M J 'P1(z,w,w 0 ). 



Applying the 22 operator of Clifford analysis we are able to write 

OT = V 1 [T](z,w,w 0 ), 

where 



(4.7) 



A/"i[T](z,w,wo) = y^A/ij [T](z, w,w 0 )ej 

3 = 2 

and Kr = ^e,-(A + i A lr . 



i=2 



9xj dyj 



Using Theorem 2 we have that if T(z, w, u; 0 ) fulfill the compatibility conditions, 
which is equivalent to a vanishing vector-part in the Borel-Pompeiu formula, we 
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can reconstruct the potential by using the scalar-part of the Borel-Pompeiu for- 
mula: 



I[T]{ z,w,w 0 ) :=T(z,w,w 0 ) 



1 



^2n-3 




2 



MP1(C,w, w 0 ) • [(Cfl - Z R ) - i( C/ - Z/)] 

(ICh-z^P + ICz-z/I 2 )"- 1 



d(R d(i 



— T( oo, w, w 0 ) = u(w 0 , w). (4.8) 



Due to the fact that |z| — > oo, also |k| — > oo and 

m(Z,x, k fl ,k/) = 1 + 0(|k| _1 ). 

Thus we have m — > 1 as |k| — > oo, 
lim T( k/*,k/,0 

|k|-+oo 

= lim f e -^(^ x , k R, k /,0 li (^, x )m(t,x,k^,k/) dxdt 

l k H°° JlRn+l 



= lim / e ^ x>w+tli;o ^ii(t,x)m(i,x, k^,k/)dx(ii 
l k H°° 7l^+i 

= T(oo, w,w 0 ) = u(w 0 ,w ) . 

Here u denotes the Fourier transform of u. Thus we have 

Theorem 3. //7/ze homogenous equation of (4.3) has no solutions and the data T 
fulfill (4.7) as well as the compatibility conditions and decay properly if\k\ — > oo, 
then the potential u can be computed by taking the inverse Fourier transform of 
I[T } (see (4.8)). 
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On Discrete Stokes and 
Navier-Stokes Equations 
in the Plane 

Klaus Giirlebeck and Angela Hommel 

ABSTRACT The main goal of the paper is to apply the theory of discrete analytic 
functions to the solution of Dirichlet problems for the Stokes and Navier-Stokes 
equations, respectively. The Cauchy-Riemann operator will be approximated by 
certain finite difference operators. Approximations of the classical T-operator as 
well as for the Bergman projections are constructed in such a way that the algebraic 
properties of the operators from complex function theory remain valid. This is used 
to approximate the solutions to the boundary value problems by adapted finite dif- 
ference schemes. 

Keywords: finite difference Cauchy-Riemann operators, discrete function theory, 
boundary value problems 



3.1 Introduction 

The idea to calculate solutions to boundary value problems by finite difference ap- 
proximations is one of the major approaches to numerical methods. There are dif- 
ferent lines to construct adapted methods. Conservation laws have to be preserved 
in the discrete model, variational principles are reflected by some finite difference 
schemes, and last but not least approximation (truncation error) and numerical 
stability play important roles. Beginning in the 1940s first attempts were made 
to consider solutions of finite difference approximations to the Cauchy-Riemann 
equations as a class of discrete holomorphic (discrete analytic) functions and to 
establish a discrete function theory ([2, 4, 5, 11]). The main problems were the 
definition of a discrete analogue to the Cauchy integral (and the equivalent defini- 
tion of the discrete analytic functions by this discrete Cauchy integral) as well as 
the property that discrete analytic functions do not form an algebra with the usual 
complex multiplication. Another problem was to find a factorization of the 2- 
dimensional real Laplacian by two (adjoint) discrete Cauchy-Riemann operators. 
So, for a long time there was no essential progress in discrete function theories in 



AMS Subject Classification: 30G25, 31C20. 




36 Klaus Giirlebeck and Angela Hommel 



the described sense. 

A series of works in Clifford analysis has shown that a commutative algebra 
is not necessary to adapt a lot of function theoretic methods to the solution of 
boundary value problems and to find representation formulas of the solution. A 
survey and a collection of examples is for instance contained in [10] and [14]. 
Also in [10], and later in [9] first ideas for the development and the application 
of discrete function theories are presented. This work was inspired by analogous 
ideas to develop a discrete potential theory ([3, 15, 17]). A first serious problem in 
the three dimensional case is to find discrete Cauchy-Riemann operators that fac- 
torize exactly the standard seven-point Laplacian. In [9] this problem was solved 
and the approximate factorization defined in [10] was improved. For the definition 
of the discrete T-operator and a discrete Cauchy integral we need a fundamental 
solution of the discrete Cauchy-Riemann operator. Using the factorization of the 
Laplacian one can work based on the fundamental solution of the Laplacian. This 
fundamental solution is well studied (see e.g., [1, 6 - 8 , 13, 16, 18-20]). An addi- 
tional problem in the two-dimensional case is related to the logarithmic behaviour 
of the fundamental solution at infinity. This asymptotic behaviour complicates the 
error and norm estimates in exterior domains and requires a careful definition of 
Cauchy integral and T-operator. We have to prove analogous algebraic properties 
of the operators. In detail we discuss in Section 2 that the discrete T-operator is 
a right-inverse of the introduced finite difference Cauchy-Riemann operator. A 
discrete Borel-Pompeiu formula is discussed and a suitable orthogonal decompo- 
sition of the / 2 -space is the basis to define discrete Bergman projections. Applying 
these tools with analogous properties as in the continuous case, we discuss finite 
difference approximations of the steady state Stokes and Navier-Stokes equations 
in Sections 3.3 and 3.4. 



3.2 Elements of discrete function theory in the plane 

Let R 2 be the Euclidean space with the unit vectors e^, i = 1,2 and G be a 
bounded domain. We consider functions f(rrth) = (fo(mh), fi(mh)) for fixed 
mesh width h > 0 on the lattice 

R 2 = {mh = (m\h,nri 2 h) : ra* G Z} 

and denote the discrete domain by Gh — (Gr 1 R 2 ). Using the notation 
K = {fci(l,0), k 2 = ( 0 , 1 ), h = (-1,0), h = ( 0 ,- 1 )} 
we define the discrete boundary 

= {rh G R£ \ G h : 3 fc* with (r + ki)h G G h , i = 1, . . . , 4}. 

Often the boundary 7 ^ is split into the parts 

= i rh e Th : ( r + ki)h G G h }, i = 1, . . . , 4. 
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At first we present an orthogonal decomposition of the space hiGk) of functions 
w(mh) = (wo(mh),wi(rnh)) and v(mh) = (vo(mh),vi(mh)) with the scalar 
product 

<wv>= T h 2 (Mmh)\ T (vo(mh)\ 

Outside of Gh we extend our functions by zero. Then the finite difference opera- 
tors are defined as 

D 3 h = h~ l {Uj — I ) , UjW l (mh) = w^nih -f hej ) , 

Df J = h~ 1 (I — U-j ) , U-jW x (mh) = w l (mh — he 3 ) , 

for j G {1,2} and z G {0, 1}. If the boundary conditions 

w(rh) — v(rh) = (0,0) , Vr/i G 7^, 
are satisfied, then we introduce the space w\(Gh) with the norm 

3 = 1 

Theorem 1. The space permits the orthogonal decomposition 

h(Gh) = ker D^ M (Gh) © D 2 M (w 2 {Gh)), 

where 



n 1 

u h,M 




and 



n 2 _ 

^h.M - 



^7 



^7 



Proof. Let w(mh) G and v{mh) G D\ ^(w^Gh))- At first we prove 

c ker D\ M {G h ). 



Suppose that 

< w, v >— 0 for all v(mh ) G D h,M(™2i G h))- 

Then, for 

s(mh) = (so(mh), si(mh)) G 

we get that < w, In D\ m s >= 0, where N denotes the number of mesh points 
mh G Gh and In denotes the A x A identity matrix. If we write the scalar 
product in the form 

0 = ^ h 2 [wo(mh) D\ So(mh) + w\(mh) Df l so(mh) 

mh£Gh 



+ (—wo(mh) D\ si(mh) + w\(mh) D h 2 si(mh))\ 
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and change the index of summation, then we obtain for functions w,s e ( Gh ) : 

53 h 2 Wi(mh)[D^ J s k (mh)\ = - ^ h 2 [D^ J Wi(mh)]s k {mh) (2.1) 

mh£Gh mh£Gh 

with j £ {1, 2} and i, k £ {0, 1}. We conclude 
0=<w,I N D 2 h M s > 

= — 53 ^ 2 [-D^ 2,u, o(m/i) + D^u;i(m/i))so(m/i) 

mh£Gh 

+ {-D^ l wo(mh) + Dlwi(mh))si(mh)] = - < I N D^ ^w, s > . 

Because the equality should be fulfilled for all s(mh) £ w\ {Gh), it follows that 
w(mh) £ ker D^ m (Gh). We prove now that also 

kel D^ M (G h ) C D h < M^2( G h)) 1 - 
For w(rnh) £ ker D^ M (Gh) we obtain analogously 

<w,v >=< w, I N D 2 hM s>=- < I N D l h M w, s >= 0. 

From both parts of the proof we get that 

terDl M (G h ) = D 2 htM (wl{G h )) x . 



We remark that the operators D\ M and M approximate the Cauchy-Riemann 
operators 



D l = H) 



d , d 
dxi dx2 



and D 2 = i 



d _ d \ 
dx\ dx 2 ) 



with x = (x\,X 2 ) £ M 2 and i 2 = —1. Consequently, functions belonging to 
the kernel of D\ M will be called discrete analytic functions. In analogy to the 
continuous case these operators have the property 

Dh,M DIm — Dh,M — ^2 A h, 

where I 2 denotes the 2 x 2 identity matrix and A h = D\ Df 1 + D\ Df 2 is a 
discrete Laplacian. 

Lemma 1. Let w,s £ w^iGh)- Then we have 

< w, I N D\ m s> = - <I N D l h M w,s> 

< w, I N Dl M s > = - < I N D 2 m w,s> . 



and 
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The first property is already shown in the proof of Theorem 1. Using formula 
( 2 . 1 ) the second property can be proved analogously. 

Let 



S h (mh) 



h 2 , for mh = ( 0 , 0 ); 
0 , for mh / ( 0 , 0 ); 



and E\(mh) with 



Df Dl\/El n (mh) El 12 (mh)\ = (6 h (mh) 0 \ 

- D f 1 D 2 h ) \El 21 (mh) E\ 22 (mh)) \ 0 &h{mh)J 



be the discrete fundamental solution of the operator D\ M . We define a discrete 
analogue to the complex T -operator 

(Tk[fo,mmh) = ((TUfoJi})(mh),(TUfoJi})(rnh)). 

The components of T £ are represented by 

(TUfoJi))(mh) = (T l h f[f 0 ,h})(mh) + {T^ [/ 0 , h])(mh) 



with 






E 

Ih^Gh 



E\ kl {mh-lh)\ T 

E hk2( mh - lh )J 



im) 



and 






E * 

lh€ 7 M U7 h4 U r u 



E hki( mh - lh)\ 
El hk2( mh ~ lh )) 



fo(lh) 



0 ) 



+ E h2 

lh€ 7 h2 u 7 h3 U r 23 



E Li( mh - lh )\ 

E hk2( mh - lh )J 



0 

fi(lh)J 



for k = 1, 2. In the last formula inner corners are counted in the union of special 
parts of 7 ^ only once. In addition the special outer corners 

r i4 = {lh G R 2 h \ (G h U 7 ^) : (/ + k 4 )h <E 7 ^ and (l + k x )h € 7 ^} 



and 



r 23 — {lh (Gh U ) : (l + ks)h G 7 h2 and (l -f k%)h G 7^3} 

are considered where f(mh ) is defined by zero. Therefore, these points are only 
important if discrete tangential derivatives are used. 

Lemma 2. For an arbitrary function f it holds that 

E h,M (Thifo, fi})(mh) = f(mh) for all mheG h . 
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The proof of this lemma is published in [12]. We remark that we can find a 
similar structure for T^[/o, fi\(mh) such that 

Dl,M (Thlfo, fi])(mh) = f(mh). 

Based on Lemma 2 a discrete Borel-Pompeiu formula can be established. We 

describe normal unit vectors by using the homeomorphism between complex 

i , u • /a -b\ J . (n\ n\\ _ . 

numbers a + ib and the matrices K I and write I • { I on y hi , i — 

1 , ... ,4 with 



n\ n\ 



n\ n\ 



nr Ur 



Tin n ; 



n i n 2 
nl n\ 



TU 71, 



ni n. 



-1 0 
0 -1 

0 1 

-1 0 



and 



Lemma 3. Let k = 1,2. Then we have a discrete Borel-Pompeiu formula 

{Thk\PZ 2 h + DUu -D^fo + Dlmrnh) 

+ A[f 0 ,mmh) = fk{mh ) 



with 



fi{mh) = ■ 


(fo{mh), 


mh G (G h Uj hl U7J; 


U 


mh <£ (Gh U 7^ U 7^ 2 ); 


f 2 (mh) = 


jfi{mh), 


mh G (Gh U 7^3 U 7 m); 


0, 


m/i £ (G h U 7^3 U 7^4); 



and 



(FL[fo,fi})(mh) = (F^ [fo, fi])(mh) + (F^‘ [/o,/i])M0, 






where 



F hk h [fo,fi])(mh) 

= l E E h 

j=1 ihe^ 



E hk\ (mh-lh) ' 
E hk2( mh ~ lh ) T 



j j 

n\ n J 2 

i j 

n 3 n 4 



fo(lh)\ 

Mlh)J 
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(^ 7 M/o,/i])(m/i) = 



-i Y, h 



E hki( mh ~ lh ) 



'^YS Ei “ (mh ~ ,h> 



E hk2( mh - lh )J \ n 3 + n 3 n\ + n 2 4 



-i Y h 

;/ l e 7 - 1 n 7 - 2 



‘ ^ - \ E M^-lh) 

/he 7h3 n 7h4 



n} + nj n\ + n£\ / fo(lh) 
n 3+ n 3 n\ + nj) \fi{lh) 

n\ + n\ nl + nl\ / f 0 (lh) 

n 3 + n 3 n 4 + n 4/\/l(^) 

„l , „2 „1 I „2\ / r //M 



Elkiimh - lh)\ (n\+n\ n\ + n|\ / f 0 (lh 






n 4 +nj n\ + n\ 
/in + ni n\ -T n\ 



This discrete Borel-Pompeiu formula was already published in [9]. Obviously, 
we get a discrete analogue of the Cauchy integral formula if 

D h,Mf( mh ) = (°>°) 

for all mh G Gh- In a similar way a discrete Borel-Pompeiu formula with respect 
to M can be obtained. 

We define finite difference approximations to the gradient and divergence as 
follows: 



[grad* w 0 \(,nh) = • 

di v^s(mh) = D^ 1 so(mh) + D^ 2 s\{mh). 

Lemma 4. For wq,s G ic 2 1 (G^ 1 ) with wo(rh) = 0 on 7^3 U 7^, so(rh) = 0 on 
7^, and s\(rh) — 0 on 7^ we have 

< In grad^ icq, s > = — ^ h 2 wo(mh) div^s(mh). 

mh£Gh 

Proof. Using the boundary conditions for w Q and s we obtain in analogy to for- 
mula (2.1): 

< /jvgrad^ w 0 ,s > = E /i 2 [(Z)^ wo(mh))so(mh) 



+ (D 2 h w 0 {mh))si(mh)\ 



/1 2 tco (m/i) div^ 5 (m/i). 
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In the following we denote by and Q £ the orthoprojectors onto 

ker D^ M (G h ) and D^ M (w^(G h )), 

respectively, that are defined by the orthogonal decomposition of the space h{Gh) 
from Theorem 1. These projectors can be represented by N x N matrices. 

Lemma 5. The orthoprojectors and Q £ have the properties 

D h,M( P h[ w O’ w i])( m/l ) = (°>°), Vm/i e G h , 

and 

Di,M^h\ w o,wi}){mh) = D l hM ((I N - P+)[w 0 ,n; 1 ])(m/i) 

= D{ M w(mh). 

The proof immediately follows from the definition of the orthoprojectors. 

3.3 Discrete Stokes problem 

As a discrete version of the Stokes problem we consider the boundary value prob- 
lem 

-A h u 0 (mh) + - D\p(mh) = - fo{mh), V mh G Gh , 

p p 

—Ah ui(mh) + — DhP(mh) = — fi(mh), Vra/i G G^, 

p p 

D~f l uo(mh) + Dfr 2 ui(mh) =0, V mh G G^, 

u(rh) = (u 0 (r/i),wi(r/i)) = (0,0), Vr/iG 7 ^, 

where £> denotes the density, p the viscosity, p the pressure, /o and /1 the vector 
components of the external forces and uo and u\ the components of the velocity 
of the fluid. To prepare the proof of an existence result for the solutions of this 
discrete boundary value problem we prove at first some decompositions. 

Lemma 6. For functions 

u(mh ) = (uo(mh),Ui(mh)) G w^Gh) fl ker div^ 
andp(mh) G h(Gh) we have < In Q^[0,p] > = 0. 

Proof From Lemma 1, Lemma 5 and the property div^ u(mh) = 0 we conclude 

< !n Dl' M u,Ql\Q,p\ > 

— < In Dfrjtfii, (In — > 

= h 2 (div^ u(mh))p(mh)+ < u,In D\,m P hi Q iP] > = °- 

mh£Gh 
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The proof of the next theorem is related to Lemma 6 . 

Theorem 2. For each f(mh) = (fo(mh), fi(mh)) E h{Gh) functions 

u(mh) = (uo(mh),ui(mh)) E w^Gh) Hker div^ and p(mh) E h{Gk) 
exist such that 

-(Qh T h[fo,fi])(mh) = —D\ M u{mh) + -(Q^[0,p})(mh), (3.1) 

/i /i 

where is the right inverse operator to the operator In D l h M . 

Proof. Obviously D\^ M u{mh) E imQ^ and (Q^[0,p])(m/i) E imQ^. In 
Lemma 6 the orthogonality of both expressions with respect to the scalar product 
is proved. A further orthogonal subspace can not exist, if from 

< -I N D\ m u , Q+ Tftfo, /i] > = 0 (3.2) 

for all u(mh) E w\(Gh) n ker div^ and 

< Qh[°MiQh T hlfoJi] > = ° ( 3 - 3 ) 

for all p(mh ) E h(Gh) it follows that 

f(mh) = (fo(mh)Ji(mh)) = (0,0). 

At first we obtain from Lemma 1 and Lemma 5 

< —In D\ m u , Q £ T^[/o, fi\ > = < u, In D\,m Qh [/o? fi\ > 

= < u,I N Dl <M Tl[fo,fi] > 

= <u,f>. 



Based on the property ker div^ = (irngrad^)- 1 , from Lemma 4, we can find 
a function go(mh) with go(rh) = 0 on 7^3 U 7 ^, such that from fo(mh) = 
D\go(mh) and/i(m/i) = go(mh) for all mh E Gh and u(mh) E ^(G/Jfl 
ker div^ it follows that 

< u, f > = < u, In grad£ g 0 > = - ^ h 2 (div^n(m/i)) go{mh) = 0 . 

TTlH^G h 

We assumed that equation (3.3) is fulfilled for any function p(mh) E h{Gk) and 
consider now the special case p(mh) = go(mh). Based on the Borel-Pompeiu 
formula from Lemma 3 we obtain 



0= <Qh[ Q ’9o],QhTh[fo,fi] > 

= < Qt[°’9o],QhThlN Dl M [0,g o ] > 

= <Qi[0, so] - Qt [0. so] - Qt Fk [0, So] > 



(3.4) 
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with 

(FtlO,g 0 })(mh) = ((F^[0,go])(mh),(F^9o])(mh)). 

The operator F £ maps the boundary values to ker D\ M and due to the orthogo- 
nal decomposition of the space h(Gh) in Theorem 1 we obtain Q £ F^[0,go] = 
I N (0, 0). Finally, from equation (3.4) it follows that 

WQh [°)5o]||/ 2 (Gh) = 0. 

Based on this result and the properties of the orthoprojectors from Lemma 5 we 
conclude 



f(mh) = [gradj g 0 \(mh) = {D l hM [Q,go}){mh) 

= D h,M(Qh{°’9o}){mh) = ( 0 , 0 ). 



Applying the operator T% in equation (3.1) we get, with the help of the Borel- 
Pompeiu formula, the representation 

- ( n Qt T h\f o, /i])M0 = -Tl I N D\ M u{mh) + - (T% Qt\0,p})(mh) 

/i p 

= - u(mh ) + - {T%Q%[0,p])(mh). 

T 

Using the operators m and D% M and their property 

D h,M = ^2 

we can show now that (u,p) is a solution of the discrete Stokes problem. In the 
next theorem we prove the uniqueness of the solution. 

Theorem 3 . The boundary value problem 

- I 2 A h u(mh ) + - [grad£ p\(mh) = - f(mh), V mh E G h , 

p p 

D^ l uo(mh) + D^ 2 ui(mh) = 0, V mh E G^, 

u(rh) = (u 0 (rh),ui(rh)) = (0,0), Vr/i E 7^ 

has for each right-hand side f (mh) = (fo(mh), f\(mh)) E h{Gh) a unique 
solution u(mh) = (uo(mh),ui(mh)). The pressure 

p(mh) E l 2 (G h U 7^3 U 7^) 



is unique up to a constant . 
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Proof From Lemma 6 and Theorem 2 it follows that 

= II hDl M uf hiah) + i||Q;[0,p]||? l(Oli) . 
p p 



(3.5) 



We assume that two solutions u l (mh) = (u^mh), u\(mh)) and u 2 (mh) = 
(' u 2 (mh ), u\(mh )) exist. The difference 

u*(mh) = u l (mh) — u 2 (mh) 



is a solution of the above boundary value problem with f(mh ) = ( 0 , 0 ) in the 
right-hand side. Without any restriction we can define f(rh ) = (0,0) for all 
rh e 7 ^. Then the left-hand side of equation (3.5) is zero and we can follow 
InDI M u*(mh) — (0,0). Applying the right inverse operator T 2 and the dis- 
crete Borel-Pompeiu formula we get 



u*(mh) — ( 0 , 0 ) 

because F 2 [uq, u\] = I^( 0, 0). We assume now that two solutions 
{u(mh)\pi(mh)} and {u(mh);p 2 {mh)} 



of the boundary value problem exist. Then for 

p*(mh) — p\{mh) — p 2 (mh) 

it follows from equation (3.5) that Q£[0,p*] = /yv(0, 0). 

That means (0 ,p*(mh)) is in the kernel of the operator and therefore in 
the kernel of the operator 1^ D ^ M . Consequently we have 

D] l p^{ r mh) — 0 and D 2 p*(mh) = 0 

for all mh G Gh and it then follows that p* (mh) = const. □ 



We consider now the discrete boundary value problem 



—I 2 Ah u(mh) + - [grad^ p](mh) 
P 



P 



D h l uo(mh) -f D h 2 ui(mh) = ip(mh), V mh G Gh, (3.6) 



u(rh) = x/j(rh) = (^o(rh), (rh)), Vr/iG y h , 

where p(mh) is a measure for the compressibility of the fluid and ^(rh) describes 
the velocity on the boundary. In the case ^(rh) = (0, 0) we have adhesion. This 
system is only solvable if the necessary condition 



^ p(mh)h 2 + ^ D h 1 uo(rh)h 2 + D h 2 u\(rh)h 

mh^Gh r h (=^~ 4 

4 



= £ E ( Mrh) )h 



i= 1 



rh ^l h , 



<*[) \^i( r h) 
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OLc 



is fulfilled. In this formula ^ a i J are the normal unit vectors on 7 ^ , i = 1 , . . . , 4 . 
In detail we have 



OLc 



OL 1 



OLc 



Ot\ 



and — 



OLc 



OLl 



OLc 



a. 



This necessary condition results from 

Y, ip(mh)h 2 + y Df 1 u 0 (rh)h 2 + ^ Df 2 ui(rh)h 2 

mheG h rh ^h3 rk ^h4 

= y uo(mh)h 2 + y D^ 2 iii(m/i)/i 2 

mh£G h U 7 ~ 3 m/iGG h U 7 ” 4 

= y uo(mh)h — y uo(mh)h 

mh£G h \J'y~ 3 mhtGhUy^ 

+ y ui(mh)h — y ui(m/i)/i 

m/i€G h U7~ 4 m/i€G/ l U7~ 2 

V* AW u °( m/l )\, v- /— A fuo(mh)\ 

^ \o) \ui(mh) ) L ( 0 ) \u\{mh) ) 

m/iG7 h3 




Theorem 4. 77i£ boundary value problem 

—I 2 A h u(mh) + — [grad^ p](mh) = — /(m/i), V mh G G^, 

p p 

Df l u$(mh) + Df 2 ui(mh) = Vra/i G G^, (3-7) 

(uo(rh),ui(rh)) = (T^[0, v?])(Wi), Vr/i G 7 ^, 

has for f(mh) = (fo(mh), f\(mh)) G I 2 (Gh) and (p(mh) G h (Gh) with 
(p(rh) = 0 for all rh G 7 ^ a unique solution {u(mh);p(mh)} in the sense 
of Theorem 3 . 

Proof The uniqueness of the solution in the case <p{mh) = 0 was already proved 
in Theorem 3. We denote this solution with {uh(mh)\ph(mh)} and describe the 
solution of the problem (3.7) by the addition of a special solution. Using the ansatz 

u(mh) = Uh{mh ) + (T 2 [0, </>])(m/i), Vm/i G G/j U 7 ^ , 
p(mh ) = ph(mh) + V mh G G^ 
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it follows that 
- I 2 Ahu(mh) 

+ -[g ra d£ p}(mh ) - I 2 & h u h (mh) - D\ M D 2 hM (T 2 [ 0, <p])(mh) 

+ — [grad^ Ph}(mh) + [grad^ ip](mh) 
t* 

= —f(mh) - (Dl M [0,<p])(mh) + [grad^ <p](mh) = -f(mh). 

H [i 

Because T% is right inverse to InD\ m the equalities 

di v^u(mh) = div^(T^[0, ip])(mh) = ip(mh) 

mdu(rh ) = (T^[0, <p])(rh) are satisfied. 

In the following the right-hand side f(mh ) and the boundary values ip(rh) of 
the system (3.6) are given. We construct a solution of the problem using the ansatz 

u(mh) = Uh(mh) + u s (mh), p(mh) = ph{mh); 

here {uh(mh)]ph{mh)} denotes the unique solution of the discrete boundary 
value problem considered in Theorem 3. With this ansatz it follows from the first 
and last equation of the problem (3.6) that the components u s o(mh ) and u s i(mh) 
of u s (mh) are solutions of the Dirichlet problems 

-A h u si (mh ) = 0 , VmheG h , 

u si (rh ) = 'ip z (rh ) , V rh € 7^, 2 = 0,1. 

A uniqueness theorem for the solution of Dirichlet problems is proved in [13]. 
We describe the solution of u s o(mh) and u s i(mh) using our operator calculus. 
From Lemma 5 we obtain 

Dl t M u s( rnh ) = ( Ph In D 2 h M u s ){mh) 

because u s is discrete harmonic. Applying the right inverse operator T% we get 
by the help of the Borel-Pompeiu formula the representation 

u s (mh) = (F 2 [iI>q, ipi])(mh) + (7 P+ I N Df ltM u a ){mh). 

Based on 

D h,M { P h['Po,i>i])(mh) = (0,0) 

and 

D h,M ( T h p h !n Dh,M u *)( mh ) = ( P h !n D 2 hM u s )(mh), 
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the right-hand side of the second equation in (3.6) is determined by 

< p(mh ) = Dpu s o(mh) + D^ 2 u s i(mh) 

= div ^ ( F hli> o , ( T h p h ! n D 2 hM u s )(mh) 

= (Ph2 J N D h,M u s)( m h), 

where P^ 2 denotes the second component of the orthoprojector at each mesh 
point mh. As a summary of these results we state the following theorem: 

Theorem 5. Let u s o(mh) and u s \ (mh) be the solutions of the Dirichlet problems 

-A h u si (mh ) = 0, V mh G G h , 

u si (rh) — f>i(rh), Vr/iG 7 ^, i = 0, 1. 

The boundary value problem 

-I 2 A h u(mh) + — [grad^ p](mh) = - f(mh ), Vm/i G Gh, 

p p 

Dpu 0 (mh) + D^ 2 ui(mh) = (P£ 2 I N Dl M u s )(m/i), Vra/i G Gh, 
u(rh) — f>(rh), Vr/i G 7 ^ 

has a unique solution in the sense of Theorem 3 and can be described by 

u(mh) = u h (mh ) 4- (F 2 [f 0 , })(mh) 4- (T 2 P£ I N D 2 hM u s )(mh ), 

p(mh) = p h (mh). 



3.4 Discrete Navier-Stokes equations 

Based on the Stokes problem we will present a possibility to solve the discrete 
Navier-Stokes equations 



- A hUo(mh) + -Dhp(mh) 



+ - ^uo(mh)D h 1 u 0 (mh) + ui(mh)D h 2 u 0 (mh) j = ^ / 0 (ra/i), 
1 



A h ui(mh) + -Dlp(mh) 



+ - yuo(mh)D h 1 ui(mh) + ui(mh)D h 2 ui(mh) j = ^/i(rah), 
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D h 1 uo(mh) + D h 2 ui(mh ) = 0 , 

u(rh) = (0,0) (4.1) 

for all mh G Gh and rh e 7 ^. Therefore the problem is rewritten. To simplify 
the expressions we use the notation 

Mh,i(mh) — — ( u$(mh)D^ 1 UQ^mh) + ui(mh)Dj 2 uo(mh) ) — — /o(m/i), 
Mh, 2 {mh) = — ( uo(mh)D^ l ui(mh) + ui(mh)D^ 2 ui(mh) ) — — f\(mh). 

mV / m 

Theorem 6. 77ze boundary value problem (4.1) is equivalent to the problem 



u(mh) = (TZQtn[M h , u M h , 2 })(mh) + -(7 Q + h [0,p))(mh) 

M 

- (Q+ = -(<##, p])(m/0 (4.2) 

M 

vv/zere denotes the second component ofQ\ at each mesh point mh G G^. 



Proo/ Let {u, p} be a solution of the problem (4.2). It follows that 
- I 2 A hu(mh) 



— — D 1 D 2 

- ^/i.M u h,M 



1 



{Th Qh Th[Mh,u Mh,2})(mh) + -(T^ [ 0 ,p])(ra/i) 

M 






1 



(Q+T h 1 [M,, 1 ,M M ])(mM) + -(g+[0,p])(mM) 

M 

1 



= - Mh(mh ) [grad^ p] (m/i) 

M 



with Mh(mh) — ( Mh,i(mh ), M^^ra/i)). In addition we have 



div-u(m/i) = div-[(T 2 Q+T, 1 [M, il) M,, 2 ])(m/ l ) + ±(T h 2 Q+[0,p])(mfc)] 

M 

= (Q+7[M M ,M M ])(m/i) + i(Q+ [0,p])(m/ l ) = 0. 

M 

The boundary condition of problem (4.1) is satisfied because, based on the prop- 
erty of the orthoprojector Q the solution of the problem (4.2) has the represen- 
tation 



{T% In D 2 h M s)(mh) = s(mh) - (F%[s 0 , si])(mh) 
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with s(mh ) G w^Gh)- Let now {u,p} be a solution of the problem (4.1). With 
the help of the Borel-Pompeiu formula we obtain 

(Tt I N A h u)(mh) = {T l h I N Dl M Dl M u)(mh) 

= D h,M U ( mh ) “ ( F h{ D l,M U \)( mh )• 

Using the orthoprojector we get 

(Q+ T l h I N A h u)(mh) = (Q+ I N Dl M u)(mh) - (Q+ 

= D 2 hM u{mh). 

Otherwise it follows from (4.1) that 

(Q+T£l N A h u){mh) = (Q+T£[M h , u M h , 2 ])(mh) 

+ ~(QtThI N D\ M [0,p})(mh) 

= {QlTl\M h<u M h Mmh) + -(QtMimh). 
From both equations we obtain 

Dl M u{mh) = (Q+T h 1 [M Ml M fc , 2 ])H + -(Q£[0,p])(mh). 

Applying the operator T% and again the discrete Borel-Pompeiu formula, the first 
equation of (4.2) follows. The second equation results from 

0 = di v^u(mh) 

= di ^[{TlQtniMh^MbMmh) + Q+[0,p])(mft)] 

r 

= (g, + 2 T‘[M M ,M M ])(m/i) + -(Qj 2 [0,p])(mft). 

A 4 

Based on the problem (4.2) the following iteration procedure can be used to 
calculate the solution of the Navier-Stokes equations (4.1): 

u n (mh) = (TjtQ+THiM^^Dimh) + ±(T h 2 Q^[0,p n })(mh), 

- {Q + h2 Tl[M^\M n h ^]){mh) = ±(Q+ 2 [0,p n })(mh) 

r 

with 

M h7j\ mh ) = ^(uZ-\mh)D^u]ll(mh) + u^-\mh)D^:l(mh)^ 

- -fj-iimh) 
A 4 
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for j = 1,2, 

(uQ(mh),Ui(mh)) £ w 2 (Gh) Hker div^ 
and n = 1, 2 , 3, — 

In this way the solution of the boundary value problem (4.1) is reduced to the 
solution (u n ,p n ) of a Stokes problem for n — 1,2,3, .. . and the uniqueness 
of this Stokes problem is proved in Theorem 3. We have to prove that the itera- 
tion procedure converges to the solution of the problem (4.1). The following two 
lemmas prepare this proof. 

Lemma 7. For 1 < p < 2 and q < or, p — 2 and q < oo, the operator 
T 2 : l p (Gh ) — > l q (Gh ) is bounded. 

Proof. Let p < q. At first we estimate the expression \(T 2 k [fo, fi])(rnh)\ for 
k = 1,2. Because T 2 acts only on functions that are defined in Gh we set f(rh) = 
(0, 0) for all rh £ 7 ^ . Based on the property 

| Elkjimh -Ih ) | < C\(\mh - lh\ + h)~ l with fc, j £ { 1 , 2 }, 

(see [7, 21]) we get 



ICr?,*[/o,/i])(rok)| 



lh£Gh 



El kl (mh-lh)\ T 

E 2 hk2 (mh-lh)J 



fo(lh)\ 

Mlh)J 



<C X Y, h\\mh-lh\ + h)-\\f 0 (lh)\ + \f x {lh)\) 

IhtzGh 



<2 C x h 2 (\mh-lh\ + h)~ l \f(lh)\. 

lh£G h 



For 1 < p < 2 it follows - > 1 — % with -L -f - = 1 and therefore, 1 = 

r — q p* p* p ’ 

| — 2/3 with /3 > 0. We get the equation 

^ h 2 (\mh — lh\ + h)~ l \f(lh)\ 
lh£Gh 

= £ h? (\f (lh)\* (\mh - lh\ 4 - 

iheG h ' J \ J 

x ({\mh — lh\ + h)~ 



1 1 

a i = -, a 2 — - 

q p 



-, and a 3 = — , 

q p 



We set 
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such that <*i + a 2 + a 3 = 1 and apply the Holder inequality to the product of the 
three functions. This leads to the estimate 



h 2 (\mh — lh\ + h)~ 1 \f(lh)\ < 

lheG h 



^h^filhWilmh-lhl+hyWo 

IhtzGh ' 



X 



J2h 2 \f(lh)\ p 

ih£G h 



^Ph 2 (\mh — lh\ + h) 2r ^ p 
lh£G h 



1 _ 1 

While we have for the second factor ( ^ h 2 \f(lh)\ p ) = 1 1 / 1 1 ^ ^ h y we 

\lh£Gh ' 

can estimate the third factor in the form 



lheG h 



/3p* 



In order to simplify the notation we define 



K( * , r S r (diam G h )^ (diam G h )^ +1 ’r 
h{p Gh) — C 2 — — C 2 



(3p* 



El + 1 _ £1 

q 2 



and obtain 
\(Th, k {foJi})(mh)\ 

<2C 1 (K(p\q,G h ))^\\f\\^J J2h 2 \fm\ p (\mh-lh\+h)-^y . 

lh£ Gh 

Using the Holder inequality again, we conclude 

ra q lq(Gh) = e h^Ti^+n^i 

mh£Gh 

< {2V2C 1 Y{K(p\q,G h ))^\\f\\l P Gh) Y,h 2 \f(lhWK(q,p*,G h ). 

lh£G h 



Finally we obtain 



\\ T h\\l q (G H ) 



<2 y/2Ci (K(p*,q,G h ))*(K(qy,G h ))'\ 



lUp(Gh)- 



We remark that the operator T £ has the same property. Important are two spe- 
cial cases, which immediately follow from Lemma 7. 

Corollary 1. The operator T% : l p (Gh) — > h{Gh) is bounded for 1 < p < 2, 
while the operator T 2 : hiGh) — > l q (Gh) is bounded for all q < oo. 
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In the next lemma we prove the boundedness of in the Sobolev norm 



IMLi(Gh) - IMktG,,) 



£ 

2=1 



\W D h w °IIj 2 (GhU 7 - 3 u 7 - 4 ) + W D h t 'iH; 2 (G h u 7 - 3 u 7 - ) 



1 

2 



Lemma 8. The operator T 2 : hiGh) n im Q\ — > w\(Gh) is bounded. 

Proof. For each 

u(mh) £ h{G h ) n im Q £ 

there exists a function v(mh) £ w\{Gh) with u(mh) = D\ M v(mh). From the 
Borel-Pompeiu formula we get v(mh) = T%u(mh). Based on equation (2.1), 
the first Green’s formula, which is proved in [13], and Poincare’s inequality we 
obtain 



IMI? 2 (G/i) ~ II ^h,M v \\ h{G h ) ~ < In In Dh,m v > 

— < —In D}i,m Dh,M v ' v > 

= ^2 h 2 [(-A h vo(mh)) vo(mh) + (-A h vi (mh))vi(mh)\ 

mh£Gh 



E E ^ [( D h l M mh )) 2 + (D h l vi ( mh)) 2 



>c 3 - x 



h(G h ) 



+ 



h(G h 



= c. 



-1 



I h(G h 



where the constant C3 is related to the smallest eigenvalue of the discrete Laplace 
operator. From this estimate we conclude 



\Th u \\w\(G h ) ~ IML^G,,) 

= IMIi2(G h ) 

2 



+ Eii^rv 12 



i= 1 



i 2 (G h U 7 h 3 U 7 - 



+ ll^lllf 2(C . 






< (C3 + !) ll u ll/ 2 (G h )- 



Theorem 7. Let 



l v (G h ) 



< c 4 



( 





16KC4Q 



2 

,1 < p < 2. For any function 



u°(mh) £ w\(Gh) H ker div h 
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with 



ll u ILj(Cf h ) ^ 






I6KC4Q 



+ n and = 






16 KC 4 q 



II l P (G h ) 



c 4 



the above iteration procedure converges to the solution of the Navier-Stokes equa- 
tions (4.1) . 



Proof Let [grad h uf](mh) = 



D h l u'>(mh) 
Dl 2 vJi(mh) 



, for i = 0, 1 and 



^ \\Th\\ [l 2 (G h )r\imQ+,w\(G h )] W^hW [lp(Gh),h(Gh)] 

based on Lemma 7 and Lemma 8. In the following we consider the difference 

u n (mh ) — u n ~ l (mh). 

Replacing (fo(mh), fi(mh)) by (Mh,i(mh), we obtain in anal- 

ogy to (3.5) 



t 1 



n— 1] 



II h(G h )nimQ+ 



< II - M”7 2 , M^ 1 - M«- 



II h(Gh)C\ im Q 



Using the property \\Qh\\{i 2 (G h )MG h )nimQ+} = 1 of the orthoprojector Q+, it 
follows from the iteration procedure that 



IK-« n 

n2 / o+T’lrA/f n - 1 



< \\ltQh T h[Kj - % Ka - ]Lj(o h) 

+ i||T 2 Q+[0,p"-p”- 1 ]|L, (G , ) 

r 

- 2 ll^ll[; 2 (G ) ,)nimQ+, iuJ(G h )] IIQfc H[t2(Gh),t2(G h )n»mQ+]ll-^ll[« P (G ! h),«2(Gh)] 

x \\M n h ~ l - - M n h f || lp(Gh) 

<2K\\Mlf - - M^ 2 ||, p(Cfc) . 

Using the inequality o 2 + b 2 < (|a| + |6|) 2 and the Minkowski inequality we 
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obtain 



- 

m v 



p \ 





E 


(u 11 krad^UQ 1 - u n 2 grad h v% 2 ) 2 


mh£G h 


L 




+ {u n ~ l grad* «" _1 - u n- 2 grad* u " -2 ) 2 


E * 


|u" _ 1 grad* u™- 1 _ u " _ 2 grad^ u%~ 2 \ 



^ \ ? 



mh£Gh 



+ | u n 1 grad /l < 1 — u n 2 grad /l u™ 2 | 



< — \\u n 1 grad^ u™ 1 - u n 2 grad h u n - 



n — 2 1 



^J =0 



lip(G fc ) 



1 r 






■«- 1 - U "- 2 )grad,- W "- 1 || ip(Gh) 



+ ll u " -2 grad^ (^- 1 - Uj "- 2 )|| lp(Gh) 



Both norms on the right-hand side have the structure || w grad h v 3 \\i (g h ) • From 
Holder’s inequality it follows that 

Ik grad* v j\\i P (G h ) < Iklli^^Gh) II grad^ vj ||/ PP1 (Gh) 

with — + — = 1. If we choose pp\ = 2, then we obtain q := pq\ = and 
p = 7 ^ . Based on the definition of the discrete Sobolev norm we get for j = 0,1 
the estimate 



Ik grad ft Vj\\ lp(Gh ) < |klli,(G h )( ^KAiS) 2 + ( D h \ ) 2 ] 

771 /iG G h 
' 2=1 

< lklll,(G h ) IML>(G h )- 

For functions w G w^Gh) we show now \\w\\i q (G h ) < C 4 H^H^i (c?^) • Because 
of the special choice of the boundary values, we obtain from the Borel-Pompeiu 
formula |kll/,(G h ) = \\ T h D l,M w \\i q (G h ) with 



D 2 hM w g l 2 (G h ) n im Q%. 
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From Corollary 1 and the equivalent formulations in the proof of Lemma 8 we 
can conclude 

\\T h Dh,M W \\l q (G h ) ^ \\ D h,M W Wh(G h ) < C4 |Hlu>J(G h )- 
On these estimates the result 



\\u n -U n 



<AKC 4 -\\u n ~ l -u n ~ 2 
is based and with the notation 



+ l|w n 2 IL*(g,,) 



L n = 4 KCA ||w" _1 |Ui ( g,) + IK- 2 !!^)) 

r 



we obtain 



U--U- ~W w l { G h )<L n \\u n 2 Li(G h )- 



i n -u n ~ l 
We assume now that the inequality 



0 < ||« n-1 ILi ( G fc ) < x&kcZq + n 



holds with 



ith f] = ( i6KC 4 ^) — Under the assumption 






\\l p (G h ) <C 4 \ l6RC4g 
it follows from the iteration procedure, in analogy to the above steps, that 

^2 I 



! U 



.;<»,> <2K;( E 

^ V mheG h 



(u n 1 grad^ Uq 1 - / 0 ) 2 

+ ( u n ~ l grad^ - /i) 2 



E x 1 
2 \ P 



1 r 



<2 K*y 



IK 1 grad/j u " + ||/||i p (G h ) 



j=0 L 
Q_ n ,n— 1 1 1 2 



< wc^\\u n -r wl(Gh) + ^\\f\\ lp{Gh) 



< 4 KC 4 - 



(_JL 



LV16 KC 4 Q 



+ ^ ) ~f 



\\l P (G h ) 

c 4 



M \ K , /rn 2 , ll/I MGh) 



Kisses J + 8 ^ + (fi) + ft 

< M I n , 1frr K m V < I ( __E__ I Q 
-64E‘C 4£ . + 2 + “2L6^C' 4 £» + 
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Because of the inequality 



Ln + l = 4KC 4 -( ||u n |Ll(G h) + ll u IL'(G h )) 



< 4K Ct 






e 3 

n 2 \I 6 KC 4 Q 



+ S2 









for all n > 1, the assertion of the theorem follows from Banach’s fixed-point 
theorem. 
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A Symmetric Functional Calculus 
for Systems of Operators of Type uj 

Brian Jefferies 



ABSTRACT For a system A = (A i, . . . , A n ) of linear operators whose real lin- 
ear combinations have spectra contained in a fixed sector in C and satisfy resol- 
vent bounds there, functions f(A) of the system A of operators can be formed for 
monogenic functions / having decay at zero and infinity in a corresponding sec- 
tor in R n+1 . In the case that the operators Ai, . . . , A n commute with each other 
and satisfy square function estimates in Hilbert space, the correspondence between 
bounded monogenic functions defined in a sector in R n+1 and bounded holomor- 
phic functions defined in a sector in C n is used to define the functional calculus 
/ /(A) for bounded holomorphic functions / in a sector of C n . The treatment 

includes the Dirac operator on a Lipschitz surface in R n+1 . 

Keywords: functional calculus, plane wave formula, Dirac operator 



4. 1 Introduction 

For a finite system A = (Ai, . . . , A n ) of bounded linear operators acting on a 
Banach space X , functions /(A) of the n-tuple A can be formed via the Cauchy 
integral formula 



f(A) = / G x (A)n(x)f(x)dfjL(x), (1.1) 

Jdn 

just under the assumption that the spectrum cr((A, £)) of the operator 

n 

(AO :=£A0 

3 = 1 

is a subset of R for every £ E R n [6]. The Cauchy kernel x h-» G X (A ) is defined 
outside a distinguished subset 7 (A) of R n , the open set Q c R n+1 contains 7(A) 
and has a smooth oriented boundary dQ with volume measure fi and outward unit 
normal n(x) at each point x of the boundary. The function / is defined in a neigh- 
bourhood of the closure Q in R n+1 , takes values in the Clifford algebra C£(C n ) 
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over C n and is left monogenic in the sense of Clifford analysis. In particular, the 
bounded linear operator /(A) : X — > X is defined for any real analytic function 
/:[/—>€ defined in a neighbourhood U of 7 {A) in R n , simply by replac- 
ing / in formula (1.1) by its left monogenic extension to an open set in R n+1 [6]. 
For a polynomial p in n real variables, p(A) is the operator formed by substituting 
symmetric products in the bounded linear operators A \ , . . . , A n for the monomial 
expressions in p that is, we have a symmetric functional calculus in the n opera- 
tors A\ , . . . , A n . In the case n — 1 , we obtain the usual Riesz-Dunford functional 
calculus for a single operator. By analogy with spectral theory in one operator, the 
set 7 (A) is called the monogenic spectrum of the n-tuple A = (Ai, . . . , A n ). 

The key ingredient of this approach is the Cauchy kernel G X (A). In [6], it is 
defined by the plane wave formula 

= sgnOro )"- 1 



x / (e 0 + is) ((x, s)I — (A, s) — x 0 sl) n ds (1.2) 

J s n ~ l 

for all x = xo&o + x with xq a nonzero real number and x G R n . Here 5 n_1 
is the unit (n — l)-sphere in R n , ds is surface measure and the inverse power 
((xl — A, s) — x 0 s)~ n is taken in the Clifford module C(X)®C( n y The spectral 
reality condition 

cr((A, £)) C R, for all £ E R n , (1.3) 

ensures the invertibility of ((xl — A, 5 ) — xo s) for all xq ^ 0 and s G S' 71-1 by 
the spectral mapping theorem. The Cauchy kernel G X (A) given by formula (1.2) 
coincides with the convergent series expansion of G X (A) for large |x|. 

If we now pass to unbounded operators, then a similar analysis holds if we 
retain the spectral reality condition (1.3), provided that we suitably account for 
operator domains. 

The question of forming functions of noncommuting systems of operators a- 
rises in quantum physics [14] and the connection with Clifford analysis is already 
apparent in [9]. Explicit calculations can be made with two hermitian matrices [7]. 
The application of ideas from Clifford analysis to Feynman’s operational calculus 
is pursued in joint work of the author with G.W. Johnson [5]. 

The examination of functional calculi for n-tuples of operators, and formula 
(1.2) in particular, was also motivated by the study of the commuting n-tuple 
-Ds = (D 1, . . . , D n ) of differentiation operators on a Lipschitz surface E in 
R n+1 , see [13]. In the case that E is just the flat surface R n , the operators 



1 d 



commute with each other and are selfadjoint, otherwise, the unbounded operators 
Dj, j = 1, . . . , n, have spectra cr(Dj) contained in a complex sector 



5 W (C) = {2 G C : z / 0, | arg(z)| < u} 
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with an angle u depending on the variation of the directions normal to the sur- 
face E. 

The existence and properties of the H 30 -functional calculus for the commuting 
n-tuple Dy are now well-understood, see for example [13]. Here we show how 
a functional calculus for an n-tuple A of unbounded operators of type cj can be 
constructed directly and to apply this to the problem of forming bounded linear 
operators /(Dy) acting on L p ( E) when / is a bounded holomorphic function 
defined on a suitable sector in C n . The construction of an H °°- functional calculus 
for Dy has applications to the solution of irregular boundary value problems. The 
spectral reality condition (1.3) needs to be replaced by a condition where the 
spectra are contained in a fixed sector of angle uo in the complex plane in order 
to treat systems of operators such as Dy . The convergence of the integral on the 
right-hand side of the plane wave formula (1.2) for the Cauchy kernel is then at 
issue. 

In Section 4.2, it is shown how formula (1.2) for the Cauchy kernel associated 
with the system A of sectorial operators still works if the spectral reality condition 
(1.3) is replaced by a sectoriality condition with the appropriate resolvent bounds. 
The system Dy of commuting sectorial operators described above is of this type. 
By this means functions f(A) of the operators A can be formed, provided that / 
is left monogenic in a sector in R n+1 and satisfies suitable decay estimates at 0 
and oo, in a fashion similar to the case of a single operator of type [12]. Because 
G X (A) is only defined for x outside a sector in R n+1 , the monogenic spectrum 
7 (A) is now contained in that sector in R n+1 . Recall that under condition (1.3), 
7(A) is a subset of R n = R n x {0}. 

A function /(Dy) of the system Dy has a natural interpretation as a multiplier 
operator acting on L p -spaces of functions defined on the Lipschitz surface S as 
well as a singular convolution operator, so the multiplier / should be a bounded 
holomorphic function defined on a sector in C n [13], rather than a bounded mono- 
genic function defined in a sector in R n+1 . 

In this work, we use the recently proven bijection [3] between bounded mono- 
genic functions defined on a sector in R n+ 1 and bounded holomorphic functions 
defined on a sector in C n to resolve this apparent dilemma about what is the ap- 
propriate function space for the symmetric functional calculus. The association 
between the two function spaces is via the Cauchy-Kowaleski extension to a sec- 
tor in R n+1 of the restriction of the holomorphic function to R n \ {0}. 

The bijection between bounded monogenic functions and bounded holomor- 
phic functions in sectors is described in Section 4.3. In Section 4.4, this enables 
us to form functions f(A) of a commuting system A of operators acting in a 
Hilbert space when the single operator Ajej satisfies square function esti- 
mates and / is bounded and holomorphic in a sector. 

Some notation and facts from Clifford analysis [1,2] follow. 

Let C£( C n ) be the Clifford algebra generated over the field C by the standard 
basis vectors eo, ei, . . . , e n of R n+1 with conjugation u u. The generalized 
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Cauchy-Riemann operator is given by 



D 






d 



' } dxj ' 



Let U C R n+1 be an open set. A function /:[/—» C£(C n ) is called left 
monogenic if Df = 0 in U and right monogenic if f D = 0 in U. The Cauchy 
kernel is given by 

G x (y) = — X V x,y e R n+1 ,x ^ y, (1.4) 

o n F - y | n+1 

with o n — 2 tt ]1 2 1 /r ( n ^) the volume of the unit n-sphere in R n+1 . So, given a 
left monogenic function /:[/—> C£( C n ) defined in an open subset U of R n+l 
and an open subset ft of U such that the closure ft of ft is contained in [/, and 
the boundary dft of ft is a smooth oriented n-manifold, then the Cauchy integral 
formula 

f{y)= f G x (y)n(x)f(x)dn(x), yeti 
Jan 

is valid. Here n(x) is the outward unit normal at x G dft and p is the volume 
measure of the oriented manifold dft. An element x — (x 0 , xi, . . . , x n ) of R n+1 
will often be written as x = x 0 e 0 + x with x = Y^= 1 x j e j • 



4.2 The plane wave decomposition 

Let A — (A \, . . . , A n ) be an n-tuple of densely defined linear operators 

Aj : V(Aj) -> X 

acting in X such that C^ =l V{Aj) is dense in X. The space C^ n ){X^ n )) of left 
module homomorphisms of X^ n ) = X 0 Cl( C n ) is identified with C(X) <g> 
C£( C n ) in the natural way and becomes a right Banach module under the uniform 
operator topology. 

If we take formula (1.2) as the definition of the Cauchy kernel G X (A), then we 
need to establish the convergence of the integral 

/ (eo + is) ((xl — A,s) — xq sl)~ n ds 

Js — 1 

for particular values of x = xoeo + x G M n+1 . Now 

((xl — A , s) — xosiy 1 — ((xl — A, s) + xo si) ((x/ — A, s ) 2 + XqI) 

if 0 ^ cr [(xl — A , s) 2 + x§) . Thus, we need to ensure the appropriate uniform 
operator bounds for 



((xI-A,s) 2 +x 2 0 I) \ seS ”- 1 
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as x — xoeo 4- x ranges over a subset of R n+1 . In the case that a ((A, s)) C R 
and (XI — (A, s)) _1 is suitably bounded for all s G S n ~ l and A G C \ R, then 
G X0 e 0 +x (A) is defined for all xq 0. First, for 0 < v < ^7r, set 

S t ,(C) = {z G C : | arg^| < v} U {0}, S°(C) = {z G C : | argz| < u}. 

The set of s G 5 n_1 with nonzero coordinates Sj for every j = 1, ... , n is 
denoted by Sq~ 1 . Then Sq' 1 is a dense open subset of S' 71-1 with full surface 
measure. 

Definition 1. Let A — (A \, . . . ,A n ) be an n-tuple of closed, densely defined 
linear operators Aj : V(Aj) — > X acting in X such that D^ =1 P( Aj ) is dense in 
X and let 0 < u < ^ . The operator 

(A,s):n] =l V(A J )-*X 



is defined by 



(A, s)u = J2 s Mi u ) . Vs e « e r\ n j=l V{Aj). 

j = l 

Then A is said to be uniformly of type to ifa((A , s)) C S u; (C)for all s G Sq _1 
and for each v > tu, J/iere exists C v > 0 such that 

||(z/- (^s))- 1 !! <C'„|zr 1 , z£S„(C), SG5”- 1 . (2.1) 

Ifn = l,the A is just said to be of type u. 

It follows that s 1 — > (A, s) is continuous on Sq _1 in the sense of strong resol- 
vent convergence [8, Theorem VIII. 1.5]. The subset Sq “ 1 of 5 n_1 is used here 
simply because n^ ==1 D(i4j) may be strictly contained in V(Ak) for k = 1, . . . , n. 

Now suppose that equation (2.1) is satisfied and let z = (x, s) 4- ix 0 . Then 
z £ S°( C) means that |xo| > tan v\(x,s)\. 

First, let 



N u = {xG R n+1 : x = xoeo 4- x, |x 0 | > tanz/|x| }, 

5^(R n+1 ) = {xG R n+1 : x = xoeo 4- x, |xo| < tanz^|x| }, 

5°( R n+1 ) = {xG R n+1 : x = x 0 eo 4- x, |xo| < tanz/|x| }. 

Note that if x 0 eo 4 -xG^, then z = (x, s) 4- zxo £ 5°(C) for every s G 5 n_1 , 

because |xo | >tan^|x| > tan^|(x, s)\. 

The proof of the following result is straightforward and appears in [4]. 

Lemma 1. Let uj < v < I 7 r. Suppose that the n-tuple A of linear operators is 
uniformly of type lu. Then for all xoeo + x E N u , the integral 




- A,s) - xq si) 



ds 

C (n) (x (n) ) 
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converges and satisfies the bound 



£(n)(X( n )) Fo| n 

Thus, if A is uniformly of type u, then 

xoeo + G Xoeo +$(A) (2.2) 

is defined by the plane wave formula (1.2) for all x 0 eo + x € N„ with u < v < 
^ 7 r. Standard arguments ensure that (2.2) is both left and right monogenic as an 
element of C^ (X^). If we denote by 7 (A) C R n+1 the set of all singularities 
of function (2.2), then 7 (A) C 5^(M n+1 ). 

Suppose that cj<i^<^7r,0<s<n and / is a left monogenic function 
defined on S°( IR n+1 ) such that for every 0 < z/ < 0 < v there exists Coy > 0 
such that 

\f(x)\ < Cey a: £ S|(R n+1 ) n N„, (2.3) 

According to Lemma 1, for every u < z/ < 9 < v, we have 

l]^r)ll.|/(x)l<er | Xo |n ( ^ N 2 s) ' x = x 0 e 0 +x 

for all x e 5^(R n+1 ) fl N v >. Now if lj < 6 < v and 

He = {x £ R n+1 : x = x Q e 0 + x, |x 0 |/M = tan0} C S°(R n+1 ) (2.4) 

it follows that 

\\G x (A)\\.\f(x)\ = 0(l/\x\ n ~ s ) as in H e . 

Hence, x i-> G x (A)n(x)f(x) is locally integrable at zero with respect to n- 
dimensional surface measure ^ on Hg. Similarly, 

\\G x (A)\\.\m\=0(l/\x\ n+s ) as |x| -> oo in H e , 

so x i — > G x (A)n(x)f(x) is integrable with respect to n-dimensional surface mea- 
sure on Hq . Therefore, we define the element f(A) of the module C^(X^) by 
the formula 

f{A) = f G x (A)n(x)f(x)dp(x). (2.5) 

J He 

If ^ : ® n \ {0} — > C has a two-sided monogenic extension ip to 5°(R n+1 ) that 
satisfies the bound (2.3) for all 0 < 0 < v, then ip (A) is written just as ip(A). 
Formula (2.5) does just what we would expect in the noncommuting situation. For 
example, let p be a polynomial of degree n with 

p(0) = 0 and b\(z) = p{z)(X — z)~ n ~ l 
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for some A ^ S°( C). Let £ E R n and set 

fa,z(x) = h{{x,0) 

for all x G K n . Denote the two-sided monogenic extension of (j)\^ to S°(M" +1 ) 
by Then <f)\^ has decay at zero and infinity and we have 

fa, ((A) = fa,t(A)=p((A,Q)(\I - (A,0)- n ~ l 

is a bounded linear operator. More generally, the next result shows that formula 
(2.5) gives the right result when for a special class of functions /, there is a repre- 
sentation of the bounded linear operator f(A) by an alternative formula, namely, 
the usual Riesz-Dunford functional calculus. Set D r = {z E C : \z\ < r} for 
each r > 0. 

Theorem 1. Let 0 < uo < v < ^7r and suppose that A = (A \ , . . . , A n ) is an n- 
tuple of operators uniformly of type uj with the property that for each lu < tu f < v, 

there exist n-tuples A^ = \ . . . , A^^j , fc = 1,2,..., of bounded linear 

operators such that, with uj' replacing to, the bound (2.1) is satisfied uniformly for 
k = 1,2,... and 

{zI-(A^ k \s))- 1 ( zI-(A,s »-\ z i S v ( C), s G S%~\ (2.6) 

as k — > oo. 

Suppose that 0<s<l, 0<u;<i/, r>0 and f : S°( C) U D r — > C is a 
holomorphic function satisfying f(0) = 0, such that for each 0 < i/ < v, there 
exists C v > > 0 with 

1^)1 -TT^F’ zGS ° AC) ■ 

Let £ E R n be a vector with nonzero components. Then the function 
x 1 — > 0((x,O), X E R n \ {0}, 

has a two-sided monogenic extension f to S°( R n+1 ) satisfying the bounds (2.3) 
and the operator f (A) E C{X) given by formula (2.5) has the representation 

f(A) = fa{A,Q) = L. J (XI - (A,®)- 1 fa A) d\. (2.7) 

The contour C can be taken to be {z E C : | Im(2)| = tan#| Re(z)| }, with 
0 < u < 8 < v. 

Proof We can suppose that £ E Sq _1 by scaling. It is routine to check that the 
two-sided monogenic function / is given by the Cauchy integral formula 

f(x) = 2— f (z - (x,0 +x 0 £)~ 1 <p(z)dz, x G S°,(R” +1 ), 

J c 




66 Brian Jefferies 



for all x = xoeo + x £ S°, (R n+1 ) and all 0 < v 1 < 6. Furthermore, as noted in 
[13, Section 5.2], the element i£ of the Clifford algebra C£( C n ) has the spectral 
decomposition 



it = x + m\-x-(m\ 



with respect to the projections 





so that / also has the representation 



f(x) = <£((£,£) + ix 0 |£|)x+(0 + <t>((x,0 -ia;o|£|)x-(0. 

for all x = xoeo + x £ 5°(R n+1 ). Hence, / satisfies the bounds (2.3). Note that 
(x, £) = 0 if x is orthogonal to £, but by assumption in (2.3), |xo| > tan v'\x\ for 
1/ > 0. 

Now suppose that u/ < 6 < v. We first observe that 
[ G x (r^ (fc) )n(x)(C - (x,0 + x 0 ^) _1 d/x(x) = ((I - r(A^ k) (2.8) 

JHo 

for every ( € C \ S t/ ( C) and r £ R. The (improper) integral is over the cone 
He = {xoeo + x £ R n+1 : |xq| = tan0|x| }, 



which can be deformed by Cauchy’s theorem in Clifford analysis to the integral 
over a ball B r of radius r > centred at zero in R n+1 when \(\ is large 

enough. Both sides of equation (2.8) are holomorphic for r in a neighbourhood 
of zero in C and have the same Taylor series about zero, because 



dP_ 

dri 



I G x (rA {k) ) n(x)(C - (x, £) + x 0 0 1 dy,(x) 

JH e 

= [ [(-A {k) ,V s ) J G x (TA ik) )}n(x)(( - (£,£,) + xoty 1 dv(x) 

JH e 

= [ (->l (fc) , Vx> J G x (r J 4 (A:) )n(x)(C - (x,£) + x 0 £) -1 d,u(x) 

J JR n ±ee 0 

= f G x (TA {k) )n(x)(A {k) ,V s ) J (( - (£,£) + x 0 0~ 1 dn{x) 

JR n ±ee o 

— > J as r — > 0. 



Here we have used the observation, apparent from the plane wave formula (1.2), 
that for each k = 1, 2, . . . , the Cauchy kernel (x, r) i— ► G x (rA^) satisfies the 
operator equation 

^-G x {tA^) + (A^,V x )G x (tA^) = 0 
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for all x £ N v . Then the equality (2.8) is obtained for all ( £ C \ S„(C) by 
analytic continuation of both sides. Next, the cone He can be deformed to an n- 
dimensional surface consisting of part of the sphere of radius 0 < r' < r together 
with Hq \ B r > and the contour C can be deformed so as to avoid the disk D r > . 
Then it follows that the equalities 



f{A {k) ) = [ G x (A w )n(x)f(x)d/j,(x) 

J He 



(C — (£) 0 + x oO V(CK 



dfi(x) 



= hL G ^ k ^ 

= X— / [ G x (A {k) )n(x)(C- (x,0 + xo0~ 1 dfi(x) 

Jc lJ Hq 

= <K(aw,q) 



hold. Now the assumption (2.6) of strong resolvent convergence guarantees that 



in the strong operator topology as k — » oo, so the convergence of the integrals 
f(A (k) )=[ G x (A {k) )n(x)f(x)dn(x) 

J He 



to f(A) as k — > oo remains to be established. According to Lemma 1, we can find 
C > 0 such that 

l|G I (^ (fc) )lk (n) (x (n) ) < 

for all x E He and k = 1,2, Moreover, the plane wave formula (1.2) and 

strong resolvent convergence (2.6) ensures that G x (A ^ ) — > G x (A) in the strong 
operator topology for each nonzero x £ He- The monogenic function / has decay 
at zero and infinity, so dominated convergence ensures that f(A^) — * f(A) as 
k —> oo. □ 

The approximation (2.6) by bounded operators is somewhat simple-minded. A 
more sophisticated approximation should yield the result for general systems of 
operators uniformly of type u. In order to form functions f(A) of the system A 
of operators for a class of monogenic functions / larger than those which sat- 
isfy a bound like (2.3), a greater understanding of function theory in the sector 
^(R 71 "^ 1 ) is needed. To this end, the simple system A = ( = ((i , . . . , ( n ) e C n 
of multiplication operators in the algebra C£(C n ) is considered in the next section. 
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4.3 Function theory in sectors 

For the particular example of the n-tuple of differentiation operators on a 
Lipschitz surface E mentioned in the Introduction, the theory described in [10] 
associates a bounded linear operator f(D^) with a Fourier multiplier /, this being 
a suitable holomorphic function of n complex variables defined in a sector in C n 
associated with D £. By means of formula (2.5), we have seen how to associate 
a bounded linear operator f(D^) with a left monogenic function / with suitable 
decay in a sector, once we verify that Dy satisfies the assumptions of Definition 1. 

This suggests that we can associate a bounded holomorphic function of n com- 
plex variables in a sector in C n with a left monogenic function / with suitable 
decay at zero and infinity in a sector in R n+1 . The association is by analytic 
continuation from R n \ {0} to a sector in C n . The precise formulation of this 
correspondence and its consequences is formulated in this section in Theorem 2. 
The proof is somewhat technical and will appear elsewhere [3]. In the next sec- 
tion, Theorem 2 and some additional estimates are used to extend the functional 
calculus defined by means of formula (2.5) to all bounded holomorphic functions 
defined in a sector C n , at least in the case when A is an n-tuple of commuting 
operators. 

Let 0 < v < ^7r and let i/^°(5°(K n+1 )) denote the set of all left monogenic 
functions 

/ : S°(R n+1 ) -+ Cl{ C n ) 

that are uniformly bounded on every subsector S°, (R" +1 ), 0 < //' < //.Endowed 
with the topology of uniform convergence on subsectors S°, (R n+1 ), 0 < z/ < i/, 
the topological vector space H^°(S°(W 1 ^ 1 )) is a Frechet space. The analogous 
space for right monogenic functions is written as R n+1 )). 

For any / G R n+1 )), the restriction 

fo : R n \ {0} -> C£(C n ) 

of / to R n \ {0} is real analytic and so it has a complex analytic extension / 
to some open neighbourhood U v of R n \ {0} in C n \ {0}. The argument above 
concerning Fourier multiplier operators on Lipschitz surfaces suggests that we 
can take U v to be some open sector in C n on which / is uniformly bounded on 
subsectors. The sectors in C n are described as follows. 

For x G R n+1 and n complex variables ( = ((i, . . . , Cn), the Cauchy kernel 
G x (C) is understood to be the maximal analytic continuation of the Cauchy kernel 
G x (y), y G R n , y ^ x, given by formula (1.4). More precisely, let 

n 

\ x - cic = x l + _ 

j = i 

and denote the positive square root of \x - Cic by \x — C|c- In the case that x 
actually lies in {0} x R n = R n , the function C »— > \x — C|c coincides with the 
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analytic extension of the modulus function £ |x - £|, £ E R n \ {x}. Then 

G x {() - — —— xzR n+1 . (3.1) 

a n \X ~ Clc 

Here we take \x — (\^ ^ (— 00 , 0] if n is an even integer and \x — Clc 7^ 0 if n is an 
odd integer. If the dimension n is even, then \x — (lc +1 has a discontinuity as the 
complex number |x — £|^ moves across (— 00 , 0] in C. If n is odd, we only require 
the denominator \x — Clc +1 t0 he nonzero. Then the function (x, 0 •— > G X (Q is 
two-sided monogenic in x and holomorphic in ( on its given domain. 

According to the point of view mentioned in the Introduction, for fixed £ E 
C n , the set of singularities of the Cauchy kernel G X (Q as a function of x E 
R n+1 is called the monogenic spectrum of ( E C n and denoted by 7 (C)- In fact, 
we are just considering ( E C n as an n-tuple of multiplication operators in the 
Clifford algebra C£(C n ) and taking the natural definition of the Cauchy kernel 
x »-> G x ((). We do not need the plane wave decomposition here — the Cauchy 
kernel G X (A) is defined by the L°° -functional calculus for any commuting family 
A = (Ai, . . . , A n ) of normal operators in a Hilbert space. 

To examine the subset 7 (C) of R n more closely, write ( = £ 4- irj for £, rj E R n 
and x = x 0 e 0 + x for x 0 E R and x E R n . Then 

n 

\x - Clc = ^0 + = x 0 + l f -C| 2 - M 2 -2i(x-Cn). (3.2) 

j = l 

Thus, |x — belongs to ( — 00 , 0] if and only if x lies in the intersection hy- 
perplane (x — £, 77 ) = 0 passing through £ and with normal 7 , and the ball 
Xq + |x — £| 2 < \r}\ 2 centred at £ with radius \rj\. If n is even, then 

7(0 = 

{x = x 0 e 0 + x E M n+1 : (x - £, 77 ) = 0, x 2 0 + |x - £| 2 < \p\ 2 }. (3.3) 

and if n is odd, then 

7(0 = 

{x = x 0 e 0 + x E R n+1 : (x - £, 77 ) = 0, Xq + |x - £| 2 = |t 7 | 2 }. (3.4) 

In particular, if Im(£) = 0, then 7 (C) = {(} C R n . 

Now suppose that / E The Cauchy integral formula gives 

/ (C) = [ G x (C)n(x)f{x)dp(x) (3.5) 

JdCt 

for a bounded open neighbourhood of 7(0 with smooth oriented boundary 
outward unit normal n(x) at x E dQ and surface measure p. The integral does not 
depend on £1 because G X (Q is right monogenic in x and / is left monogenic [1, 
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Corollary 9.3]. Differentiating under the integral sign shows that / is holomorphic 
on the open subset 

{7GC n : 7 (C)c5°(r +1 )} of C n , 
which we now describe [4, Proposition 2.1]. 

Proposition 1. Let ( G C n \ {0} and 0 < to < ^ ir . Then 7 (C) C S u ( R n+1 ) if 
and only if 

|C|ct^(-oo, 0] and | Im(C)| < Re(|C|c) tanw. (3.6) 

For each 0 < w < let S^(C n ) denote the set of all ( € C" satisfying 
condition (3.6) and let 5° (C n ) be its interior. The sector S' u; (C n ) arose in [10] as 
the set of ( E C n for which the exponential functions 

= e* (x,c> e _Xo|clc x+(C), x = x 0 e 0 + x, 

have a decay at infinity for all x E R n+1 with (x, m) > 0 and all unit vectors 
m = raoeo + fh E R n+1 satisfying mo > cotcu|m|. 

For each 0 < v < ^7r, let H°°(S°(C n )) denote the set of all holomorphic 
functions 

/ : S°(C n ) -* C£( R n ) 

which are uniformly bounded on every subsector 

S°,(C n ), 0 <v' <v. (3.7) 

Endowed with the topology of uniform convergence on subsectors (3.7) the topo- 
logical vector space C n )) is a Frechet space and a (nonabelian) Frechet 

algebra under pointwise multiplication. The subalgebra of C- valued functions is, 
of course, an abelian Frechet algebra. The proof of the following result is given 
in [3]. 

Theorem 2. The mapping f i— > f given by the Cauchy integral formula (3.5) is 
an isomorphism between the Frechet spaces R n+1 )) and H°°(S°(C n )). 

Given two functions /, g E R n+1 )), the restrictions 

fo : R n \ {0} -> C£(C n ) and g 0 :R n \ {0} — > C£(C n ) 

of / and g to R n \ {0} are real analytic, so their product fogo has a left monogenic 
extension / •£ g to an open neighbourhood of R n \ {0} in R n+1 . The analogous 
product for right monogenic functions is written as / > g. 

Corollary 1. For each f,g E R n+1 )), the function f -e g has a left 

monogenic extension to S°( R n+1 ), denoted by the same symbol, and f •£ g E 
H£°(S°(M n + 1 )). Moreover, R n+1 )) is a Frechet algebra with the prod- 

uct 

(f,9) f-i 9- 

The mapping / »— > / given by the Cauchy integral formula (3.5) is an isomor- 
phism of the Frechet algebras H^°(S°( R n+1 )) and H°°(S°(C n )). 
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If the restrictions /o, go of / and g take values in C rather than C£( R n ), then 
both / and g are two-sided monogenic and 

/ > g = / -eg = f.g = gj- 

Hence, / -e 9 = 9 't f and the subalgebra of R n+1 )) consisting of such 

functions is abelian. 



4.4 Joint spectral theory of operators of type uj 

Suppose that T : V(T) — > H is a single densely defined linear operator acting in 
the Hilbert space H. If 0 < lj < is a. number for which T is of type uj (see 
Definition 1), then the one-dimensional version of formula (2.5) gives a bounded 
linear operator f(A) defined by 

/(T)= iX (A/ “ T) " 1/(A)<iA (4-i) 

for any function / satisfying the bounds (2.3) in S U (C) in the case n = 1. The 
contour C can be taken to be 

{z G C : | Im(z)| = tan^| Re(z)| }, 

with u < 6 < v. The operator T of type uj is said to have a bounded H°°- 
functional calculus if for each uj < v < ^7r, there exists an algebra homomor- 
phism / i — > f(T ) from H°°(S°(C)) to C(H) agreeing with (4.1) and a positive 
number C v such that 

||/(T)||<a||/|| 00 for all /Gff°°(S“(C)). 

The following result is from [11, Theorem 6.2.2] 

Theorem 3. Suppose that T is a one-to-one operator of type uj in H. Then T has 
a bounded H 00 -functional calculus if and only if for every uj < v < ^7r, there 
exists c v > 0 such that T and its adjoint T* satisfy the square function estimates 
r°° 

/ ||V' t (TH| 2 -<c I ,||u|| 2 , ueH , (4.2) 

Jo 1 

f°° dt 

/ \\ipt(T*)u\\ 2 — < cJluf , ueH, (4.3) 

Jo 1 

for some function C)), which satisfies 

nOO It pOO Jj. 

/ i>\t) - = / - = 1, and (4.4) 

Jo 1 Jo 1 

\4>(z) zeS°( C), 

for some s > 0. Here = U( iz ) f or z € 5°(C). 



(4.5) 
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We now use formula (2.5) to generalise this result to n-tuples of commuting 
operators acting in a Hilbert space TL. 

Theorem 4. Let A = . . . , A n ) be an n-tuple of densely defined commuting 

linear operators 

A r .V{Ai)^H 

acting in a Hilbert space H such that ^j = {D(Aj) is dense in H. Suppose that 
0 < to < ~7T and A is uniformly of type to. 

IfT — i(A\ei + • • • + A n e n ) is a one-to-one operator of type uj acting in 
H( n ) and T has an H 00 -functional calculus , then the n-tuple A has a bounded 
H 00 -functional calculus on S°(C n ) for any lo < v < that is, there ex- 
ists a homomorphism b \ — > b(A ), b £ H°° (S°(C n )) , from H°°(S°(C n )) into 
£( n ) (H{ n )) and there exists C v > 0 such that 

\\b(A)W <aiH|oo for all 6 £ ff~(S°(<C n )). 

Moreover, if f is the unique two-sided monogenic function defined on S°( R n+1 ) 
such that b = /, as in Theorem 2, and f satisfies the bound (2.3), then b(A) = 
f(A) is given by formula (2.5). 

Proof. By assumption, the operator T has an H°°- functional calculus, so there 
exists a function f £ H°°(S°(C)) satisfying conditions (4.2) - (4.5). Our aim is 
to define b(A) for b £ H°°(S°( C n )), by the formula 

(b(A)u, v) = J™ ({ b<t> t )(A)MT)u , i) t (T)*vj j (4.6) 

for all u,v £ W( n ). The function f : S°( C n ) — > C is constructed from f by 
setting 

<KC) = ^ 2 R} = V ; 2 (ICIc)x+(C) + V’ 2 (-|C|c)x-(C), 

for all C G S°( C n ). Then for a particular choice of the function f [3], the holo- 
morphic function f t defined for t > 0 by ft(() — f{t() has the property that 

f h f), f£M n \{0}, (4.7) 

has a left (and right) monogenic extension b •£ ft to S'°(R n+1 ) with decay at zero 
and infinity [3]. Hence (bf t )(A) := ( b •£ ft) {A) is defined by formula (2.5) and 
satisfies 



sup \\{b-ef t )(A)\\ <C||6||oo 
t> o 



By normalising f so that 




the desired functional calculus b b(A ), b £ H°°(S°( C n )), is obtained. □ 
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Remark, (i) In [13], it is shown that the n-tuple Du of differentiation operators 
on a Lipschitz surface E C R n+1 mentioned in the Introduction satisfies the 
conditions of Theorem 4. In particular (Du, s) is the differentiation operator on 
the Lipschitz graph determined by the slice of E in the direction 5 £ S' 71-1 , so for 
some 0 < cj < ^7r, the n-tuple Du is uniformly of type uj. 

(ii) A result similar to Theorem 4 is obtained in [1 1, Theorem 6.4.3] using the 
idea that A generates a bounded monogenic semigroup rather than the assumption 
that A is uniformly of type u. I do not know the connection between the two 
concepts. 
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Poincare Series in Clifford Analysis 

Rolf Soren Krausshar 



ABSTRACT In this paper we deal with Clifford-valued generalizations of several 
families of classical complex-analytic Eisenstein series and Poincare series for dis- 
crete subgroups of Vahlen’s group in the framework of Clifford analysis. 

Keywords: Eisenstein series, Poincare series, Clifford analysis, hypercomplex mod- 
ular groups, automorphic forms 



5.1 Introduction 

One very central and important part of complex analysis is the classical theory 
of elliptic modular forms which deals with holomorphic functions on the upper 
half-plane iT + (C) := {z e C; Im(z) > 0} that satisfy for all 2 £ H + (C) the 
transformation law 

f(z) = {f\ k M){z), for all M = (“ ^ £ 5L( 2,Z), (1.1) 

where (f\kM)(z) := (cz + d)~ k f^ One way to construct examples of 
holomorphic functions satisfying (1.1) is to start with a holomorphic generating 
function / : iT + (C) — > C that is bounded on iJ + (C) and that is furthermore 
already totally invariant under the translation group T, the group generated by 

the matrix ^ inducing the translation 2 1— > 2 4- 1 . If we sum the expressions 

(f\kM)(z) over a complete system of representatives of right cosets of 5L (2, Z) 
modulo the translation invariance group T, that is, 

f( z ) E u-2) 

M:T\SL{ 2,Z) 

then (provided k is sufficiently large, so that this series converges) the limit func- 
tion / has the requested properties. These function series are often called Poincare 
series in a wider sense. The simplest non-trivial examples can be obtained by 
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putting / = 1, thus yielding — up to a normalization factor — the classical Eisen- 
stein series 



Gm(z) ~ Y 



(c,d)G Z 2 \{(0,0)} 



(cz -f d) r> 



m > 4 



(1-3) 



which is one of the prototypes for elliptic modular forms. It should be mentioned 
that the Eisenstein series (1.3) were generated originally from the Laurent coeffi- 
cients of Weierstrass’ p-function. 

Already in the first half of the twentieth century an extensive theory of higher- 
dimensional (but complex- valued) extensions of the theory of classical elliptic 
modular forms was developed in the framework of several complex variables the- 
ory, starting with Blumenthal [1] in 1904, continuing with Siegel and Maass in 
the 1930s and early 1940s [24, 31]. 

One of the first authors who considered automorphic forms in a hypercom- 
plex variable was Fueter. In 1927 Fueter constructed in his paper [8] automorphic 
forms and functions related to Picard’s group in a three-dimensional hypercom- 
plex variable. However, the functions treated in his paper are not in kernels of 
Dirac type operators. 

In 1949 Maass proceeded to introduce in [25] a further different kind of au- 
tomorphic forms in a hypercomplex variable. The class of automorphic forms 
introduced by him consists of complex-valued non-analytic eigenfunctions of the 
higher-dimensional hyperbolic Laplace-Beltrami operator. 

Many authors extended these fundamental higher-dimensional theories from 
a large variety of viewpoints. Just to give a few examples, see works of El- 
strodt, Grunewald and Mennicke, [5,6], Freitag and Hermann [7], Krieg [19-21], 
O. Richter and H. Skogman [26, 27] which are some of many authors that pro- 
vided important contributions in this direction. 

If one looks for analogous extensions within the framework of Clifford analysis 
and its regularity concepts, then, as far as we know, one primarily finds only very 
few contributions to higher-dimensional versions of the doubly periodic holomor- 
phic Weierstrass functions as e.g., [4, 9-11, 28] before we started in the period of 
1998 - 2002 to fill this gap step by step, and to develop systematically a theory of 
monogenic and polymonogenic hypercomplex-valued modular forms in the set- 
ting of real and complexified Clifford analysis (see e.g., [13-18]). 

In this paper we explain how a Poincare summation process can be applied 
to some special generating functions from Clifford analysis in order to construct 
non-trivial monogenic and polymonogenic Eisenstein and Poincare series, thus 
providing many non-trivial examples of hypercomplex- valued modular forms re- 
lated to discrete subgroups of Vahlen’s group for the half-space model as well as 
for the unit ball model within the Clifford analysis setting. We further discuss the 
problem of the construction of cusp forms in the setting of two monogenic vector 
variables in which a certain biregular variant of the exponential function plays a 
crucial role. 
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5.2 Notation and preliminaries 



For details about Clifford algebras and basic Clifford analysis see e.g., [3]. Let 
ei , e 2 , . . . , e n be the standard basis of the Euclidean space R 0,n and C^ 0 , n be the 
associated real Clifford algebra in which e^j +e J e 2 = holds, 5 tJ standing 
for the Kronecker function. We write ca •= • • • e ik where A = {zi, . . . , z/J C 

{1, . . . , n}, and in particular e$ = eo = 1. 

By < a, 6 > := Sc(ab) we denote the usual scalar product between two 

Clifford numbers a, b £ C¥ 0 ,n- The Clifford norm of an a = Y cla^a £ C^o,n is 

A 

then ||a|| = (Y |a^| 2 )2 .The reversion is defined by (ab)* = b*ci* and e* = 

A 

for all i = 1, . . . , n and the conjugation by ab — b a and e7 = — e t for all 
i — 1, . . . , n. For a natural p G {1, . . . , n — 1}, 



e p \ (0 -1 

Y'Yj- 

■T r =:Q 

stands for the hypercomplex modular group (compare with [6]) which acts dis- 
continuous^ on the upper half-space 

7/ + (M n ) := {z = x + x n e n : x G M n_1 , x n > 0} (2.2) 



r p := 



e l 

1 



— '■T\ 



( 2 . 1 ) 



via its associated Mobius transformation 



M < z >:= (az + b)(cz + d) 1 = ( zc * -f d*) l (za* + b*) 

where M = ^ ^ G T p . The subgroup generated by the first p translation 

matrices 7\, . . . , T p will be denoted by T p throughout this paper. The Dirac op- 

n 

erator in R 0,n will be denoted by D z := Y Clifford-valued functions, 

i= 1 

defined in open subsets of R n , that are annihilated from the left by the Dirac op- 
erator are called left monogenic and functions annihilated from the right are con- 
sequently called right monogenic. More generally, one speaks of left (resp. right) 
fc-monogenic functions, if D k f — 0 or fD k = 0, where D k means the k - th it- 
erate of the Dirac operator. Note that the Dirac operator factorizes the Laplacian 
viz: D 2 = - A®n . If n is even, then Ker DA n / 2-1 is the class of the Fueter-Sce 
solutions (see e.g., [22, 30]). Clifford- valued functions of two vector variables that 
are left monogenic in the first one and right monogenic in the other one are called 
biregular. For the fundamental solution of D and its partial derivatives (general- 
ized negative powers) we write 

z <9l m l 

qo{z) = q m (z) := g^qo{z), (2.3) 

where m = (mp . . . , m n ) G Nq using the usual multi-index notation as in [16]. 
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5.3 The monogenic Weierstrass p-function and 
associated Eisenstein series 

The Eisenstein series (1.3) were originally introduced by means of the Laurent 
coefficients of the classical Weierstrass p-function associated to a lattice of the 
form Q = Zr -f Z, r £ # + (C) (method of generating functions). Without loss 
of generality one can restrict to such types of lattices, since every general lattice 
can be transformed into such a one by a simple rotation. 

In this chapter we discuss a generalization of (1.3) that is obtained when we 
apply the same method on the higher-dimensional monogenic variants of the 
Weierstrass p-function. Contributions to higher-dimensional generalizations of 
the Weierstrass p-function in the Clifford analysis setting can be found e.g., in [4, 
9-11, 28] and more recently in [12-14, 23]. We recall: 

Definition 1. Suppose that Q n := Zc^i + • • • + Ztj n is a non-degenerate 
n-dimensional lattice in R n . Let r(i) be the multi-index m := (mi , . . . , m n ) with 
mj = Sij , for j = 1 The associated generalized monogenic p-function is 

then defined as 

p(z) ■= <lT(i)(z+u) -9r(i)H • (3.1) 

As in the complex plane, let us restrict to lattices of the form il n _i + Zr 
where r £ H+( R n ) and where O n -i is an (n — l)-dimensional lattice from 
spa%{ei, . . . , e n _i}. Notice that every general lattice can be transformed into 
such a one by a simple rotation. Generalizing the classical idea, the Laurent coef- 
ficients generate thus the following function series (when regarding r as a hyper- 
complex variable of H+( R n )) as a generalization of (1.3): 

Definition 2. (cf [13, 15]) 

Fora multi-index m £ Ng of the form (mi, . . . , m n _i, 0) of odd length , i.e., with 
l m l = E?= ;1 X rrii = 1 (mod 2), the monogenic Eisenstein series associated to 
the generalized monogenic Weierstrass p-function is defined by 

G m (z):= £ q m (az + co), zeH + (R n ), |m|>3, (3.2) 

(a,u;)eZxft n _i\{(0,0)} 

where Q n -i denotes a non-degenerate lattice in span K {ei, . . . , e n _i}. 

One can restrict to indices with m n = 0 as a consequence of the Cauchy- 
Riemann system. The series G m is absolutely convergent and thus monogenic on 
H+( R n ) whenever |m| > 3 (cf. [15]). Additionally we have 

G m (T < z >) = G m (z) 

for all T £ T n -i so that it can be represented by a normally convergent Fourier 
series on iT + ( R n ). Writing z £ H + ( R n ) in the form z — x + x n e n where 
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x G span R {ei, . . . ,e n _i} and x n > 0 then we get, restricting without loss of 
generality to an orthonormal lattice, that its Fourier expansion on H+ (R n ) reads: 



G m {z) = 2CM -1 (m) 

+ A n {2m)^ V a m (s)(ie„ + ^-)e 2Ti<s ’ x> e- 2,r l |s ll a: ". 

^ ^ S 

sGZ n ~ 1 \{0} 11 11 

Here, A n is the surface of the unit ball in R n , Cm" 1 the generalized Riemann 
zeta function from [13] (see [2, 15] for more details) and 

<7 m (s) = 5]r" 

r|s 

where r|s means that there is an a G N such that ar = s. 

The proof can be done in a similar way as we presented in [13, 15] in the 
paravector formalism. The main difference is that the Fourier transform of qo in 
the vector formalism has the different form 

(ie„ + ^.) e 27ri < s . x > e -27r||s||3: n ^ 

INI 

All the other technical steps can be adapted easily. 

One observes a very close similarity between the form of the Fourier expan- 
sion of the classical Eisenstein series (1.3) and that of the higher-dimensional 
variant defined in (3.2). As a consequence of the monogenicity, the monogenic 
plane wave function appears as a generalization of the classical exponential func- 
tion. Infinitely many Fourier coefficients do not vanish, so that the series G m are 
actually non-trivial functions. 

Summarizing: from the structure of the Fourier expansion of the series G m and 
from their connection to the Laurent coefficients of the generalized p-function, 
the series (3.2) provide canonical generalizations of (1.3) in the Clifford analysis 
setting. However, one important aspect is not yet included in this generalization. 
Although the set of singularities of (3.2), i.e., Qei H b Qe n _i is totally invari- 

ant under the full hypercomplex modular group T n _i (like their complex coun- 
terparts), on the half-space we do only observe the invariance under the subgroup 
T p and not a quasi-invariance under the whole group T n _ i . 

Nevertheless, as we shall see in the next section, the functions (3.2) serve as 
useful building blocks for generating non-trivial families of monogenic and poly- 
monogenic vector- valued or Clifford- valued modular forms with respect to larger 
hypercomplex modular groups. This provides in a certain sense still an analogy to 
the complex case. 



5.4 Poincare series in Clifford analysis 

From the conformal invariance formula for the iterated Dirac operators D h (see 
for example [29]) it follows that if / : H + ( R n ) — > Cl§ n is a left monogenic 
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function, then 

D l z {J h (cz + d)*f{M <z>)} = 0 
for every positive integer l with l > h. Here 

J h (cz + d) := ( cz + d)/\\cz + d\\ n ~ hJrl 

whenever h is odd, and 



J h (cz + d) := || cz + d\\ h ~ n 

whenever h is an even positive integer. Starting now with an appropriate mono- 
genic generating function / from Clifford analysis that is totally invariant under 
the action of the translation group T p , and summing then expressions of the form 

Cf\hM)(z) := J h (cz + d)*f(M <z>) 

over a complete set of representatives of right cosets in modulo T p , will lead 
to several classes of non-trivial Clifford- valued modular forms with respect to 
the whole group T p within the function classes Ker D l for all positive integers l 
satisfying l > h whenever h is an even positive integer. We are going to explain 
this now in detail. In the sequel the notation M : T P \T P means that M runs 
through a system of representatives TZ p of the right cosets of T p with respect to 
T p , i.e., U Men p T p M = T p and T P M ^ T P N for M,N e1Z p with M / N. We 
showed in [16]: 

Construction Theorem 1. Let n £ N, h e 2 N, h < n, p E {1, . . . , n — 1} 
and p < n — h — 1. Let f : — > Clo n be a bounded and left monogenic 

function on H+ (R n ) that is totally invariant under the translation group T p . Then 

m= Cf\hM){z), zeH+(R n ) (4.1) 

M:T P \T P 

is a Clifford-valued function which is bounded in any compact subset ofH + (M n ). 
Moreover, it satisfies on the whole upper half-space D l f(z) = 0 for all l > h and 
f{z ) = {f\hM)(z)forallMeT p . 

Short sketch of the proof As / is bounded on H + ( R n ) it suffices to make clear 

that \\cz + d\\ h ~ n is normally convergent on i/ + (R n ) forp < n — h — 1. 

M:T p \r p 

To this end, let e > 0 and consider the vertical strip 

V € (F + (K n )):={2 = x + x„e n Gi/+(M"), ||x|| < - x n > e}. (4.2) 

With a standard compactification argument one concludes that one can find for 
every e > 0 a real p > 0 such that 



\\cz + fif|| > p||ce n +d\\ 
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for all z G V e (H+(R n )) and for all M eT p . The series 

X \\ce n + d\\~ a 
M-.T p \r p 

has the same convergence abscissa as the series 

Y ll ce p-(-i T d\\ 

M:T p \r p 

which is - according to [6] - precisely p + 1. Thus the series (4.1) converges 
normally whenever p < n — h — Iso that / is in Ker D l whenever l > h as 
a consequence of Weierstrass convergence theorem. To prove the automorphic 
behavior it suffices to verify it for the generators of the group. This can be done 
by applying direct rearrangement arguments in the series and using furthermore 
the property Jh{ab) = Jh(b)Jh{a) for elements a, b from the Clifford group. For 
the detailed proof, see [16]. 

By means of the Eisenstein series introduced in the previous section we can 
generate a couple of non-trivial vector- valued examples: 

Theorem. Let n G N, h G 2 N, h < n, p G {1, . . . , n — 1}, p < n — h — 1 and 
suppose that m G Nq _1 is a multi-index where |m| > 3 is an odd number such 
that (^ n-1 (m) / 0. Then , putting G m (z) := G m (z + e n ), the series 

Eh,m(z) = XI (G™\hM)(z) (4.3) 

M:T p \r p 

does not vanish and represents a non-trivial transcendental vector-valued auto- 
morphic form with respect to T p of weight (n — h) on the upper half-space in the 
classes Ker D l for every positive integer l with l > h. 

Notice, (G m \hM)(z) stands for Jh{cz + d)G m (M < z > +e n ). To verify the 
non-triviality, consider lim Eh, m {z) which equals 

x n —>oo 

X li m G m (a[x + x n e n ]dT l + bd" 1 + e n ) 

M:T p \r p ,CM=0 

= 2 P+1 lim G m (z) = 2 P+1 (^/‘ _1 (m). 

X rl — MX) 

It suffices to show that there are indices m for which Cm 1 (m) fz 0. To this 
end notice that Cm" - 1 ( m ) 0 m l °dd) are the Laurent coefficients of the series 

e r (z) = X 9r i z + oj) |r| > 1. 

wGfi n -i 

If Cm 1-1 (m) vanished for all m, then one would obtain that e r (z) ^ q r (z) within 
a whole ball which is a contradiction. 
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As a concrete example consider the class of Clifford- valued functions in Ker D A 4 
defined in R 10 (Fueter-Sce equation). Let us consider the group T2 which acts 
discontinuously on H+( R 10 ). If we take a multi-index m G Nq of odd length 
|m| > 3 for which (^ n-1 (m) ^ 0, then the families of function series 
^ 4 ,m and provide non-trivial examples of vector-valued hypercomplex 
modular forms with respect to T2 in the function class of holomorphic Cliffor- 
dian functions when expressing them in the paravector formalism. 

Notice that if we substituted G m by a general (n - l)-fold periodic function 
/ E Ker DA 4 , then in general f\2aM < z > would not be anymore an element 
of Ker DA 4 , since a has to be chosen out of the set {1, 2, 3} in order to preserve 
the convergence of the associated function series. 

Notice that the construction (4.1) generates only the zero function whenever h 
is odd. See [16] for details. More generally, whenever we have a group that in- 
cludes the negative identity matrix, then all functions satisfying 

f(z) = ( f\ h M)(z ) 

for all M from that group vanish identically, whenever h is odd. 

If we want to construct non-trivial /i-monogenic functions that satisfy the trans- 
formation law f(z) = (f\hM)(z) for odd /i, then we need to restrict on smaller 
groups. Just to give one very simple example of /i-monogenic functions with such 
a transformation behavior involving odd /i, let us take for G the trans version group 

W - : = (* = (« l)" ■ V ’ : =(l !)) (44) 

where p < n — 1. This group does not contain the negative identity matrix. All 
matrices from this group induce transformations of the form 

V(z) = z(vz + l )- 1 

with ai)€ Z p . The expressions V(z) do not have singularities in H + (W l ) and it 
is not difficult to show (compare with [17]) that the series 

Vh(z) := ^2 J h (cz + d) (4.5) 

Mew p 

converges normally on D + (R n ) whenever^ < n—h. Furthermore, a computation 
leads to lim Vh{e n x n ) = 1, providing us with the non-triviality argument of 

X j\ ^ OQ 

this series. This example is indeed one of the simplest ones, since the transversion 
group is conjugated to the translation group T p . 

One possibility to construct also non-trivial modular forms in Ker D l (/ > h) 
where h is odd (in particular for h = 1), with respect to the whole modular group 
r p and even with respect to the full modular group T n _ i , is to consider two weight 
factors (one from the left, another one from the right) and to introduce a second 
auxiliary variable. We exploited this idea recently in [18] in order to construct 
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monogenic and polymonogenic Hilbert modular forms on Cartesian products of 
the upper half-space. 

In what follows we focus on the construction of non-trivial monogenic and 
polymonogenic modular forms on 

tf+(M n ) - H+( R n ) © H + (W l ) 



that transform 



f{z,w) = (f\\h,iM)(z, w) for all M eT p (4.6) 

and whose restriction to the diagonal (z = w) is still a non-constant C°°-automor- 
phic form satisfying 

f(z,z) = ( f\\h,iM)(z,z ) 

for all M e r p . Here 

(f\\h,iM)(z,w) = J h (cz + d)*f(M < z >,M < w >)Ji(wc * +d*). 

The following construction theorem provides many non-trivial examples: 

Construction Theorem 2. Let p < 2n — (h + 1) — 1 and leth + l = 0 (mod 2). 
Suppose f : H £ (R n ) — > C£ o n is a bounded function that satisfies for all ( 2 , w) G 
H 2 (K n ) the relation 

D h z [f(z,w)} = [f(z,w)}D l w = 0 

and furthermore 

f(T < z >,T < w >) = f(z,w) 

for all T €T p . Then 



f(z,w):= ( f\\ h ,iM)(z,w ) (4.7) 

M:T p \r„ 

is left h-monogenic in the variable z and right l-monogenic in w and satisfies for 
all (z,w) £ R n ), 

f(z,w) = (f\\ h jM)(z,w), for all M G T p . (4.8) 

Sketch of the proof (for details see [18]). Let JC be a compact subset of 
H 2 (M n ). Choose an £ > 0 such that 

V £ (H+(W 1 )) X V £ (H' ¥ (W 1 )) 

covers /C where V £ (H+(R n )) is the vertical strip domain defined in (4.2). With 
the argumentation of the convergence proof of the first construction theorem, we 
know that there is a real p > 0 such that for all M G T p : 

\\Jk(cz + d)\\ < p h - n \\ce n + d\\ h ~ n (4.9) 
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for all z that belong to V £ (H+( R n )). Since / is bounded, we obtain that for every 
(z,w) E /C : 

| Y II J h (cz + d)*f(M < z >,M < w >)J t (wc* + d*)\ 

M:T p \r p 



< Lp h+l ~ 2n y 



M:T p \r p 



1 1 
||ce n + d\\ n ~ h ||e„c* 4-d*|| n_J 



< 



M:T p \r p 



1 



ce n + < 



|2n-h-Z ’ 



where L and L are non-negative real constants. The series in the previous line 
is absolutely convergent if and only if p < 2n — h — l — 1 (see [6, 16]). From 
Weierstrass’ convergence theorem and the quasi-conformal invariance formulae 
for the iterated Dirac operators from [29], we infer that the limit function / is left 
/i-monogenic in z and right /-monogenic in w. To show the automorphy relation 
it is again sufficient to verify the relation for the generators of the group. With 
respect to the inversion we observe that 

W> -1 ) = Y [jh{-c z ~ 1 + d)* 

M:T p \r p 

x /((— az~ l + b)(—cz ~ 1 + d) -1 ), (— w~ l c* + d*) -1 (—w~ l a* + 6*)) 
x Ji(— w~ l c* -h d*) 

= J h {c - dz)* x 

M:T p \r p 

f ((a — bu)(c — du )~ l , (c* — wd*)~ 1 (a* — J/(c* — wd*) Ji(—w ~ 1 )~ 1 

By a similar rearrangement procedure we can show that the automorphy relation 
holds for simultaneous translations (z, w) (z + m, w + m) where m = e* for 
z = l,...,p. 

We get even convergence for p — n — 1 whenever h + l < n — 1 . Even for 
the monogenic case, i.e., h = l — 1, we get non-trivial Clifford- valued modular 
forms for the full modular group provided we are in a space of dimension n> 3. 
The simplest examples are the biregular Eisenstein series 

£{z,w)= Y ( 1 lll,l M )( Z > U ’) 

M:T p \r p 

= £ J fc (cz + d)*Jj(«;c* +d*), (z,u>) G H+(R n ). (4.10) 

M:T p \r p 
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To verify that these series are non-trivial transcendental functions, consider 



lim S(g n x n ^Q n y n ) 

x n ,y n -+oc 



E 



M :T P \r p , cm 7^0 v 



lim Ji(ce n x n + d)*Ji(e n y n c* + d*) 

iVn ^OO 



=0 



+ 53 Ji(d)*J x (d*) 

M:Tp\r p ,CM=0 



= 2 e A eZ = 2r + \ 

AC{l,...,p} 

where we put z = Sj=i x j e j an d 1(7 — ]Cj=i Vj e j • Also their restriction to the 
diagonal £(z, z) is thus a non-constant C°°-automorphic form with respect to T p . 
Once again, the series G m can be used to generate classes of non-trivial examples 
for hypercomplex modular forms transforming as (4.6) — in this setting, even for 
obtaining families of non-trivial monogenic modular forms for the full modular 
group r n _i — namely by applying (4.7) on 



f{ z i w ) — G m (z + e n )G m (w + ©n) (4.11) 

which is bounded on H ^ (K n ). To show the non-triviality we consider again 



lim Y Ji(ce n x n + d)*G m (M < e n x n > +e n ) 

x n ,y n -+o o z ' 

M:T p \F p 

X (Af < G nVn ^ "b^n) ~b d ) 

= (-8) £ ^||C(m)|| 2 eI*. 

AC{l,...,p} 

As mentioned above, there must be at least some indices m G Ng ~ 1 of odd length, 
for which C( m ) and hence ||C(m)|| 2 is different from zero. For these indices we 
have then that the limit considered is different from zero, proving the non-triviality 
of the associated series 

53 (f\\h,iM)(z,w) 

M:T p \r p 

and in particular that it is non-constant on the diagonal (z, z). 

In this paper we proceed now to treat also biregular and polybiregular mono- 
genic modular forms / that satisfy additionally 

lim f(z,w) = lim f(z,w) — 0. (4.12) 

x n — >• oo y n — ► oo 

To this end it is suggestive to use the following variant of a biregular exponential 
function on (M n ) : 
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Definition 3. For p, q € R" J \{0} and (z, w ) G H£ (l n ) define 

E pq (w,z) := (ten + ^)e- x " | l p|l e i<x ' p> e Ky ’ q> e- y ” l|q|l (ie n + ^L). 

This function is left monogenic in 2 and right monogenic in w on the whole 
domain H 2 (K n ). Notice that left and right monogenicity is provided by the con- 
stant factors standing left and right from the exponential expressions. It is crucial 
that 

lim E pq (z,w) = lim E pq (z,w) = 0 

x n ^00 y n >00 

for any parameter p and q from R n_1 \{0}. In order to construct hypercomplex 
modular cusp forms that are left /i-monogenic in 2 and right /- monogenic in w, 
the following construction is thus plausible: 

rVz(z,iu;p,q)= E (£p,ql \h,iM)(z,w). (4.13) 

M:T p \r p 

The function E p , q is bounded on H 2 (M n ); thus the convergence abscissa of the 
series is p < 2n — (h + l) — 1. However, until now we have not managed yet to 
find an argument to show that these series do not vanish. Note that by construction 

lim 7\/(2,w;p,q) = lim Vh,i( z , w \ P,q)=0. 

x n — ► 00 y n — >■ 00 

At the moment it is not clear how to extend to Clifford analysis the non-triviality 
arguments that are applied for cusp forms in the classical theory, as for instance 
in the one complex variable case. 

Further possible candidates for modular forms with (4.12) that are left h- mono- 
genic in 2 and right /-monogenic in w might be obtained by inserting for / in the 
expression (4.7) the function 

f{z, w) := e m (z + e n )e r (w + e n ) 

where |rj and |m| are assumed to be greater than or equal to 1, since 

lim e m (z + e n )=0 

X n — too 



whenever |m| > 1. 

Finally, we also want to give some examples of non-trivial hypercomplex 
monogenic and polymonogenic modular forms on the Cartesian product of two 
unit balls in this paper. 

Let 



9 = 



/ 1 

y/2 

C72 e 



V2 en 

1 

x/2 



Then the associated Cayley transformation 



(z,w) -> (g < z >,g <w >) 
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maps H .2 (M n ) conformally onto 

fi 2 (0,l) := 5(0,1) x 5(0,1) = {(z, w) G M n x R n ;||z|| < 1 and ||w|| < 1}. 

In order to construct monogenic and polymonogenic automorphic forms on 
B 2 ( 0, 1) it is suggestive to make a construction of the form 

52 (f\\h,iM)(z,‘w) (4.14) 

M:gG'g-'\gGg-i 

where G < G < T n _i and where / is a left /i-monogenic function in z and 
/-monogenic in w that is bounded on B 2 (0, 1) and totally invariant under gG g~ l . 
Here, we observe that the transfer to the unit ball is not as direct, as one might 
logically think at the very first: Note that if / is monogenic in z, then in general 
f(M < z >) will only remain monogenic, when M is a translation matrix. 

Putting G — T n -i (which is natural), then not every monogenic / has the 
property that f((gGg~ l ) < z >) remains monogenic. The same holds also for 
/- monogenic functions. The simplest non-trivial (poly-) biregular examples that 
we can construct to the whole modular group G := gT p g~ l on B 2 (0, 1) and for 
which we can avoid this problem, are the following (poly-) biregular Eisenstein 
series on the Cartesian product of the two unit balls. 

We introduce: 

Definition 4. Let p <2n — (h + l) — 1 where h + l = 0 (mod 2) and let 
Gp 9^p9 \ G' p :=gT p g 1 

where g stands for the Cayley transformation. Then, the (poly-) biregular Eisen- 
stein series on the Cartesian product of the two unit balls B 2 ( 0, 1) are defined 
by 

H(z,w):= 52 (111 h,iM){z,w). (4.15) 

M:G’ p \G„ 

The function /(z, w) = 1 is of course left /i-monogenic and right /-monogenic 
and totally invariant under G p . Thus, the limit function is left /i-monogenic and 
right /-monogenic. The series has the same convergence abscissa as its analogue 
on the half-space, and its automorphy behavior 

H(z,w) — (Tt\\h,iM)(z,w) forall M e G p 

can again be shown similarly as we did for the previous examples. To show 
that the series does not vanish identically on the diagonal whenever h -f / = 0 
(mod 2), evaluate the function at the point (0,0), which of course is a (poly-) 
regular point as it is inside the unit balls, and all the singularities are concentrated 
on the boundary of B 2 ( 0, 1). The series is therefore normally convergent in any 
compact set around (0, 0) inside B 2 ( 0, 1). Thus, 

H(0,0)= £ J h {d)*Md*)= 52 Nil 2 ^ 0 ( 4 -!6) 

M:G' P \G P M:G' P \G P 

and the non-triviality of the series is proven. 
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Harmonic Analysis for General 
First Order Differential 
Operators in Lipschitz Domains 

Emilio Marmolejo-Olea and Marius Mitrea 



ABSTRACT We develop a function theory associated with non-elliptic, variable co- 
efficient operators of Dirac type on Lipschitz domains. Boundary behavior, global 
regularity, integral representation formulas, are studied by means of tools originat- 
ing in PDE and harmonic analysis. 

Keywords: Cauchy type integral operators, Lipschitz domains, Hardy spaces 



6.1 Introduction 

In this paper we discuss some basic aspects of the function theory associated 
with a general first-order differential operator D , such as: Pompeiu type inte- 
gral representation formulas, Cauchy type integral operators, and Hardy spaces 
(with emphasis on boundary behavior and global regularity), in a subdomain Cl 
of a manifold M. For maximum applicability, it is important to allow minimal 
smoothness assumptions on dCt , variable coefficients and vector character for D, 
as well as operators with non-invertible (principal) symbols. 

While our primary source of inspiration in pursuing the program sketched 
above has been the function theory of the classical Dirac operator associated with 
C£(R rn ) — the standard Clifford algebraic structure in R 771 — our goal here is to 
explore the extent to which such a theory can be essentially worked out by real 
variable methods , in the absence of the rich algebraic-geometric structural prop- 
erties of C£(R m ). 

Thus, loosely speaking, the aim is to do Clifford analysis without (extensive 
use of) Clifford algebras, whenever possible. On the other hand, results which are 
fundamentally dependent on the underlying algebraic structure are recognized as 
such. 

This point of view is rather useful in many applications, particularly to PDE 
problems. Specific problems arising in mathematical physics have been treated in 
[ 7 , 14 , 15 , 20 ]. 
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The layout of the paper is as follows. In §2 we collect some basic definitions, 
and discuss an important integration by parts identity. In §3 we introduce the 
class of Dirac type operators — via a nonstandard definition — and present some 
relevant examples. Section 4 is devoted to discussing the role of unique continu- 
ation. Cauchy type operators are introduced in §5; here we also prove a general 
Pompeiu type integral identity. Some special cases have been previously treated 
in [3] and [24]. 

The interior regularity and boundary behavior of the aforementioned Cauchy 
integral operators make the object of §6. Here we develop the main tools needed 
in §7, where Hardy spaces in Lipschitz domains are studied. Our results in this 
part of the paper extend those of [6, 19]. In the next three sections, §8-10, we spe- 
cialize some of our previous results to specific cases of particular interest: in §8 
we look at Clifford separately monogenic (or multi-monogenic) functions, while 
in §9 we consider a similar theory in the context of octonions. In §10 we work 
out the details of Pompeiu’s formula for operators such as D = d + d* , and its 
complex counterpart, D = d 4-9* (rather in matrix form). Among other things, 
these are shown to contain the classical Bochner-Martinelli-Koppelman formu- 
las as particular cases. In §11, we list several open problems closely related to 
the theory outlined in this paper. Finally, there is also an appendix, containing a 
self-contained proof of a Hartogs type theorem in the Clifford algebra context. 

Acknowledgments. The work of M. M. has been supported in part by an NSF 
grant and a University of Missouri Research Board grant. The authors would like 
to thank Rafal Ablamowicz and John Ryan for their efforts to make the Cookeville 
meeting a success. The research described here was initiated when E. M.-O. vis- 
ited M. M. at the University of Missouri-Columbia in 1999-2000. Both authors 
would also like to thank Dorina Mitrea for her role in bringing them together. 



6.2 Background material 

Let M be a smooth, compact, boundaryless, orientable, Riemannian manifold. 
Depending on the context, we may allow M to be either a real or a complex man- 
ifold and occasionally we may take M to be the entire (real or complex) Euclidean 
space. We denote by dV the volume element and by d the usual exterior differen- 
tial operator on M. Other standard notation will be adopted without any special 
mention. Next, consider £ , T — > M two smooth, Hermitian vector bundles and 
let the linear map 

D : C°°(M,£) — > C°°(M, J 7 ), (2.1) 

be a differential operator of order k > 1 . According to a well-known character- 
ization of J. Petree, this is equivalent to the property that D is local , i.e., D does 
not ‘increase’ the support. In the sequel, the superscripts ‘star’ and ‘t’ are used 
to denote (complex) adjunction and (real) transposition, respectively. Also, ‘bar ’ 
stands for complex conjugation. Define Du := Du, so that D* = (D) 1 = ( D *). 




Harmonic Analysis for Dirac Operators 93 



Finally, ( •, • ) will denote the bilinear (rather than sesqui-linear) pointwise pairing 
in the fibers of a Hermitian vector bundle (each time, its specific meaning should 
be clear from the context). In particular, \u \ 2 = (u, u). 

The principal symbol of D in (2.1) is the map 

a{D\ •) : T*M \ 0 — -> Hom(£, Pj (2.2) 

defined as 

(a(D-Oe) a = (2.3) 

| 7 |=/c 

if ( Du) a = Z\ y \=k a^d y up + lower order terms, in local coordinates. 

Examples of first-order differential operators and their symbols: 

• If ‘wedge’ A stands for the usual exterior product of differential forms on 
M, then a(d ; £) = A • for any £ e T*M = K l TM. 

• If the ‘upside-down wedge’ V denotes the interior product of forms (i.e., 

the adjoint of wedge), then cr(d*; £) = V£ £ A l TM. 

• Denote by C£(M) the exterior power bundle on M, equipped with the Clif- 
ford algebra structure induced by the Riemannian metric on M. Then, if 
D = d + d*, it follows that cr(D; £) = A • - V • =: z£-, the Clifford 
algebra multiplication with in C£(M). 

• More generally, let £ — > M be a Hermitian vector bundle with the prop- 

erty that each fiber £ x is a Hermitian C£(M) x -modu\e. Assume that £ is 
equipped with a connection V : C°°(M,£) — > C°°(M,T*M © £), and 
denote by m : ©> £) — > C°°(M,£), m(£ 0 e) = £ • e, the 

Clifford multiplication map. Finally, introducing D := m o V, we have 
a(D]£) = z£-, the Clifford algebra multiplication with in for each 
£ £ T*M. See [25] for more details. 

• In the case when M = R m , the flat, Euclidean space, and with ei, e 2 , . . . , 
e m denoting the standard orthonormal basis in R m , the above construc- 
tion yields D := J^e a <9 a , the canonical Dirac operator in R ni (see [2]). 
Once again, a(D ;£) = z£-, the Clifford algebra multiplication with = 

^ • 

• Assume now that M is a complex manifold and consider the decomposition 
A l TM = A 1 ’°TM © A °’ 1 TM. For a generic 1-form $ we write f - $ 1 «° + 
£°T. Accordingly, the exterior differential operator d splits as d = d + d. 
Then a(d; £) = i £ 1,0 A • and cr(<9; £) = i £ 0,1 A •. 

Lemma 1. For any differential operator D : £ — > T of order /c, we have 
a(D t ; £) = ( — l) fc a(D; £)* and ct(jD; £) = ( — l) k o(D; £). In particular, 

a(D*-,Z)=a(D;Zr. 



Finally, a(DiD 2 ;£) = a(D 1 ] ^)a(D 2 ;^) . 




94 Emilio Marmolejo-Olea and Marius Mitrea 



An integration by parts formula for first-order differential operators, very useful in 
the sequel, is recorded below. Before stating it, recall that a Lipschitz subdomain 
of M is an open set with the property that, in appropriate local coordinates, 
dD is given by graphs of (Euclidean) Lipschitz functions. We denote by ds the 
surface element on dQ (inherited from the metric on M), and by v E T*M the 
outward unit conormal to Q defined <is-a.e. on dQ. 

Proposition 1. LetQ be a Lipschitz subdomain of M. Then, for any u G C l (Q,£) 
and v G C l (fl,J : ),we have 

[ (Du, v)dV = l ( u,D t v)dV — ( (io(D\v)u,v) ds. (2.4) 
Jn Jn Jan 

See also [25]. In closing, we would like to point out that (2.4) is the departure 
point in setting up areolar type derivatives, in the spirit of D. Pompeiu’s original 
derivee areolaire , as well as Morera type theorems for the operator D. 



6.3 Dirac type operators 

The Laplace operator associated with a first-order differential operator as in (2.1), 
in short the D-Laplacian , is the self-adjoint, nonpositive definite, second order 
differential operator -D*D : £ — > £. Call D of Dirac type if its associated 
Laplacian has the property that 

o(D*D ;£) is scalar, real and invertible, V£ G T*M \ 0. (3.1) 

The reader should be aware that this definition is non-standard (the ones com- 
monly used in the literature are more restrictive). In particular, we do not require 
D to be elliptic (i.e., have an invertible symbol); this latter case has been consid- 
ered in [20]. 

It is clear from definitions that our class of Dirac type operators is stable un- 
der conjugation and transposition and, hence, under adjunction also. A typical 
Dirac type operator is contained in the fourth example in §2; more examples are 
provided below. 

6.3.1 Examples of Dirac type operators 

• Consider the operator D d + d*. Upon recalling that d 2 = 0 and that 
A = - dd * — d*d is the Hodge-Laplacian , it follows that D = D* and 
—D 2 = A. Hence, D is of Dirac type. A variant is to take D := (d,d*) 
and the same conclusion holds. 

• A discussion similar in spirit applies to the case when D := d + d*. This 
time, d 2 = ( d *) 2 = 0 and the complex Laplacian 



□ : = - dd * - d*d 



(3.2) 
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has a real, scalar, invertible principal symbol; see Theorem 3.16, p. 154 in 
Kodaira’s book [12]. Since D* = D and - D 2 = □, we may conclude that 
D is of Dirac type. In particular, so is D = 3 + d*. Another related Dirac 
operator is D := (d,d*). 

Assume that 

D ■■=(!• *)• < 33 > 
Then —D*D = —DD* — A/ 2 X 2 , so that D is of Dirac type. 

Consider next the case when M is a complex manifold and set 



D := 



d 

d* 



d* 

d 



(3.4) 



It follows that —D*D = —DD* = 0 / 2 x 2 , so that D is of Dirac type. In 
particular, so are 




d* 

d 



D* = ( d ; £ ) and D* 



a* d 

8 a* 



When M = R m and D := then D* = D and - D 2 = A, the 

flat-space Laplacian. That D is self-adjoint is a consequence of the identity 
(a-f,g) = (/,a c -^f), where the superscript c indicates standard conjugation 
in the Clifford algebra C£(R rn ) (see [2]). Thus, D is of Dirac type. 

In the case when M = M m+ 1 and D is the Cauchy-Riemann operator 

m 

D'=d 0 + ^2e a d a , ( 3 . 5 ) 

a = l 

we have D* = —do -f Nonetheless, ultimately, -D*D = A, the 

Laplacian in M — R m+1 , so D in (3.5) is of Dirac type. 

Let M — R m and take D := i e ad a + fce m+ i, where k e R. Then 
D* — e ocd a ~ kem^i so —D*D = A + /c 2 , the Helmholtz operator 

in R m . Thus, D is a Dirac type. A similar conclusion holds for D := ik + 
YJ2=i e 0ida With k <E R. 

Analogous constructions (as in the last two examples above) can be carried 
out employing the embedding R 4 H, the skew field of quaternions. That 
is, if z, j, k are the standard anticommuting imaginary units in H, consider 
D := d X0 +id Xl +jd X2 + kd X3 , sothatD* = -d Xo +id Xl +jd X2 + kd X3 . 
It follows that -D*D = -DD* = A, so D is of Dirac type. 

Analogous considerations apply to the case of the Dirac operator associated 
with octonions\ see §9 for more details. 
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• A simple, yet important, example of a Dirac type operator whose symbol 
is not invertible is D := grad, on scalar functions. Then D* = -div, so 
—D*D = A, the Laplace-Beltrami operator on M. 

• Suppose that several manifolds Mj and vector bundles Tj — > Mj, j = 

1 have been given along with a vector space A (identified with a 

trivial vector bundle over each Mj). It is assumed that on each manifold Mj 
a Dirac type operator Dj (mapping .4- valued functions into sections in Tj) 
has been fixed. Consider then the product manifold M := Mi x • • • x M n 
together with bundle F := @Tj and, finally, the operator D := {Dj)j. That 
is, for each u : M — > A, we set Du = ( Dju)j , where DjU means that Dj 
acts with respect to the variables on Mj . 

It can then be checked that D*(uj)j = ^ D*Uj so that, ultimately, 

-rio> = £(-£*£,). 

It follows that D is of Dirac type. It should be pointed out that, for n > 1, 
the operator D is not elliptic (in the sense that its symbol is not invertible). 

This last example contains, in particular, the Dirac operator relevant in the case 
of separately monogenic (or multi-monogenic) functions. More specifically, con- 
sider the situation when M\ — M 2 = • • • = M n = R m , A = C£(R m ), and 
Dj := the standard Dirac operator with respect to the j-th copy of R m . Set 

© := (Dj)j. (3.6) 

Then u : R m x • • • x R m — » C^(R m ) is separately monogenic, i.e., DjU — 0 for 
each j, if and only if Du = 0. When each Dj is as in (3.5) and m — 1, then, of 
course, we are talking about holomorphic functions of several complex variables. 



6.4 Unique continuation property 

Suppose we are interested in inverting a first-order differential operator D as 
in (2.1) between appropriate Sobolev spaces. To this end, it is convenient to 
consider first the case of its associated Laplacian, —D*D. Formally, D~ l = 
— (—D*D)~ 1 D*, so it suffices to invert the operator —D*D. Since 

a(-fl*D;0 = -a(D;0V(D;a 

it follows that -D*D is strongly elliptic , i.e., the (symmetric part of its) principal 
symbol is negative definite, if and only if 



o(D ; £) is injective £ T*M \ 0. 



(4.1) 
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Assuming that (4.1) holds, it follows from the classical theory that —D*D has in- 
dex zero and a finite-dimensional kernel. In particular, the existence of an inverse 
(-D*D)~ l hinges on whether 

Ker (—D*D) — Ker D — 0. (4.2) 

Many times in practice, however, the above condition fails to hold (although there 
are remarkable exceptions). For example, if D = d + d*, then -D*D — A and 
the triviality of Ker A ultimately depends on certain topological characteristics of 
the manifold M. 

An alternative is to require the weaker condition that 

Du = 0 on M, u = 0 on an open subset of M => u = 0 in M. (4.3) 

In the sequel, the operator D is said to have the unique continuation property 
(abbreviated D G UCP henceforth) if (4.3) holds. The deep result of N. Aron- 
szjan [1] asserts that all elliptic, second-order differential operators with a real, 
scalar principal part have the unique continuation property. This immediately im- 
plies that all Dirac type operators have the unique continuation property. The case 
of (rough) lower-order perturbations has been considered more recently in [8]. 

Consider now a fixed, scalar, real, nonnegative, smooth potential V on M, and 
introduce the formally self-adjoint, nonpositive, second-order differential opera- 
tor 

L:=-D*D-V, L:£ — >£. (4.4) 

Thus, L is a zero-order perturbation of the D-Laplacian. To discuss conditions 
under which L~ l exists, recall that for a Lipschitz domain (7, the (zero-trace) 
Sobolev space £) is the closure of Cff(ft, £) in the norm 

IMIl2 + ||V-|| L 2. 

Proposition 2. Let D : £ — > T be a first-order differential operator for which 
( 4 . 1 ) holds and such that D G UCP. Then for each Lipschitz domain ft C M 
(possibly the whole manifold M ), 

u G H^ 2 (ft, £ ) & Lu = 0 in ft u — 0 in £l. (4.5) 

In particular, L~ l exists ifV > 0 on some open subset of M. 

Proof To justify (4.5), if Lu = 0 in £1 for some u G Hq ,2 (Q, £), then 

0 = f (Lu, u) dV — - [ (\Du\ 2 + V\u\ 2 )dV. (4.6) 

Jn Jn 

Hence, Du = 0 in £l and u = 0 in supp V H U. In particular, with tilde denoting 
the extension by zero outside we see that u G H l ' 2 (M,£) and Du — 0 
on M. Also, u vanishes in some open subset of M. Indeed, this is clear when 
ft is a proper subset of M, whereas when ft is the whole manifold, the desired 
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conclusion follows from the fact that u = 0 in supp V. At any rate, since D 6 
UCP, these considerations imply u = 0 in M and, ultimately, u = 0 in Cl, proving 
(4.5). 

That Ker L = {0} on M also proves that L -1 exists globally. □ 



6.5 Cauchy type operators 

Assuming that the conditions (4.1) and (4.3) hold, let the distribution 

EeV'(MxM,£ ®£) 
be the Schwartz kernel of L~ l . Then 

L x E(x,y) = LyE(x,y) — S x (y ), the Dirac distribution on M. (5.1) 

The self-adjointness of L then translates into E(x , y )* = E(x, y). Set 

T(x, y) := io (D*; v y )D y E(x , y), T G V'(M x M, £ 0 5). (5.2) 

Loosely speaking, E(x,y) and T(x,y) play the roles of the Newtonian kernel 
and the Cauchy kernel, respectively. It is then natural to introduce the Cauchy 
type integral operator 

Cf(x):=[ ( T(x,y),f(y))ds y , x <£ <90, (5.3) 

JdQ 

where / : dQ, —> £ is an arbitrary section. Another way of understanding the 
operator C, is to observe that the first-order differential operator 

d y := ia(D *; v)D (5.4) 

is the conormal derivative corresponding to the factorization: top part of L — 
—D*D. Consequently, dy := v)D is the complex conjugate of d v . Thus, 

on account of (5.2 ), C can be thought of as the natural double layer potential 
operator associated with the D-Laplacian. Let us also introduce the single layer 
potential operator 

Sf(x) := [ { E(x,y),f(y))ds y , xeM\dCl. (5.5) 

JdQ 

Our main result in this section is a generalization of the classical Pompeiu formula 
(from complex analysis, to the current setting). 

Theorem 1. Suppose that (4.1) and (4.3) hold for D. Then , for each xGd and 
u G C 2 (Q , £) we have 

u{x) =C(u\ d n)(x) - S(Vu\ dn )(x) - [ (D y E(x, y), (Du)(y)) dV y 

Jn 

= C(u\an) - S(dyu\ dn ) - [ (E(x,y),(Lu)(y)) dV y . 

Jn 



(5.6) 
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In particular, 



V = 0 on 0 and Du — 0 in Cl => u = C(u\dn) onCl, (5.7) 
Lu — 0 in 0 => u — C(u\dn) — S(d ]y u\dn) on Cl. (5.8) 



Proof. Since L t = —D l D — V, both identities in (5.6) follow starting with 

u(x)= [ ([-D t y 5 y -V(y)]E(x ) y),u(y))dV y , 

Jn 

then invoking (2.4) twice. □ 



Note that a more ‘standard’ form of Pompeiu’s formula is obtained from the 
first line of (5.6) by requiring that V = 0 on Cl. Also, (5.7) plays the role of 
Cauchy's reproducing formula in this setting, whereas (5.8) is Green’s represen- 
tation formula for null-solutions of L (in terms of double and single layer poten- 
tials). It should be mentioned that the conditions on the boundary behavior of u 
can be considerably relaxed. E.g., a proper control of the nontangential maximal 
function of u will do; see §6. 

An intriguing aspect is that while the reproducing formula (5.7) holds for any 
null solution of D (which is reasonably well behaved near the boundary), in gen- 
eral it is not true that DC = 0. It is therefore natural to try to find necessary 
conditions for this latter condition to hold. 

Proposition 3. Under the additional assumption that 

-DL~ l D* = I in T, (5.9) 

the Cauchy operator (5.3) has the property that DC — 0 in Cl. 



Proof. At the level of Schwartz kernels, (5.9) entails — D x D y E(x,y ) = 5 x (y ), 
the Dirac distribution with mass at x. Then, the desired property is an immediate 
consequence of this identity and the definition (5.2). □ 



Remark 1. It is worth pointing out that (5.9) can only hold when D is elliptic 
and has a trivial kernel (globally, in L 2 ). When Ker D = 0, the potential V can 
be chosen to be zero on M, without affecting the existence of L~ l . If in addition, 
C — T and D*D — Z9D*, then (5.9) automatically holds. This holds for most 
of the operators in §5 when M is, e.g., the entire space (the exceptions are the 
operators in the last two examples). 
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6.6 Global regularity and boundary behavior 

For a fixed, Lipschitz domain (IcM, let • |ar2 denote the nontangential boundary 
trace operator. That is, 



uLdz) := lim u(y), x G 90, (6.1) 

u ye i{x) 

where *y(x) C Ll is a suitable nontangential approach region; see [17]. Also, let 
M stand for the nontangential maximal operator, i.e., 

Mu(x) := sup{|u(y)|; y G 7 ( 2 ;)}, x G dLl. (6.2) 

We shall also work with the boundary version of (5.3) and (5.5) given, respec- 
tively, by 



Cf(x) := p.v. f {T{x,y),f(y))ds y , x G 90, (6.3) 

Jaa 

Sf{x):= [ ( E(x,y),f(y))ds y , x G 90. (6.4) 

Jaa 

In (6.3), ‘p.v.’ indicates that the integral is taken in the Cauchy principal value 
sense, i.e., by removing small geodesic balls and passing to the limit. Finally, we 
set $2+ := Lt and := M \ Ll. 

The next two theorems collect some of the basic properties of the operators 
(5.3)-(5.5) and their boundary versions (6.3)-(6.4). To state our first result, denote 
by L^(dLl,£) the Sobolev space of sections / G L p (dLl,£) with |Vtan/| £ 
LP{d 0) and set L p _ x (dVL,£) := (Lf (90,£))* if J + J = 1 . 

Theorem 2. Let LI be a Lipschitz subdomain of M and assume that V is an ar- 
bitrary, smooth, nonnegative, non- identically zero, scalar potential. Also, assume 
that the first-order differential operator D as in (2.1) satisfies (4.1) and (4.3). 
Then, for each 1 < p < oo, the following are true: 

(1) \\fif{V s Cf)\\ LP (dn) < K\\f\\ L r(dn,£)if or s £ { 0 , 1 }; 

(2) C is bounded on L p (dLl,£) and on L\{dLl,£)\ 

(3) Cf\dn± = (± 5 / + C)f for any f G L p (90,£); 

(4) 9„S/|en ± = (^I + C*)fandSf\an ± = Sfforanyfe L p (9 0,5); 

(5) C o S = S oC* on 90, and L(Cf) = 0 , L(Sf) = 0 in 0± ; 

(6) d v C : L\(dLl,£) — > L p (dfl,£) is bounded, and d v Cf\dn + = dyCf |an_ 
for any f G L\(dLl,£)\ 

(7) The adjoint of the operator above is d u C : L q (dLl,£) — > £), 

where j | = 1; 
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(8) S : L p (9f2,£) — > L\(dVt,£) is bounded and its adjoint is the operator 

S : -> L<>(m,£), where ± + ± = 1; 

(9J ||A/ r (V s <S/)|| iP (an) < K ll/llLP_ 1 (an,£)/ or s e {Oil}; 

(10) II |V 1+s C/|dist(-,afi) 2 || L 2 (n) < K||/|| L 2 (an)£) /ors G {0,1}. 

Next we consider the action of C and C on other types of Sobolev and Besov 
spaces, denoted by L p and B p ' q , respectively. These are considered either on 0 
or 90; see [9] for detailed definitions. 

Theorem 3. With the same hypotheses as above, the following operators are 
bounded on the indicated spaces: 

(1) C on L p (d£l, £) for each 1 < p < oo and 0 < s < 1; 

(2) C on £)for each 1 < p, q < oo and 0 < 5 < 1; 

(3) S : L p _ s (d£l, £) — > L p _ s (d£l, £)for 1 < p < oo and 0 < s < 1; 

(4) S : jB^(90, £) — > B p 'f s (d£l 1 £)for 1 < p, q < oo a/zd 0 < s < 1; 

(5) 9*,C : L p (dQ, £) — > L^_ x (90, £)/br 1 < p < oo and 0 < 5 < 1; 

(6) : B p ' q (d£l,£) -> 5^(90, £)/or 1 <p,g < oo and 0 < 8 < 1; 

(7) 5 : £) — > S^ s+1//p (0, £ )/or 1 < p, q < oo, and 0 < s < 1; 

(S) C : L p (dQ.,£) -» L^ +1/p (0,£) n B p f P l/p (Q,£) for 1 < p < oo and 
0 < s < 1; 

(9) C : B p ' q (dQ, £) —> £^^(0, £)/or 1 < p, q < oo, ami 0 < 8 < 1; 

(10) C : L p (dQ,£) — > B^ # ^(0,£) f/1 < p < oo, with p# := max{p, 2}, 
and 8 G {0, 1}; 

(ID S : VL a (jKl,£) -> L^_ s+1/p (0,£) n B™ a + l/p (tl,£) for 1 < p < oo 
and 0 < 8 < 1; 

(72) 5 : L^Vd^S) -> B p s f* /p (V,£)for 1 < p < oo and s e {0,1}. 

The proof of Theorems 2-3 can be carried out using the techniques developed 
in [17] and [20]. These build on the work of many authors, including [4, 9, 23], 
(see also the references in [11]). For the case of constant coefficient operators see 
[19]. 

As pointed out in the statement of Theorem 1, null-solutions of D are ‘repro- 
duced’ by C. Next, we consider the converse issue (see also [13]), i.e., whether 
null-solutions of L which are reproduced by the Cauchy operator C are in fact 
null-solutions of D. 
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Theorem 4. Retain the same hypotheses as in Theorem 2 and assume that V = 0 
nearLl. Ifu £ C^ft,^) is such that AT (u) £ L 2 (<9ft), then 

ff(Du) £ L 2 (<9ft) & u = C(u\on) in ft Du — 0 in ft. (6.5) 

Proof. The right-to-left implication is a consequence of Cauchy’s reproducing 
formula (5.7). As for the opposite implication, first notice that u — C{u\qq) => 
Lu = 0 in ft. With this at hand, (2.4) gives 

[ \Du\ 2 dV — f ( Du,Du)dV — f ( d y u,u)ds , (6.6) 

Jq Jq Jon 

since V = 0 near ft. Thus, in order to conclude that Du = 0 in ft, it suffices to 
show that d v u = 0 on <9ft. 

With this goal in mind, note that our Green type identity (5.8), in concert with 
the assumption u = C(u\dn), forces S(d u u) = 0 in Cl. Taking the boundary 
trace entails S(d„u) = 0 on dQ. Assuming for a moment that S is injective on 
L 2 (dD,£), we may finally conclude that d v u — 0, as intended. Thus, the proof 
is finished, modulo showing that 

S:L 2 (dn,£)^L 2 (dD,£) (6.7) 



is one-to-one. 

To see this last claim, let / G L 2 (dQ , £) be such that Sf = 0 on dfl, and set 
u := Sf in £l±. It follows that u £ Hq i 2 (Q±,£) and Lu = 0 in D±. Thus, u = 0 
in by virtue of Proposition 2. Let 9 £ T*M be transversal to dft. Then, 

0 = (Vfltt) |an + - (Vou)|af 2 _ = -i(0, v)a(L~ l \ v)f ( 6 . 8 ) 

pointwise a.e. on dft. Since ( 0 , v) 7 ^ 0 on dCt and since cr(L -1 ; v) is an isomor- 
phism, the identity ( 6 . 8 ) forces / = 0. This justifies the claim made about (6.7) 
and finishes the proof of the theorem. □ 



6.7 Hardy spaces in Lipschitz domains 

For a first-order differential operator D as in (2.1) and for 0 < p < 00, define the 
Hardy space associated with D in a domain ftcMas 

D) := {u : ft £ ; Du - 0 in ft, and AAi G L p (dD)}. (7.1) 

The modern theory of Hardy spaces of holomorphic functions in Lipschitz sub- 
domains of the complex plane originates in [4, 10]. Here we continue this line of 
work by proving the following. 

Theorem 5. Let ft be a Lipschitz subdomain of M and assume that the first- 
order differential operator D satisfies (4.1), (4.3). Then, for each u £ TL p (£l, D), 
1 < p < 00, the following hold. 
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(1) u has anon-tangential limit u\dna.e. and u\dn G L p (d£l,£). Furthermore, 

IManll Lp(dQ) ~ \\Nu\\ LP(dn)- 

(2) u\qq G L\(dUt,£) if and only if N (V u) G L p (d£l). Moreover, there holds 
IMddl L\{dVL,E) ~ \\N(Vu)\\LP(dn)- 

(3) u\qq G B p,q (d£l,£) if and only if u G £), with equivalence of 

norms. 

(4) u\oq G L 2 (dfl,£) if and only if u G L 2 ^ 2 (0,£) if and only if 

|Vu|dist(-,0ft)5 G L 2 (0); 

(5) u | G L\(dQ,£) if and only if u G L 2 / 2 (0,£) if and only if 

|V 2 u|dist(.,0O)^ G L 2 (0); 

(6) Under the additional assumption (5.9), the Plemelj-Calderon decomposi- 
tion L p (d£L,£) = H p (£l+,D)\dn ® H p (£l-, D)\qq holds ; here the direct 
sum is both algebraic and topological. 

In the case when M — C and D = , (4)-(5) go back to Kenig’s area-function 

estimate in [10]; see also [4]. Regularity results in the same spirit for M = C m 
and D = d have also recently been considered in [16]. 

Proof of Theorem 5. The idea is to choose V = 0 near 0, and then use the fact 
that any u G H P (Q,D) is the Cauchy integral extension of its boundary trace, 
i.e., u = C(u\qq) in 0. At this stage, we may rely on Theorems 2-3 in order 
to prove the claims (l)-(5). As for (6), the fact that / = C f\dn+ + Cf |ao_ for 
any / G L p (d£l,£) shows that H p (£l±,D ) |an generate L p (d£l,£). To see that 
the sum is direct, assume that u± G H p (£l±, D) are such that u + |an = ^-Ian- 
Then, for every x G Q, u+( x ) = C(u+\dn){x) = C(u_\on){x) — 0, integrating 
by parts. Thus, u± = 0, and the desired conclusion follows. □ 

For (m - l)/m < p < 1, similar properties hold provided one employs (lo- 
cal) atomic Hardy spaces f)£ t (9Jl) (see [11] for definitions) in place of L p . We 
now proceed to discuss a compensated compactness type result, extending work 
from [18]. 

Proposition 4. Let 0 < p, q < oo and (m — l)/m < r < oo be such that 
1 + 1 = 1. Then 

p q r 

(u|ao,a(DV)w|an) £ 
for every u G H P (V, D ) and v G D l ). 



(7.2) 
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Proof. This can be proved by mimicking the approach in [26]. There, the con- 
text is that of (generalized) Cauchy-Riemann systems in the upperhalf-space but 
since the approach utilizes only nontangential maximal function estimates and 
cancelations based on integrations by parts, it can be extended to the present set- 
ting. □ 



6.8 Clifford analysis for functions of several variables 



In this section we consider in more detail the last example in §3, corresponding 
to Clifford separately monogenic functions. For a fixed integer n, let (R m ) n be 
embedded in A := ©C£(R m ), the sum of n copies of C¥(R m ). That is, if x E 
(R m ) n = R mn we write x = (x u • • . ,x n ) with Xj = Yl=i x< j e <* € C£(R m ) 
for 1 < j < n. In this setting, the Dirac operator corresponding to the j - th copy 
of C£(R m ) is 



Dj = XI e a 

Q!=l 




( 8 . 1 ) 



while the global Dirac operator reads B := {Dj)j. It follows that B is real (i.e., 
6 = 0) and -0*0 = ^ . A j — A. Here A j is the Laplacian in the j - th factor 
of the Cartesian product (R m ) n , for each 1 < j < n, and A is the Laplacian in 
the whole space R mn . Call a C^(R m )-valued function u defined in a region Q of 
R mn separately monogenic if DjU = 0 in ft for each j = 1, . . . , n. Note that this 
amounts to the requirement that Oi£ = 0; in particular, a separately monogenic 
function is harmonic. 

Recall the standard fundamental solution for the Laplacian in (R m ) n , i.e., 



. J-5F lo gM. if mn = 2; 

E(x) = < 

{ E (^) = - (mn-L. > assuming mn> 3; 

hereafter, u> mn stands for the area of the unit sphere in R mn . 

Starting from this, we set 

E(x, y ) := ^2 E ( x _ y) ei ® 

/ 

where e/ := if I = (u, . . . , ^), and the sum is performed over strictly 

increasing multi-indices. This corresponds precisely to the setup discussed at the 
beginning of §5, adapted to the current context. 

Going further, fix a bounded Lipschitz domain Q C R mn and denote by 

v = : dQ -4 (R m ) n A 

the outward unit normal to dQ (with Vj = J2 ^“e a ). In particular, 

dy = y] Vj * Dj , 
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where each product is to be understood in the Clifford algebra sense. One can then 
easily check that, with (•, •) denoting the natural (bilinear) pairing in C£( R m ), we 
have (d l/y E(x,y)J(y)) = T(x,y) ■ f(y) where 



T(x, y ) := D yj E(x , y ) ■ 

i = i 



LOr\ 



E X J 
It — 



Xj - y 3 



3 = 1 



y\ r 



Vj(y), ( 8 . 2 ) 



can be thought of as the natural Cauchy kernel in this setting. Then, in complete 
analogy to (5.3), the corresponding Cauchy integral operator is defined by 

Cf(x) = f r (x,y)'f(y)ds y , xeR mn \dO, (8.3) 

Jdfi 

for / : dd — > C£{R m ). Let us point out that while T(x, y) is harmonic in x, it is 
not separately monogenic when n > 1. Thus, the Cauchy operator is separately 
monogenic if and only if n = 1. The fact that C reproduces separately mono- 
genic functions (see Theorem 1) can then be viewed as a version of the classical 
Bochner-Martinelli formula in C m ; we leave the easy details of this verification to 
the interested reader. At this stage, having completed the parallel with our earlier, 
general theory, we may conclude that 



all main results in §6-§7 hold in the current setting . (8.4) 



Other, more specific results, are dealt with separately below. 

A distinguished feature of the setting we are currently considering is the fact 
that a Hartogs type theorem is valid in this context. Specifically, we have the 
following. 

Theorem 6. Let ft be an open set in R mn , n > 1, and let K be a compact 
subset of ft with Cl \ K connected. Then each separately monogenic function 
F : Tt\ K — > C£(W n ) extends to a separately monogenic function in Q. 

The proof is presented in Appendix A. This result has many remarkable corollar- 
ies. Here we only want to single out a few. 

Corollary 1. Assume that n > 1 and fix f) C M mn a bounded Lipschitz domain 
with a connected complement. Then 

u is separately monogenic in R mn \ Cl » ^ 

, . , , j . - . u = 0. (8.5) 

and u is bounded near infinity 



Proof. This is a direct consequence of Theorem 6 and Liouville’s theorem (for 
harmonic functions). □ 



We now discuss conditions guaranteeing that a function / defined on dfl ex- 
tends to a separately monogenic function F in Cl. 

Corollary 2. Let fl C R mn , n > 1, be a bounded Lipschitz domain with a 
connected complement. Then for each f £ L\(dCl,C£( R m )), the following are 
equivalent: 
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(i) There exists a (unique) separately monogenic function F in Q such that 
A f(F) G L 2 (dQ) and for which F\oq = /; 

(ii) Same as above, except that also N (V F) G L 2 (d£t)\ 

(iii) Cf = i/; 

Ov) C 2 f = \f. 

Proof That (i) and (ii) are equivalent follows from (2) in Theorem 5. The impli- 
cation (i) =4> (iii) is a direct consequence of (5.7) and the point (3) in Theorem 2. 
Also, clearly, (iii) => (iv). There remains to show that (iv) => (i). To this end, 
let F± := Cf in fi±, so J\f(S7F±) G L 2 (dQ). Next, so we claim, 

F± = ±C(F±\en) in (8.6) 

Indeed, by uniqueness in the Dirichlet problem for harmonic functions in Dahl- 
berg’s sense (see [11]), it suffices to show that both sides in (8.6) have identical 
traces on dQ. In turn, this is a straightforward consequence of definitions, the 
point (3) in Theorem 2, and our working assumption, (iv). 

Granted (8.6), it follows now from Theorem 4 that are separately 

monogenic in and respectively. In particular, Corollary 1 forces F_ = 0 
in Consequently, / = F+\on — F-\oq = so (i) is satisfied by choos- 
ing F := F + . □ 

It is important to point out that the above corollaries above fail for n = 1. 
The main difference, in this latter case, is that (iv) in Corollary 2 holds for any 
f G L 2 (dQ, C^(R m )); see [19]. 

In closing, let us mention that several other variants are possible. For instance, 
a completely analogous theory, indeed more akin to the classical theory of several 
complex variables , can be developed in connection with the Cauchy-Riemann 
operator (3.5). 



6.9 Octonionic analysis 

Let us first discuss a general algebraic construction. Assume that A = ( A , +, •) 
is a real, unitary, associative algebra endowed with a linear involution OHa c . 
The latter is called a conjugation if (ab) c = b c a c for every a,b e A. In particular, 
if 1^ is the multiplicative unit in A , then 1 C A = 1^. For later reference, define 
the trace of a £ A as Tr^(a) := ^(a + a c ). Next, consider the functor A h-> 
K(A) = (A x A, +, •) where the addition in K(A) is componentwise while the 
multiplication and conjugation are, respectively, 

(a, b) • (a, /3) := (aa — (3 c b , (3a + 6a c ), (a, b) c := (a c , —b). (9.1) 

The following claims are straightforward (see also [5]). 
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Proposition 5 . Assume that the algebra A is as above. Then: 

( 1 ) K(A) is a real, unitary algebra, equipped with a conjugation. 

( 2 ) (1.4,0) is the multiplicative unit in K(A). Also, if(e 3 ) J( z j are imaginary 
units which anticommute in A, then (0,1^), (e^, 0), (0,e ; ), j £ J, are 
imaginary units which anticommute in K(A). 

( 3 ) If A is commutative, then K(A) is associative. 

( 4 ) If A is associative and if any element has a scalar trace, or, Tr^ (a) £ R 
for any a £ A, then K(A) is an alternative algebra, i.e., 

[x,y, z \:= {xy)z -x(yz) ( 9 . 2 ) 

is trilinear alternate, and any element in K (A) has a scalar trace. 

( 5 ) IfTv^(ab) = Tr ^(ba)for any a,b £ A, then we have 

T *K(A){xy) = T *K(A){y x ) 

foranyx,y £ K(A). 

( 6 ) If A is associative and Tr^(afe) = Tr^(ba) for any a, b £ A, then 

T *K(A)(( x y) z ) = r ^K(A){y{ zx )) 
or, alternatively, [x,y,z] = 0 for any x,y,z £ K(A). 

( 7 ) If A is a normed algebra and aa c = |a| 2 , then K(A ) becomes a tiormed 
algebra with \x \ 2 = xx c = x c x = \a \ 2 + \b \ 2 for every x — (a, b) £ K(A). 
In particular, the multiplicative inverse of any element x £ K(A) \ 0 is 
x~ l — x c /\x\ 2 . 

( 8 ) Assume that A is a normed algebra such that aa c — |a| 2 , and so that any 
element in A has a scalar trace. Then K ( A ) turns into a real Hilbert space 
with respect to the pairing 

(x,y) :=Tr K(A){xy c ) = \{xy c + yx c ). ( 9 . 3 ) 

( 9 ) K(R) = C, the complex numbers, K(C ) = HI, the quaternions, and 
K( H) = O, the Cayley algebra of octonions. 

Specializing these to the case of the Cayley algebra of octonions gives: 

Proposition 6. The following are true: 

( 1 ) Ifi,j,k are the standard anticommuting imaginary units in H, then £1 : = 

MU2 := (j,0),6 := (k, 0),& := (0,1), & := := (OJ), 

£7 (0, k) are anticommuting imaginary units which, along with £o •— 

( 1 , 0 ), the multiplicative unit in O .form a basis for O. Accordingly, one 
can embed R 8 <— > O by identifying each point x = (x a ) a £ R 8 with the 
Cayley number ^ x a £ a £ O. 
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(2) Equipped with the pairing (9.3), the Cayley algebra O becomes a real, 
eight-dimensional Hilbert space, with {£ a } a an orthonormal basis. 

(3) Ifx = ^2 x a£a € vv/f/i x a G R, then x c — xq — — • • • — x-j^j. In 

particular, Tro(x) = xq. 

(4) x(xy) = x 2 y, (xy)y = xy 2 and {xy)y~ l = x = y^ 1 {yx). 

The first-order differential operator we intend to study in this section is 

< 9 - 4 > 

a = 0 

which acts on O- valued functions (defined in R 8 ) in a natural fashion. Set D c := 
2£od Xo — D. An incisive observation is recorded below. 

Proposition 7. The adjoint of the operator ( 9.4) with respect to the pairing ( 9.3) 

is —do + £i<9i H b ^ 7 ^ 7 . That is, D * = —D c . In particular, —D*D = A, the 

Laplacian in R 8 , so D is of Dirac type. 



Proof. The crucial observation is that 



(x • y,z) = {y,x c • z), Vx,y,zG O. (9.5) 

Indeed, by (9.3), the above is equivalent to TYo ((xy)z c ) — Tro (y(z c x)) which 
holds true, by virtue of ( 6 ) in Proposition 5. With this at hand, it is then straight* 
forward to finish the proof of the proposition. □ 

Once the Dirac character of D in (9.4) has been established, all our general 
results developed in §4-§7 will work when specialized to the present context as 
well. 

Thus, one can talk about analysis in several Cayley variables in (R 8 ) n ; i.e., 
the function theory for the differential operator D := ( Dj)j , where Dj is the 
octonionic Dirac operator (9.4) in the j- th factor in (R 8 ) n , for 1 < j < n. Here 
we only want to point out that, in this latter context, the natural Cauchy operator 
in a Lipschitz domain fl C (R 8 ) n , n > 1 , is 

c f(x) = f- f 3k ‘{ ,/ M-f(y)) ds v> x ? d n- ( 9 - 6 ) 

^8 n JdQ~i \ x ~y I 

Above, x = (xi,...,x n ) G (R 8 ) n , and v : dCt — > (R 8 ) n , v = (yf)j, is the 
outward unit normal to dCl. That (9.6) is a particular case of (5.3) can be seen 
much as we have done for (8.3) in the Clifford analysis context. 
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6.10 Integral representation formulas 

Here we work out Pompeiu’s integral representation formula in some important 
special situations. First, let D be as in (3.3) and set u := (/, 0), with / an arbitrary 
differential form. In this context, the first line in (5.6) gives 

f{x) = i f \(cr{d*\v y )d y E(x,y)J(y)) + (a(d-,i/y)d*E{x,y)J{y))\ds y 
JdCi 1 > 

~ [ (d y E(x,y),df(y))dV y + f (d*E(x,y),d*f(y))dV y , (10.1) 
Jn Jn 

for each x G 0, where E(x,y) is the Schwartz kernel of (A - V)~ l . If we 
integrate by parts in the last solid integrand, the identity (10.1) becomes 

f{x)=i [ ( a(d*]iSy)dyE(x,y)J(y))ds y (10.2) 

Jdn 

- J {{d y E(x, y),df(y)) + (d y d* y E(x,y), f{y))}dVy, x £ O. 

Next, the identity d*( A - V) = (A - V)d* - dV V • translates, at the level of 
Schwartz kernels, into d*E(x,y) = d x E(x,y) + R{x,y) where the (residual) 
double form R(x , y) is 0( \x - y |~( m_4 )). In fact, R = 0 when V is a constant. 
Accordingly, (10.2) further becomes 



f{x)= (d y E(x, y), [v A f )(y)) ds y 
Jon 



I (d y E(x,y),df(y)) dV y 
n 



-d 



X 



I (d y E(x,y)J(y))dVy 
n 



4 - a residual term. 



(10.3) 



The residual term (a weakly singular integral operator acting on /) disappears 
when V is constant. E.g., when certain Betti numbers vanish for the manifold M 
(so that A -1 exists), or when M — R m , we can take V — 0. 

The same program, carried out above, works equally well in the case when M 
is a complex manifold, and D is the complex conjugate of the operator (3.4). In 
this setting, the analogue of (10.3) reads 



fi x ) 



I (d y E(x,y),(d’°Af)(y))ds y - (d y E(x,y),(df)(y))dV y 
on J n 




I (d y E(x,y)J(y))dV y 
n 



a residual term. 



(10.4) 



Much as before, the residual term (in general, a weakly singular integral operator) 
is not present when V is constant. 

This is a version of the Bochner-Martinelli-Koppelman formula (see [21] for 
the case when M is C m ). Martinelli’s original version was proved in 1938 for the 




